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ABSTRACT 

Dynamic  positioning  of  a  moving  platform  is  analyzed  using  data  from  the 
Global  Positioning  System  (GPS)  acquired  in  Phase  II  of  the  Seafloor  Benchmark 
Experiment  on  R/V  Point  Sur  in  August  1986.  GPS  position  determinations  are 
compared  to  simultaneous  Mini- Ranger  fixes. 

The  GPS  positions  computed  using  only  broadcast  ephemeris  data  were  within 
20  m  from  the  Mini-Ranger  fixes  when  data  from  four  satellites  were  used  and  within 
30  m  when  data  from  three  satellites  and  a  geoidal  height  constraint  were  used.  It  was 
found  that  the  position  accuracy  is  degraded  when  data  from  a  satellite  reaching 
culmination  is  used. 


THESIS  DISCLAIMER 


0 


The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may 
not  have  been  tested  for  all  cases  of  interest.  While  every  effort  has  been  made,  within 
the  time  available,  to  ensure  that  the  programs  are  free  of  computational  and  logical 
errors,  they  cannot  be  considered  validated.  Any  application  of  these  programs 
without  additional  verification  is  at  the  risk  of  the  user. 
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I.  INTRODUCTION 


The  Seafloor  Benchmark  Experiment,  a  project  of  the  Hydrographic  Sciences 
Group  of  the  Oceanography  Department  at  the  Naval  Postgraduate  School  (NPS),  was 
initiated  in  1984  with  a  goal  to  establish  stations  on  the  seafloor  in  real-time  with 
geodetic  accuracy  (Saxena,  1982).  The  first  experiment  (Phase  I),  conducted  in  May 

1985,  included  a  configuration  of  four  benchmarks  on  the  seafloor  (Spielvogel  et  al., 

1987).  Further  study  indicated  a  sizeable  improvement  in  the  positional  accuracy  of  the 
inner  stations  as  compared  to  the  accuracy  of  the  outer  stations  of  a  larger  net 
configuration  (Kumar  and  Saxena,  1985).  For  this  reason,  during  the  Phase  II 
experiment  in  August  1986,  a  nine-station  net  configuration  was  used  (Kumar,  1986). 

Tlnee  types  of  GPS  receivers  and  a  Mini-Ranger  (MR)  Falcon  system  were  used  to 
determine  the  position  of  the  R,'V  Point  Sur. 

Keeping  the  goal  of  establishing  and  continuing  GPS  research  at  NPS,  Brown 
(1986)  established  the  Texas  Instruments  4100  GEOSTAR  GPS  software,  written  by 
the  Naval  Surface  Weapons  Center  for  their  CDC  Cyber  865  computer,  on  the  NPS 
IBM-3033  computer.  Brown  was  able  to  validate  the  modified  software  using  static 
position  data,  but  not  the  KALMN2  program  which  uses  dynamic  position  data. 

The  main  objective  of  this  thesis  is  to  determine  the  positions  of  a  dynamic 
platform  using  GPS  and  to  compare  them  with  MR  Falcon-determined  positions. 

Other  objectives  were  to  install  software  for  processing  the  raw  data  from  cassettes  to 
FIC  (Floating- Integer-Character)  format,  modify  the  KALMN2  program  and  validate  j 

it  using  Phase  II  data.  Programs  to  plot  and  analyze  the  continuous  ship  positions  by  i 

GPS  and  by  MR  Falcon  were  also  written. 

The  data  used  in  this  thesis  were  acquired  during  Phase  II  of  the  Experiment. 

The  GPS  data  are  limited  to  those  collected  by  the  TI-4100  receiver  from  NSWC  i 

installed  on  board  of  the  R/V  Point  Sur.  The  day  chosen  for  analysis  was  16  August 

1986,  since  adequate  pitch  and  roll  data  are  available  for  that  day. 
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II.  INSTRUMENTATION 


I 


A.  POSITIONING  SYSTEMS 

1.  GPS  receivers 

During  the  experiment,  the  following  GPS  receivers  belonging  to  several 
agencies  and  companies  were  installed  on  the  R/V  Point  Sur: 

-  TI-4100  receiver  from  Naval  Surface  Weapons  Centers  (NSWC) 

-  TI-4100  receiver  from  National  Ocean  Survey  (NOS) 

-  Eagle  Mini-Ranger  receiver  from  Motorola,  Inc. 

-  Trimble  4000A  from  Trimble  Navigation 

There  were  shore  based  GPS  receivers  at  the  following  stations: 

Beach  Lab 

-  TI-4100  from  Pacific  Missile  Test  Center  (PMTC) 

-  Trimble  4000A  from  Trimble  Navigation 
Ferrier 

•  Eagle  Mini- Ranger  from  Motorola  Inc. 

Dome 

-  Magnavox  Manpack  from  GPS  Joint  Program  Office. 

2.  Other  positioning  System 

For  navigation  and  positioning  on  board  were  a  MR  Falcon  system  from 
Motorola,  Inc.,  and  a  Loran-C  receiver. 

B.  OTHER  SYSTEMS 

An  acoustic  positioning  system  from  Oceano  Instruments  to  acquire  acoustic 
ranges  (Kuo,  1985)  was  installed  on  the  R/V  Point  Sur;  it  also  gives  the  position  of  the 
ship  in  relation  to  the  bottom  benchmarks  (acoustic  transponders)  and  can  be 
integrated  with  a  satellite  receiver  for  computation  of  their  positions  using  the  GPS 
coordinates  of  the  ship's  acoustic  transducer. 

In  order  to  allow  for  corrections  due  to  the  pitch  and  roll  of  the  ship,  there  was  a 
data  acquisition  system,  based  on  an  HP9826  computer,  to  collect  pitch,  roll  and 
heading  data  every  second.  Due  to  interfacing  problems  there  was  no  possibility  of 
logging  heading  data.  The  ship’s  Data  Acquisition  System  (SDAS)  collected,  among 
other  data,  the  ship's  heading  with  a  sampling  period  of  about  19  s. 
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C.  ANTENNAS  AND  TRANSDUCER  POSITIONS 

The  antenna  Gf  the  MR  Falcon  was  mounted  on  the  mast  of  the  ship  in  order  to 
increase  its  range.  The  antenna  of  the  TI-4100  from  NOS  was  installed  close  to  the 
MR  Falcon  antenna.  All  the  other  GPS  antennas  were  mounted  on  an  elevated 
wooden  table  above  the  ship's  laboratories  and  close  to  the  transducer  position  (Figure 
2.1). 


Figure  2.1  Relative  positions  of  antennas  and  transducer. 
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III.  DATA  PROCESSING 


A.  COORDINATE  SYSTEM 
1.  Parameters 

In  order  to  compare  the  positional  accuracies  of  the  TI-4100  GPS  receiver 
with  the  Mini-Ranger  (MR)  Falcon,  computations  were  made  in  a  single  coordinate 
system. 

Although  there  was  a  Doppler  survey  to  find  the  WGS  72  coordinates  of  all 
control  stations  used,  due  to  the  impossibility  of  setting  a  station  over  Dome  Ecc,  the 
site  for  one  of  the  MR  Falcon  shore  stations,  the  WGS  72  coordinates  of  nearby 
station  Dome,  68  m  away,  were  obtained.  During  the  Phase  II  experiment,  GPS 
provided  coordinates  in  the  WGS  72  system  (Bomford,  1980,  Sepplin,  1974),  although 
it  converted  to  WGS  84  on  1  January  1987  (Decker,  1986). 

The  coordinates  of  Dome  were  computed  in  NAD  27  and  then  transferred  to 
Dome  Ecc.  The  coordinates  of  Dome  Ecc  were  converted  from  NAD  27  to  WGS  72  in 
order  to  compute  directly  the  ship  positions  in  this  system,  avoiding  the  problem  of 
converting  each  position  individually. 

The  subroutines  used  for  converting  geodetic  coordinates  from  WGS  72  to 
Universal  Transverse  Mercator  (UTM)  and  vice-versa  wer  written  in  Fortran  and  are 
based  on  Basic  subroutines  used  by  DMA  for  an  HP9826  desktop  computer.  A  cross 
check  was  made  which  revealed  a  maximum  conversion  error  of  0.01  m,  which  is  well 
within  the  accuracy  limits  of  the  system  and,  hence,  does  not  affect  the  analysis  of  the 
data  processed. 

The  following  parameters  were  used: 

WGS72 

Semi  major  axis :  6,378,135  m 

Flattening  :  1  '298.26 

UTM  projection 

Central  meridian:  123°  00'  00”  W 

If  this  package  of  subroutines  is  used  for  any  location  other  than  the  area  of  the 
Seafloor  Benchmark  Phase  II.  the  UTM  parameters  must  be  changed  in  order  to  get 
the  correct  conversions. 
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Figure  3.1  Area  of  the  experiment. 

2.  Coordinates  of  stations 

The  two  shore  stations  were  located  on  high  points  along  the  coast  (Figure 
3.1).  Their  locations  and  coordinates  in  NVGS  72  are  as  follows: 
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Ferrier 

Latitude  e  36°  33'  53". 748  N 
Longitude- 121°  53'  52".939  W 
Dome  Ecc 

Latitude  36°  18'  20".897  N 
Longitude  121°  54'  00". 594  W 

The  elevations  of  the  different  antennas  were  set  in  the  Mini-Ranger  receiver  and  the 
measured  slant  distances  were  converted 

to  appropriate  horizontal  distances  by  the  receiver  system  during  data  acquisition. 

B.  REDUCTIONS  OF  POSITIONS  TO  THE  TRANSDUCER 
1.  Offsets  of  transducer  relative  to  antennas 

Since  the  two  antennas,  namely  the  MR  Falcon  antenna  near  the  top  of  the 
mast  and  the  TI-4100  antenna  on  the  elevated  table  above  the  lab  on  the  R/V  Point 
Sur,  were  at  different  locations,  a  common  point  had  to  be  chosen  to  compare  both 
systems.  The  acoustic  transducer  mounted  beneath  the  hull  was  chosen  as  a  common 
point,  which  was  used  for  further  integration  with  the  acoustic  data,  to  compute  the 
position  of  the  Seafloor  Benchmarks. 

For  data  reduction  there  was  a  need  for  pitch,  roll  and  ship's  heading  data. 
While  the  existing  pitch  and  roll  data  were  sampled  with  a  period  of  1  s,  the  ship's 
heading  was  sampled  with  a  period  of  about  19  s.  Computed  headings  were  tried  to 
improve  the  rate  but  tests  indicated  that  the  real  data,  even  with  a  low  sampling  rate 
was  better  than  the  computed.  The  method  used  to  compute  the  heading  is  described  in 
Section  C.3. 

The  geometric  offsets  between  the  MR  Falcon,  TI-4100  antenna  and  the 
acoustic  underwater  transducer  were  computed  using  a  right-handed  coordinate  system 
centered  in  each  antenna,  and  having  the  Y  axis  oriented  towards  the  bow,  the  X  axis 
to  starboard  and  the  Z  axis  towards  the  zenith.  Using  this  convention,  all  offsets  of 
the  transducer  are  negative.  Using  the  formulas  developed  in  the  next  section,  the 
corrections  for  pitch,  roll  and  heading  data  were  computed  and  applied. 

For  the  MR  Falcon  positions  only  horizontal  corrections  were  applied,  since 
the  positions  are  computed  in  a  two-dimensional  system.  For  the  GPS  positions  the 
corrections  were  applied  to  all  three  coordinates. 
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For  both  types  of  positions  the  coordinates  are  converted  to  UTM,  corrected 
and  converted  baA  to  WGS  72.  For  the  GPS  positions  the  Z  correction  is  applied  to 
the  geoidal  height. 

The  following  offsets  were  measured  and  computed  by  Prof.  Tucker  and  Mr 
James  Cherry,  and  changed  to  the  above  format. 

Transducer  offsets  relative  to  the  GPS  antenna  : 

XOFF  =  -4.268  m 
YOFF  -  -2.955  m 
ZOFF  -  -9.505  m 

Transducer  offsets  relative  to  the  MR  FALCON  antenna  : 

XOFF  -  -5.345  m 
YOFF  -  -9.453  m 
ZOFF  -  -16.152  m 

The  course  was  computed  in  a  range  0°  to  360°,  and  the  pitch  and  roll  data 
were  stored  in  a  range  -90a  to  +■  90°;  following  this  convention  the  roll  is  positive  when 
the  port  is  up  and  the  pitch  positive  when  the  bow  is  down.  The  sign  of  roll  was 
changed  to  conform  with  the  mathematical  model  described  next. 

2.  Math  model 

The  method  to  compute  the  correction  to  the  coordinates  is  based  on  a  seven- 
parameter  transformation  described  in  MofTitt  and  Mikhail  (1980),  so  successive 
rotations  of  the  vector  defined  by  XOFF,  YOFF  and  ZOFF,  give  the  corrections  to  the 
coordinates  illustrated  in  Figure  3.2. 

The  rotation  matrices  are  defined  as  follows: 


cos(PITCH) 

-sin(PITCH) 


sin(PITCH) 

cos(PITCH) 
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Figure  3.2  Rotation  of  axis. 
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The  orientation  matrix  is  defined  by  the  following  equation: 

Mq  *»  Mp  Mr  Mh  (eqn  3.1) 

and  the  corrections  are  computed  by  multiplying  the  orientation  matrix  by  the  offset 
vector,  to  get  the  following  equations: 

DX  =  XOFF  cos(ROLL)  cos(HEADING)  4-  (eqn  3.2) 

YOFF  {cos( PITCH )  sin( HEADING)  + 
sin(PlTCH)  sin(ROLL)  cos( HEADING)}  + 

ZOFF  (sin(PITCH)  sin(HEADING)  - 
cosfPITCH)  sin(ROLL)  cos(IIEADlNG)} 


DY  -  XOFF  cos(ROLL)  sin( HEADING)  (-1.)  +  (eqn  3.3) 

YOFF  {tos(PITCH)  cos(HEADING)  - 
sin(PITCH)  sin(ROLL)  sin(HEADIN'G)}  + 

ZOFF  {sin(PITCH)  cos( HEADING)  + 
cos(PITCH)  sin(ROLL)  sin( HEADING)} 

DZ  -  XOFF  sin(ROLL)  -  (eqn  3.4) 

YOFF  sin(PITCH)  cos(ROLL)  + 

ZOFF  cos(PITCH)  cos(ROLL) 

C.  FALCON  MINI-RANGER  DATA  PROCESSING 

1.  Error  of  system 

The  MR  Falcon  was  used  during  the  cruise  for  navigation  and  positioning.  In 
post-processing  it  is  used  to  evaluate  GPS  data.  Since  the  Mini-Ranger  is  a  system 
frequently  tested  and  used,  its  error  is  known  to  be  ±  3  m  (Laurila,  1976:  Munson, 
1977). 

It  would  be  better  to  have  had  more  than  two  lines  of  position  (LOP)  to  allow 
a  minimization  of  the  error  in  the  computation  of  the  Mini-Ranger  position  and  also 
to  allow  a  smoothing  of  the  data  in  places  where  the  data  gets  noisy  as  was  seen  in 
post  processing.  In  Figure  3.3  it  is  possible  to  see  the  areas  where  the  data  is  noisy  and 
to  see  the  small  jumps  characteristic  of  the  system. 

2.  Antenna  positions 

For  computing  the  antenna  positions  the  coordinates  of  the  shore  stations 
were  converted  to  WGS  72  UTM,  and  then  each  position  was  computed  in  the  WGS 
72  L'TM  projection.  The  subroutine  used  to  compute  the  positions  was  picked  from 
program  UCOMPS,  a  utility  package  used  in  the  Hydrographic  Sciences  Group  of  the 
Oceanography  Department. 

3.  Heading  Calculation 

As  mentioned  previously,  the  ship's  heading  was  needed  to  compute  the 
transducer's  position,  but  there  was  no  data  available  with  the  1-s  periodicity  needed, 
so  headings  were  calculated  from  Mini-Ranger  data.  The  calculation  of  the  ship's 
headings  was  made  during  the  computation  of  the  two  antenna  positions  and  kept  on  a 
separate  file.  These  computed  headings  were  affected  by  the  pitch  and  roll  of  the  ship 
and  also  by  the  noisy  data  of  the  Mini- Ranger. 

IS 


Figure  3.3  Plot  of  ship's  track,  using  MR  Falcon  data. 
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Some  trials  were  made  using  different  running  averages  for  smoothing  the 
computed  heading^ in  order  to  get  a  course  close  to  reality.  It  was  found  that  the 
courses  obtained,  even  if  filtered  by  a  running  average,  oscillated  and  hence  were  not 
suitable  for  this  analysis. 

4.  Transducer  positions 

The  transducer  positions  were  computed  using  the  model  described  in 
Section  B.  The  pitch  and  roll  data  used  were  sampled  with  a  periodicity  of  1  s  most  of 
the  time  when  there  was  positioning  data.  The  ship's  heading  was  sampled  with  a 
periodicity  of  from  19  s  to  1  minute.  Although  tests  were  made  in  order  to  compute  the 
heading,  as  stated  in  the  previous  section  it  was  decided  to  use  recorded  data, 
constraining  the  analysis  to  periods  of  constant  course. 

5.  Programs.  Inputs  and  outputs 

During  the  development  of  the  programs  many  more  routines  were  written 
than  mentioned  here.  In  order  to  make  the  use  of  the  programs  easy,  the  total  number 
of  programs  to  process  the  MR  Falcon  data  was  reduced  to  four. 

A  unique  format  for  the  different  position  files  was  defined,  allowing  only  one 
plotting  program  in  which  the  user  must  change  the  title  according  to  the  type  of 
positions  plotted. 

Program:  FALCON 

Input:  Month,  day,  year,  hour,  minutes,  seconds  (date  time  tag), 

codel,  range  1,  strength  L  codel,  range2,  strength2 
Output:  Listing  with  date  time  tag  and  geographic  positions 

File  with  date  time  tag  and  geographic  positions 
File  with  date  time  tag  and  rough  course 
Source:  See  Appendix  A 

Program:  COURSE  SMOOTH 

Input:  File  with  date  time  tag  and  rough  course 

Output:  File  with  date  time  tag  and  smooth  course 

Listing  with  date  time  tag  and  smooth  course 
Source:  See  Appendix  B 
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Program:  TRANSDUC  FALCON 

Input:  eFile  date  time  tag  and  geographic  positions  of  antenna 

Output:  Tile  with  date  time  tag  and  geographic  positions  of  transducer 

Optional  listing  with  date  time  tag  and  geographic  positions 
of  transducer 

Source:  See  Appendix  C 

Program:  PLOT 

Input:  File  with  date  time  tag  and  geographic  positions 

Output:  Plot  with  track  of  the  ship  in  UTM  projection 

Source:  See  Appendix  D 

All  files  are  in  free  format,  and,  if  for  by  any  reason  a  fixed  format  is  needed,  the  user 
should  change  the  programs  to  get  the  output  desired.  The  MR  Falcon  data  were 
originally  on  floppy  disks  in  HP9816  format.  A  program  was  written  in  Basic  in  order 
to  convert  and  write  the  data  onto  a  9-track  magnetic  tape.  The  transfer  of  data  to 
mass  storage  on  the  \PS  IBM  mainframe  computer  and  data  manipulation  were  made 
using  existing  procedures  for  the  mainframe  (Favorite,  1986;  Mar,  1984). 

The  scale  and  the  area  of  the  plot  are  fixed  for  the  particular  area  of  the 
Benchmark  Experiment  but  may  be  changed  in  the  program;  thus  any  proportional  size 
can  be  obtained  by  changing  the  control  card  as  is  explained  in  the  manual  (W.R. 
Church  Computer  Center,  1981). 

As  the  plot  program  can  be  used  for  plotting  the  outputs  of  different 
programs,  the  title  should  be  changed  accordingly,  leaving  the  appropriate  lines 
uncommented.  An  hour-minute  tag  is  plotted  every  5  minutes. 

D.  GPS  DATA  PROCESSING 

1.  Reading  the  TI-41 00  cassettes 
a.  System  configuration 

The  data  from  the  TI-4100  receiver  is  recorded  on  digital  cassettes  in  a 
format  not  directly  compatible  with  the  existing  programs  at  NPS.  To  read  and  decode 
the  data,  it  was  necessary  to  purchase  a  MEMTEC  cassette  terminal  (Model  5450XL) 
and  obtain  the  appropriate  software  from  the  Applied  Research  Laboratories  (ARL)  of 
the  University  of  Texas  at  Austin.  The  MEMTEC  cassette  reader  was  connected  to  an 
IBM  PC  having  a  hard  disk  of  20  Mb  (Megabytes),  a  modem  and  a  math  co-processor. 


Preliminary  tests  made  to  verify  the  software  and  hardware  connections  of 
the  tape  deck  witlLthe  IBM  PC  indicated  good  results.  During  the  tests  it  was  found 
that  the  format  of  the  final  output  of  this  package  was  not  the  same  as  defined 
previously  by  Scott  and  Peters  (1983)  but  was  in  a  new  format,  FIC,  where  each  block 
is  organized  in  a  structured  way  with  control  records  followed  in  order  by  all  floating 
type  data,  integer  type  data  and  character  type  data.  The  first  program,  CON9TR,  a 
pan  of  the  library-  installed  by  Brown  (1986)  in  the  SEF  (Standard  Exchange  Format), 
was  replaced  by  a  new’  program  using  the  FIC  format. 

b.  Reading  cassettes  to  disk 

Each  cassette  is  dumped  to  disk  using  the  program  MFERD.  The  operation 
takes  about  20  minutes,  and  the  file  created  fits  on  one  floppy  disk  which  can  be  kept 
for  raw  data  backup. 

c.  Converting  the  data 

The  binary  image  of  the  cassette  is  then  converted  to  a  FIC  format  file 
using  the  program  GS2FIC,  which  takes  about  20  minutes.  This  format  is  suitable  for 
processing  in  the  IBM -PC  with  appropriate  software  to  be  developed. 

In  order  to  transfer  these  files  to  the  main  frame  or  other  computer,  it  is 
necessary  to  run  a  program  FICFICA  that  converts  the  FIC  files  to  an  ASCII,  format 
where  the  data  is  organized  in  SO-column  records.  This  operation  takes  about  30 
minutes,  and  the  file  created  by  this  program  can  be  dumped  to  a  magnetic  tape  or 
sent  to  a  host  computer. 

The  files  created  by  these  two  programs  cannot  be  stored  on  a  normal 
floppy  disk  and  were  erased  as  soon  as  the  data  was  on  the  mass  storage  of  the 
mainframe  and  tested. 

d.  Transfer  to  the  mainframe 

The  transfer  of  the  files  to  the  mainframe  was  done  using  a  micro-computer 
connected  via  modem  to  the  main  frame.  A  terminal  emulator  distributed  by  the  NPS 
Computer  Center,  SIMPC,  simplifies  the  transfer  of  files  (Simware  Inc,  1984). 

Each  file,  corresponding  to  one  cassette,  was  sent  to  a  disk.  This  operation 
took  about  3.5  hours,  and  sometimes  the  transfer  was  stopped  in  the  middle.  When  this 
happened  the  transfer  had  to  be  started  from  the  beginning. 

Each  file  has  around  15000  lines  and  occupies  about  65%  of  an  A-disk. 
Using  an  extra  disk  every  file  w'as  converted  to  a  job  format  in  order  to  transfer  it  to 
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mass  storage.  The  data  were  sent  to  a  member  of  a  partition  data  set  where  ail 
cassettes  were  stored. 

2.  Conversion  of  formats 

After  a  study  of  the  FIC  format,  the  input  of  KALMN2  program,  and  the 
program  CON9TR,  a  program  was  written  called  CVFICA,  which  had,  after  some 
debuging  and  improvements,  a  final  form  that  fulfilled  the  needs  for  processing  the 
data  for  this  thesis. 

CVFICA  decodes  each  FICA  block,  producing  two  files  and  a  listing  of 
warnings  and  general  information.  One  file  contains  the  data  to  be  processed  by  the 
KALMN2  program  and  the  other  has  the  positions  computed  by  the  GESAR  software 
installed  in  the  receiver. 

In  the  present  version  of  the  program  the  pseudo-ranges  and  other  tracking 
information  are  stored  only  after  the  navigational  data  for  the  desired  number  of 
satellites  is  in  the  output  file,  since  this  data  is  needed  to  compute  the  positions  in  the 
Navigational  Mode.  It  is  very  easy  to  change  this  program  to  process  the  data  using 
the  precise  ephemeris.  The  program  flags  as  bad  the  data  that  belong  to  a  Space 
Vehicle  (SV)  whose  data  does  not  correspond  to  the  existing  SV  identification  in  that 
tracer.  The  data  are  also  flagged  as  bad  when  the  status  vector  has  any  value  other 
than  zero. 

The  program  is  able  to  convert  the  following  data  block  types: 

101  -  GESAR  Versions  1.0+  Input  data, 

3  -  GESAR  solutions, 

6  -  Tracking  data, 

S  -  Tracking  configuration, 

109  -  Navigation  Message  record  as  transmited, 

9  -  De-blocked  subframes  1-3  from  block  109 
1 1  -  Receiver  error  block, 

13  -  Tape  header/ trailor. 


According  to  the  CVFICA  program  any  other  block  is  read  and  dumped  to 
the  listing  with  a  warning  saying  that  the  program  was  not  able  to  handle  it.  No  such 
warning  was  found  after  using  the  program  with  data  from  10  cassettes. 

The  listing  also  indicates  when  the  navigational  data  for  each  satellite  is 
received  in  order  to  eive  to  the  user  the  status  of  the  constellation  of  satellites. 
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As  this  program  outputs  the  GESAR  position  solutions  and  the  data  to  the 
KALMN2  program  into  data  set  files,  these  must  be  created  prior  to  running  the 
program. 

A  control  file  is  used  to  pass  to  the  program  information  such  as  two  lines  for 
the  title  of  the  listing  and  the  starting  date  of  the  GPS  week.  The  starting  date  allows 
the  conversion  of  the  time  tag  in  seconds  to  a  date  time  tag.  The  program  does  not 
work  with  data  that  crosses  the  one-week  limit. 


3.  Program  information 
Program:  CVFICA 

Input:  FICA  files  without  the  first  two  comment  lines 

Control  file  with  the  title  and  starting  date  of  GPS  week 
(Appendix  J) 

Output:  File  with  GESAR  solutions 

File  in  NSWC  format  (input  to  KALMN  program) 


Listing  with  warnings  and  general  information 
Source:  See  Appendix  E 

4.  Computation  of  positions 
a.  Problems  found 

When  staning  to  use  the  Kalman  filter  program  (Brown,  1986)  for 
processing  the  GPS  data  some  difficulties  were  encountered,  either  due  to  the 
replacement  of  the  CON9TR  program  by  an  equivalent  one  or  by  errors  existing  in  the 
program  KALMN2.  The  first  problems  were  caused  by  the  different  units  used  in  the 
ARL  and  NSWC  software  which  were  a  consequence  of  the  lack  of  documentation 
during  the  development  of  the  program  CVFICA;  these  problems  were  easily  found 
and  solved.  The  last  problems  took  a  long  time  to  find,  not  only  due  to  the  complexity 
of  the  program,  but  also  due  to  the  nonexistence  of  processed  data  to  check  the  results. 
A  big  effort  was  made  in  order  to  correct  all  errors,  but  only  when  two  cassettes  were 
sent  to  NSWC  and  processed  there  were  all  errors  removed  from  the  processing.  It  was 
found  that  the  last  two  major  errors  were  due  to  causes  external  to  the  KALMN2 
program,  one  an  error  in  the  ARL  software  that  affected  the  ambient  temperature  by  a 
factor  of  10  and  caused  an  excessive  tropospheric  correction  and  the  other  a  bad 
constant  in  the  control  file  of  the  KALMN2  program. 
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b.  Changes  made 

In  addition  to  the  corrections,  some  changes  were  made  to  improve  the 
outputs.  Now  the  program  prints  a  date  time  tag  just  prior  to  all  outputs  labeled  with 
time  tag  in  seconds  of  GPS  week.  A  title  is  printed  as  a  header  of  the  output  listing. 
Other  small  changes  were  made  to  improve  the  code.  All  changes  were  made  carefully, 
and  the  programs  were  throughly  tested. 

c.  Data  processing 

After  inserting  data  from  ten  TI-4100  cassettes  into  the  mass  storage,  a 
control  file  for  the  KALMN2  program  was  made  (Brown,  1986).  For  an  initial 
approximate  position,  the  position  of  the  receiver  based  on  the  GESAR  solutions  with 
the  time  tag  closest  to  the  starting  time  was  used  for  the  data  processing.  For  the 
receiver  time  bias,  an  iterative  procedure  was  used;  thus,  the  data  were  processed  using 
a  value,  starting  with  zero,  that  was  replaced  by  one  listed  in  the  outputs  of  the 
KALMN2  program  when  the  computed  positions  were  close  to  the  predicted  ones. 

Due  to  the  good  initial  results  of  the  test  runs  with  certain  satellite 
configurations,  further  tests  were  investigated  for  the  different  types  of  solutions. 

The  value  of  -37  m  was  used  as  geoidal  height  in  order  to  constrain  the 
positions  to  the  geoid  when  using  data  from  less  than  four  satellites  or  when  forced  by 
the  control  file.  This  value  was  obtained  by  interpolation  on  a  chart  (Blaha  et  al.,  1986) 
and  corrected  for  the  antenna  height  above  sea  level,  giving  a  corrected  value  of -31  m. 
The  control  file  was  set  up  in  such  a  way  as  to  save  the  positions  on  another  file.  The 
broadcast  ephemeris  was  used,  as  the  precise  ephemeris  was  not  available  when  the 
data  were  processed. 

d.  Math  model 

To  compute  position,  the  program  uses  the  pseudo-ranges  corrected  for  the 
various  factors  affecting  them  and  an  eight-state  Kalman  filter  (Brown,  1986). 
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e.  Program  information 
Program :  eJCALMN2 

Input:  File  in  NSWC  format  produced  by  the  program  CVFICA 

Control  file  as  defined  by  Brown  (1986)(Appendix  I) 

Control  file  with  the  title  and  starting  date  of  GPS  week 
Output:  Listing  with  positions,  errors  and  other  information 

as  selected  by  the  control  file 
File  with  time  tag  in  seconds  of  GPS  week,  X,  Y,  Z 
and  Lat,  Lon,  HT,  both  in  WGS  72 
and  the  number  of  satellites  used  in  the  solution 
Optional  file  with  time  tag 

and  Geometric  Dilution  of  Precision  (GDOP)  every  minute 

5.  Transducer  positions. 

a.  Procedure 

The  positions  of  the  transducer  were  computed  using  an  algorithm  similar 
to  the  one  used  with  the  MR  Falcon  data.  The  only  difference  is  that  the  GPS 
positions  are  three-dimensional,  so  the  corrections  for  all  dimensions  were  computed 
and  applied. 

The  offsets  between  the  transducer  and  the  antenna  and  also  the  format  of 
the  input  files  were  different.  The  program  uses  the  same  control  file  as  the  program 
CVFICA  in  order  to  compute  the  date  time  tag  of  each  position. 

b.  Program  information 

Program:  TRANSDUC  GPS 

Input:  File  with  time  tag  in  seconds  and  WGS  72  coordinates 

Control  file  with  the  title  and  starting  date  of  GPS  week 
Output:  File  with  date,  time  tag,  latitude,  longitude  and  ellipsoid  height 

Optional  listing  with  above  information 
Source:  See  Appendix  F 
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It  is  possible  to  plot  the  track  of  the  ship  using  either  the  GESAR  solutions  or 
the  transducer  positions  as  described  before  in  Section  C.5.  Track  plotting  was  done  in 
order  to  give  the  GPS  data  coverage  and  data  status.  The  program  to  plot  is  the  same 
one  used  to  plot  the  positions  computed  from  the  MR  Falcon  data  but  with  the 
appropriate  title  (Figure  3.4). 


Figure  3.4  Plot  of  ship's  track  using  Tl-4100  data. 
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IV.  DATA  ANALYSIS 

A.  COMPUTATION  OF  DIFFERENCES 

1.  Approach 

To  find  the  accuracy  of  the  GPS  positions  in  a  two-dimensional  coordinate 
system,  the  observed  positions  in  WGS  72  coordinates  are  transformed  to  UTM 
coordinates  and  compared  with  the  UTM  coordinates  of  the  corresponding  MR  Falcon 
positions  interpolated  for  the  same  time  tag.  The  differences  in  the  coordinates  of 
Mini- Ranger  and  GPS  positions  are  expressed  in  meters. 

Due  to  the  low  sampling  rate  of  the  existing  course  data,  only  data  from 
steady,  straight  courses  were  used  for  analysis.  The  data  were  divided  into  eleven 
periods  of  constant  headings,  while  the  data  collected  during  the  turns  of  the  ship  were 
deleted.  This  procedure  made  it  possible  to  analyze  the  differences  between  GPS  and 
Mini-Ranger  data  during  the  different  courses. 

Due  to  unexpected  small  differences  found  for  positions  computed  using  data 
from  three  satellites,  several  tests  were  made  using  various  combinations  of  solutions. 
Although  for  data  analysis  many  combinations  were  studied,  three  representative  cases 
are  discussed. 

Different  starting  times  during  processing  of  the  GPS  data  using  the 
KALMN2  program  were  used  in  order  to  see  the  effects  of  the  propagation  of  errors 
due  to  the  initial  noisy  data,  but  the  discrepancies  were  inconsistent.  The  term  "noisy 
data"  is  used  to  refer  to  the  oscillatory  behavior  resulting  from  the  Kalman  filter  when 
the  number  of  satellites  changes  from  three  to  four  or  vice  versa.  The  effect  of  ship's 
motion  on  the  transfer  of  antenna  positions  to  two  chosen  common  points,  i.e.  the 
acoustic  transducer  and  the  GPS  antenna,  was  found  to  be  0.2  m,  which  falls  within 
the  accuracy  limits  of  GPS  and  hence  is  considered  negligible  for  this  study. 

2.  Case  studies 

The  following  cases  were  analyzed: 

Case  A  :  Data  from  four  satellites  were  used  or  data  from  only  three  available 

satellites  and  geoidal  height  constraint  for  position  computation  were 
used. 
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Case  B  : 


Same  as  case  A  but  using  the  geoidal  height  constraint  when  both 
three  or  four  satellites  were  used. 

Same  as  B  but  not  using  data  from  SV1 1. 


Case  C  : 

3.  Differences  obtained 

a.  Period  1.  19:14  to  19:20  ( Cases  A  and  B ) 


TABLE  1 

PERIOD  1.  CASES  A  AND  B. 

DIFFERENCES  BETWEEN  GPS  AND  MINI-RANGER  POSITIONS 


pos 


Svs 


Rangesfm) 

Avg  Min-Max 

A  ®B  A  B 


179 

6,8.11 

26 

26 

24-32 

24-32 

80 

6,8,9,11 

24 

30 

21-31 

29-31 

1 

6,8. ll 

27 

30 

27-27 

30-30 

* 

1 

6,8.9.11 

34 

42 

34-34 

42-42 

* 

1 

6.8.11 

34 

42 

34-34 

42-42 

* 

23 

6,8,9.11 

27 

31 

25-29 

30-33 

l 

6,8.11 

28 

31 

28-28 

31-31 

* 

74 

6,8,9.11 

23 

30 

20-30 

28-32 

*  noisy  data 

Case  A  -  using  four  satellites,  or 

four  satellites  with  geoidal  height  constraint 
Case  B  -  geoidal  height  constraint  in  case  A 

pos  *  number  of  consecutive  positions  computed  every  second 


During  the  first  part  of  this  period,  all  positions  were  computed  using  data 
from  three  satellites.  In  both  cases  (Table  1)  the  differences  found  are  in  the  same 
range,  24  to  32  m.  For  the  next  80  positions,  computed  using  data  from  four  satellites, 
Case  A  lead  to  values  lower  than  in  Case  B.  The  upper  limit  of  the  ranges  in  cases.  31 
m,  is  due  to  the  initial  oscillation  from  the  Kalman  filter. 

The  next  part  of  the  period  is  one  of  noisy  data,  i.e.  changes  from  four 
satellites  to  three  and  vice-versa.  In  both  cases  the  positional  accuracy  is  degraded, 
giving  larger  differences  between  GPS  and  Mini- Ranger  positions. 

At  the  end  of  the  period,  with  continuous  data  from  four  satellites  and  with 
the  exception  of  one  position,  the  differences  are  low,  reaching  20  m  in  Case  A  and  28 
m  in  Case  B. 


It  was  found  that  the  use  of  geoidal  height  as  a  constraint  for  position 
computations  wit£i  data  from  four  satellites  degrades  the  computed  position.  The 
reason  for  this  is-  explained  in  Section  4.B.  Noisy  data  causes  oscillations  in  the 
observed  differences  which  are  due  to  the  Kalman  filter. 


TABLE  2 

SATELLITE  VISIBILITY 
16  AUGUST  1986 


Time 

SV6 

El  Az 

SV8 

El 

Az 

SV9 

El  Az 

SV11 

El  Az 

SV13 

El  Az 

19:00 

59 

2 

51 

55 

24 

315 

57 

156 

- 

-- 

19:20 

63 

21 

44 

44 

32 

321 

67 

147 

6 

198 

19:40 

66 

45 

36 

38 

39 

328 

75 

124 

15 

200 

20:00 

66 

72 

28 

34 

48 

335 

77 

78 

24 

203 

20:20 

62 

96 

19 

33 

57 

342 

71 

74 

33 

208 

20:40 

56 

115 

11 

34 

67 

350 

62 

36 

43 

214 

The  offset  of  the  GPS  positions  in  relation  to  the  MR  Falcon  positions  is 
in  a  southeasterly  direction  (Figure  4.1  on  page  40).  Also  the  ship's  course  is  close  to 
the  azimuth  of  SV9  (Table  2),  which  must  be  the  reason  for  the  successive  tracking 
losses.  The  GPS  antenna  is  in  the  shadow'  of  the  ship's  mast  when  the  satellite  is  low. 
b.  Period  2.  19:20  to  19:25  (Cases  A  and  B) 

All  positions  were  computed  using  data  from  four  satellites.  There  is  a  nine- 
second  interruption  of  data  in  the  beginning  of  the  period  corresponding  to  the  time 
needed  to  change  tapes  in  the  receiver  during  logging. 

In  Case  A  the  differences  range  from  17  to  26  m  (Table  3).  The  maximum 
value  was  found  in  the  first  position  after  the  lack  of  data.  Most  of  the  difference*  are 
around  19  m.  In  Case  B  the  differences  range  from  23  to  38  m.  The  minimum  value  is 
found  after  the  gap  in  data.  Most  of  the  differences  are  around  29  m.  Visual 
comparison  is  possible  by  referring  to  Figures  4.2  and  4.3.  As  before  the  olTset  of  the 
GPS  positions  is  in  the  southeasterly  direction. 
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TABLE  3 


PERIOD  2.  CASES  A  AND  B. 

DIFFERENCES  BETWEEN  GPS  AND  MINT-RANGER  POSITIONS 

Ranges(m) 

Avg  Mm-Max 

pos  SVs  ABA  B 

295  6,8,9,11  19  29  17-26  23-38 

Case  A  -  using  four  satellites,  or 

four  satellites  with  geoidal  height  constraint 
Case  B  -  geoidal  height  constraint  in  Case  A 

pos  -  number  of  consecutive  positions  computed  every  second 


c.  Period  3.  19:25  to  19:30  (  Cases  A  and  B  ) 


TABLE  4 

PERIOD  3.  CASES  A  AND  B. 

DIFFERENCES  BETWEEN  GPS  AND  MINT-RANGER  POSITIONS 


Ranges(m) 

Avg  Mm-Max 

A  6B  A  B 


152  6,8,9,11  18  30  17-19  29-31 

1  6,8. 1 1  18  32  18-18  32-32 

l  6,8,9.11  13  28  13-13  28-28 

9  6,8, 1\  30  31  17-37  28-33 

118  6,8,9,11  26  33  24-37  32-36 

11  6,8,1 1  46  36  31-49  35-37 

4  6,8,9,11  39  37  36-43  37-38 

3  6,8,1 1  37  38  37-37  38-38 

*  noisy  data 

Case  A  -  using  four  satellites,  or 

three  satellites  with  geoidal  height  constraint 

Case  B  -  geoidal  height  constraint  in  Case  A 

pos  -  number  of  consecutive  positions  computed  every  second 


29-31 

32-32 

28-28 

28-33 

32-36 

35-37 

37- 38 

38- 38 


This  period  is  characterized  by  having  noise  in  different  parts  of  the  data 
and  shows  great  variation  in  the  differences.  In  Case  A,  at  the  beginning  of  the  period 
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the  average  differences  are  18  m,  after  the  first  period  of  noisy  data  they  reach  values 
of  26  m,  and  at  tlje  end  of  the  period  the  differences  are  up  to  46  m.  In  Case  B  at  the 
begining  of  the  period  the  average  differences  are  30  m.  This  values  changes  to  33  m 
after  the  noisy  data,  and  at  end  of  the  period  the  differences  are  up  to  38  m.  It  is 
possible  to  distinguish  the  oscillation  in  the  positions  due  to  noisy  data,  much  more 
emphasised  in  Case  A  than  in  Case  B  (Figures  4.4  and  4.5).  The  offset  of  the  GPS 
positions  is  as  before  in  the  southeasterly  direction. 

The  successive  periods  lacking  data  from  one  satellite  are  as  before  due  to 
SV9,  which  is  still  in  an  azimuth  close  to  the  ship's  course.  The  fact  that  use  of  geoidal 
height  as  a  constraint  with  data  from  four  satellites  for  computation  of  positions  docs 
not  improve  the  solution  is  evident  in  Table  4. 

d.  Period  4.  19:30  to  19:35  ( Cases  A  and  B) 


TABLE  5 

PERIOD  4,  CASES  A  AND  B. 

DIFFERENCES  BETWEEN  GPS  AND  MINT-RANGER  POSITIONS 


pos 

7 

40 

7 

1 

1 

49 

10 

130 

1 

11 

9 

2 
1 

81 


SVs 

6,8,11 

mi 

6.8.9.11 

6.8.11 

6.8.9.11 

6.8,  ll 

6.8.9.11 

6.8,  ll 

6.8.9.11 
6,8  ll 

6  8  9,11 


Ranges(m) 

Avg  Mm- Max 

ABA  B 


38 

38 

38-40 

38-39 

3 1 

35 

29-40 

34-39 

41 

38 

32-43 

37-39 

♦ 

38 

37 

38-38 

37-37 

* 

38 

38 

38-38 

38-38 

* 

30 

35 

27-38 

33-39 

40 

37 

31-44 

35-3S 

* 

31 

38 

28-41 

37-40 

34 

38 

34-34 

38-38 

* 

33 

39 

29-38 

36-45 

* 

41 

38 

33-46 

36-40 

* 

39 

37 

38-40 

37-38 

* 

38 

38 

38-38 

38-38 

* 

38 

41 

35-44 

40-45 

*  noisy  data 


Case  A  -  using  four  satellites,  or 

three  satellites  with  geoidal  height  constraint 
Case  B  -  geoidal  height  constraint  in  Case  A 

pos  -  number  of  consecutive  positions  computed  every  second 


This  period  is  characterized  by  successive  changes  in  the  number  of 
satellites  used  in  the  computation  of  positions;  the  differences  are  larger  than  observed 
in  the  previous  periods.Successive  losses  of  lock  on  SV9  were  the  cause  for  this.  In 
Case  A  (Figure  4.6)  it  is  possible  to  see  the  oscillations  of  the  computed  positions  when 
the  data  are  noisy’.  Improvement  in  the  solutions  is  seen  (Table  5)  when  we  change 
from  three  to  four  satellites;  the  average  differences  go  from  41  m  to  30  m,  becoming 
higher  values  when  only  three  satellites  are  available. 

In  Case  B  the  differences  have  higher  values  than  in  Case  A,  but  the  ranges 
are  smaller  due  to  the  smaller,  noisy  oscillations.  As  before  the  noisy  is  caused  by  the 
loss  of  lock  on  SV9,  and  the  offset  in  the  GPS  positions  is  in  the  southeasterly 
direction  from  the  Mini- Ranger  positions. 

e.  Period  5.  19:35  to  1940  ( Cases  A  and  B) 

This  period  has  the  noisiest  data  processed,  so  the  differences  have  greater 
spread.  Although  the  range  between  differences  is  wide,  20  to  68  m  in  Case  A  and  29 
to  67  m  in  Case  B,  the  differences  are  mostly  smaller  in  Case  A  than  in  Case  B. 

In  Figure  4.7  we  see  successive  jumps  from  three  satellites  to  four  and  vice 
versa.  The  GPS  positions  are  south  and  east  of  the  Mini-Ranger  positions.  The 
differences  in  the  northings  are  almost  double  of  the  ones  observed  in  the  eastings. 
The  reason  for  the  frequent  changes  in  the  number  of  satellites  used  is  the  relatively 
low  altitude  of  SV9  and  an  azimuth  close  to  ship's  heading,  which  together  with  the 
pitch  and  roll  led  frequently  to  shadowing  of  the  TI-4100  antenna. 

/.  Period  6.  19:47  to  19:53  ( Cases  A,  B  and  C) 

This  period  is  characterized  also  by  noisy  data,  but  it  is  for  a  different 
reason.  Here  data  from  SV8  is  rejected  during  the  phase  of  smoothing  it  with  Doppler 
data.  The  observed  roll  is  higher  than  in  previous  periods,  caused  perhaps  by  the 
almost  90°  change  in  course.  As  Table  2  shows,  SV8  is  setting,  having  a  lower 
elevation  than  any  other  satellite  used  during  this  period. 

There  are  no  significant  differences  between  Cases  A  and  B,  and  average 
differences  range  from  26  to  30  m.  In  Case  C  we  see  that  even  lowering  the  number  of 
satellites  to  two  or  three  and  using  the  geoidal  height  as  a  constraint,  smaller 
differences  are  obtained  than  the  ones  computed  with  data  from  S VI I  (Table  6).  The 
offset  in  the  positions  has  no  relation  to  the  course  of  the  ship,  i.e.  it  is  still  in  the 
southeasterly  direction,  and  a  sma’l  or  no  time  lag  exists  between  the  data  of  the  two 
systems  (Figure  4.8). 
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TABLE  6 


PERIOD  6,  CASES  A,  B  AND  C. 

DIFFERENCES  BETWEEN  GPS  AND  MINI-RANGER  POSITIONS 


Rangcs(m) 
A  BV8  C  A 


8.9.11.13 

8.9,  i  1.13 

9.11.15 

8.9,  il  ,13 

9:11,15 

8, 9  A  1.1 3 
9.tl,l5 

8. 9,  { 1.1 3 
9  Il,l5 

8.9,  f  1.13 
9  I 1,1 5 
8, 9  1.13 
9  1 1,1 5 

8. 9.1 1.1 3 

9.11.15 

8  9,1 1.13 

9  ll,l5 

8.9.11.13 
9  1 1,1 5 

8.9.11.13 

9.11.15 

8.9.11.13 
9  11.15 

8.9.11.13 


28  l9  l9 

^8  29  2^ 

27  28  26 

28  28  26 
28  28  26 
27  28  26 

27  28  26 

29  30  30 

28  29  28 

29  29  29 

28  29  27 

29  29  28 
28  29  27 
28  29  27 

28  30  27 

29  30  28 
27  29  27 
29  30  28 
27  28  25 
27  27  21 
27  28  24 
29  30  27 
26  27  23 
26  26  21 


28-29 

27- 27 

28- 29 
27-27 

27- 28 

28- 28 

26- 29 

27- 27 

27- 32 

28- 28 
28-30 
28-28 

28- 31 
28-28 

26- 31 
28-28 

29- 29 

27- 27 
29-29 
27-27 

26- 27 

27- 27 
29-30 
26-26 
26-28 


*  noisv  data 


Min- Max 
B 

28- 30 

29- 29 
29-30 
28-28 

27- 30 

28- 28 

26- 30 
58-28 

28- 33 

29- 29 

28- 30 

29- 29 

28- 31 

29- 29 

27- 31 

30- 30 
30-30 

29- 29 

30- 30 

28- 28 
26-27 
28-28 

m 

25-28 


28-30 

28-28 

28-29 

26-26 

25- 27 

26- 26 
24-28 
26-26 

27- 33 

28- 28 
27-30 
27-27 

26- 32 

27- 27 

24- 31 

27- 27 

28- 28 

27- 27 

28- 28 

25- 25 
20-23 
24-24 

13-23 

20-24 


Case  A  -  using  four  satellites,  or 

three  satellites  with  geoidal  height  constraint 
Case  B  -  geoidal  height  constraint  in  Case  A 
Case  C  -  ignoring  S V  J 1  but  constraining 
with  geoidal  height 

pos  -  number  of  consecutive  positions  computed  everv  second 


g.  Period  7.  19:57  to  20:05  ( Cases  A,  B  and  Q 

In  Case  B  there  is  an  improvement  in  the  positions  in  relation  to  Case  A.  If 
we  refer  to  previous  periods,  we  see  that  what  is  happening  is  the  reverse  of  what  was 
happening  there,  i.e.  the  use  of  the  geoidal  height  in  positions  computed  using  data 
from  four  satellites  improves  the  solution.  As  before,  we  see  that  noisy  data  causes  an 
oscillation  in  the  observed  differences  but  with  small  amplitudes.  Noteworthy  are  the 
highly  improved  positions  indicated  by  a  13-m  difference  (Table  7). 


TABLE  7 


PERIOD  7,  CASES  A,  B  AND  C 

DIFFERENCES  BETWEEN  GPS  AND  MINI-RANGER  POSITIONS 


pos 

SVs 

A 

Ranges(m) 

Avg 

B  C  A 

Min- Max 

B 

C 

185 

8,9.11.13 

25 

23 

21 

20-27 

19-25 

17-23 

1 

8  n.is 

24 

23 

20 

24-24 

23-23 

20-20 

* 

22 

8, 9  J  1.13 

22 

20 

18 

20-24 

18-22 

16-20 

I 

8  lU3 

22 

21 

18 

22-22 

21-21 

18-18 

* 

17 

8,9.  f 1.1 3 

22 

21 

18 

20-25 

19-23 

17-21 

1 

8  l3,l3 

23 

22 

19 

23-23 

22-22 

19-19 

4 

12 

8,9.  f 1.1 3 

22 

20 

17 

21-22 

19-21 

17-21 

1 

8  ll.l3 

19 

19 

16 

19-19 

19-19 

16-16 

4 

72 

8,9.11.13 

20 

18 

16 

18-22 

17-19 

13-18 

1 

8  ll.lS 

19 

19 

16 

19-19 

19-19 

16-16 

4 

13 

8, 9. f 1.1 3 

19 

18 

14 

17-20 

17-19 

13-16 

1 

8  ll.li 

18 

17 

14 

18-18 

17-17 

14-14 

4 

81 

8, 9. { 1,1 3 

21 

19 

16 

19-26 

17-24 

14-21 

1 

9  11.15 

17 

17 

14 

17-17 

17-17 

14-14 

4 

35 

8  9.11.13 

20 

18 

16 

19-22 

16-20 

13-18 

1 

9.1 1.1 3 

17 

16 

13 

17-17 

16-16 

13-13 

4 

2 

8,9,  f  1.13 

29 

25 

23 

28-29 

25-26 

22-23 

4 

1 

9,1 3.1 3 

27 

24 

21 

27-27 

24-24 

21-21 

4 

27 

8,9, f 1 , 1 3 

25 

23 

21 

18-36 

17-34 

15-31 

*  noisy  data 

Case  A  -  using  four  satellites,  or  .... 

three  satellites  with  geoidal  height  constraint 
Case  B  -  geoidal  height  constraint  in  Case  A 
Case  C  -  ignoring  Svil  but  constraining 
with  geoidal  height 

j  pos  -  number  of  consecutive  positions  computed  every  second 

I 

i  .  _  _  _  _  -  -  _ 

The  GPS  positions  are  also  south  and  east  of  the  MR  Falcon  positions,  but 
the  fact  of  ignoring  SVII  causes  a  decrease  in  the  differences  m  easting  (Figure  4.8). 
h.  Period  8.  20:07  to  20:13  ( Cases  A,  B  and  C) 

During  this  period  all  differences  in  all  cases  are  smaller  than  21  m  (Table 
8).  The  higher  values  are  obtained  after  the  computation  of  a  position  using  fewer 
than  four  satellites  in  Cases  A  and  B,  and  three  satellites  in  Case  C,  and  are  due  to  the 
oscillation  caused  by  the  Kalman  filter.  Table  8  shows  that  constraining  the  positions 
computed  using  data  from  four  satellites  improves  the  solution;  the  minimum  difference 
reaches  only  10  m.  Also,  ignoring  SV11  causes  a  reduction  of  the  differences,  mainly  in 
the  easting,  and  we  see  that  the  differences  are  smaller  than  in  the  other  cases  (Table 
8).  Figures  4.10  and  4.11  show  the  improvement  in  the  positions  when  SVII  is 
ignored. 
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TABLE  8 


PERIOD  8,  CASES  A,  B  AND  C. 

DIFFERENCES  BETWEEN  GPS  AND  MINT-RANGER  POSITIONS 


pos 

SVs 

A 

Range  s(m) 

Avg 

B  *C  A 

Min-Max 

B 

C 

37 

8,9.11.13 

19 

19 

14 

19-21 

16-18 

13-15 

1 

8.11.15 

18 

17 

14 

18-18 

17-17 

14-14  * 

153 

8,9.11.13 

20 

17 

15 

17-21 

16-19 

12-17 

1 

8A\a5 

17 

17 

13 

17-17 

17-17 

13-13  * 

28 

8, 9. f 1.13 

20 

17 

15 

19-21 

16-19 

13-16 

l 

8,11.1$ 

17 

17 

13 

17-17 

17-17 

13-13  * 

142 

S,9,  f 1 , 1 3 

19 

16 

12 

17-20 

15-17 

10-14 

*  noisy  data 

Case  A  -  using  four  satellites,  or 

three  satellites  with  geoidal  height  constraint 
Case  B  •  geoidal  height  constraint  in  Case  A 
Case  C  -  ignoring  S v  1 1  but  constraining 
with  geoidal  height 

pos  -  number  of  consecutive  positions  computed  every  second 


i.  Period  9.  20:13  to  20:20  ( Cases  A ,  B  and  Q 


TABLE  9 

PERIOD  9,  CASES  A,  B  AND  C. 

DIFFERENCES  BETWEEN  GPS  AND  MINI-RANGER  POSITIONS 


pos 

SVs 

A 

Avs 

B 

Ranges(m) 
;C  A 

Min- Max 

B 

C 

391 

8,9,11,13 

20 

17 

13 

19-22 

15-18 

11-14 

1 

no  data 

18 

15 

11 

18-18 

15-15 

11-11 

31 

8,9,11,13 

20 

15 

11 

18-21 

15-16 

10-11 

*  noisy  data 

Case  A  -  using  four  satellites,  or 

three  satellites  with  geoidal  height  constraint 
Case  B  -  geoidal  height  constraint  in  Case  A 
Case  C  -  ignoring  S V 1 1  but  constraining 
with  geoidal  height 

pos  -  number  of  consecutive  positions  computed  everv  second 


The  differences  found  during  this  period  are  similar  to  the  ones  obtained  in 
period  8.  Table  9  mdicates  that  the  differences  in  all  cases  are  correlated  as  before.  The 
oscillations  are  due  to  the  noise  in  the  data  from  the  MR  Falcon  and  not  to  the  GPS 
positions,  which  define  a  smooth  and  continuous  course  (Figure  4.12). 
j.  Period  10.  20:20  to  20:25  ( Cases  A,  B  and  C) 


During  the  first  tests,  the  differences  were  smaller  when  the  positions  where 
computed  from  three  satellites,  SV  1 1  excluded  and  the  ship  constrained  to  the  geoidal 
height,  than  the  ones  using  four  satellites,  SV11  or  SV6  included.  This  happened 
because  during  data  acquisition  SV6  was  selected  to  replace  SV11,  and  while  the 
receiver  was  not  locked  on  the  SV6,  only  data  collected  from  three  satellites  was  good 
(Figure  4.13). 

In  Case  B  (Figure  4.14  or  Table  10)  constraint  of  the  position  to  the 
geoidal  height  improves  the  solution  where  SV11  is  used.  In  Case  C  (Figure  4.15)  we 
notice  a  reduction  in  the  differences  in  the  positions  where  data  from  SV11  was  used. 
When  SV6  replaces  SV11,  in  all  cases  the  differences  are  larger  but  of  same 
approximate  magnitude. 


k.  Period  11.  20:30  to  20:37  ( Cases  A ,  B  and  C) 


TABLE  11 

PERIOD  11,  CASES  A,  B  AND  C. 

DIFFERENCES  BETWEEN  GPS  AND  MINT-RANGER  POSITIONS 


pos 

130 

287 


SVs 

6,8,9.13 

8,9: 1$ 


Ranges(m) 

Avg  Min-Max 

A  B  SC  A  B 


35  38  38 

36  36  36 


35-36 

30-37 


37-39 

30-38 


37-39 

30-38 


Case  A  -  using  four  satellites,  or 

three  satellites  with  geoidal  height  constraint 
Case  B  -  geoidal  height  constraint  in  Case  A 
Case  C  -  ignoring  SV11  but  constraining 
witn  geoidal  height 

pos  -  number  of  consecutive  positions  computed  every  second 


The  initial  positions  were  computed  using  data  from  four  satellites,  and 
after  20h32ml0.7s  the  receiver  started  to  track  SV12  in  place  of  SV8,  but  no  broadcast 
ephemeris  data  were  received  during  this  period.  The  dilFerences  computed  using  data 
from  three  satellites  are  the  same  in  all  cases,  but  those  computed  using  data  from  four 
satellites  have  different  values.  Thus  in  Case  A  the  differences  are  smaller  than  in  Cases 
B  and  C,  where  the  positions  are  constrained  to  the  geoidal  height  (Table  11).  In  all 
cases  the  differences  in  eastings  are  larger  than  in  northings,  contrary'  to  what  was 
happening  before;  this  offset  is  similar  for  the  other  two  cases  (Figure  4.16). 
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Figure  4.2  Period  2  (Case  A).  Using  all  available  satellites 
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TJt. 


Figure  4.6  Period  4  (Case  A).  Using  all  available  satellites. 


Figure  4.7  Period  5  (Case  A).  Using  all  available  satellites. 


SEAFLOOR  BENCHMARK  EXPERIMENT 
PHASE  I  I 
16  RUGUST  1986 


4  HR  Ffll  CON  rRRNOUCFR  CUSHIONS 
*  r,t*S  1 AANODCEA  POSH  IONS 


Satellites  used; 
SV8,  SV9,  SV13 


s1innTiTi?^rmrizi ;  rrr  1"  rzzn~zziiriirJ3r3tz~f 

HE1ER3 


Figure  4.8  Period  6  (Case  C).  Ignoring  SVl  1. 
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v.v 


>VV-'^%^'S*rv  -~s-~ 


Figure  4.10  Period  8  (Case  B).  Constraining  to  geoidal  height. 
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Figure  4.11  Period  8  (Case  C).  Ignoring  SVU 
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aVjI: 


SEAFLOOR  BENCHMARK  EXPER I  MEN! 
PHASE  IJ 
1$  AUGUST  1986 
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KGPS  TWWOUCE*  POSITIONS 
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Figure  4.12  Period  9  (Case  C).  Ignoring  SVi  1. 
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seafloor  benchmark  experiment 

PHASE  II 
16  AUGUST  1986 


lost  tracking  on  SV11 


-  start  tracking  on  SV6 


.20?  S 
\  2023 


2024 
^  2024 


Satellites  used: 

SV6,  SV8,  SV9,  SVU,  SV13 


S.202S 
\  2025 


+  «n  falcon  tmnoucfo  positions 
■t  CPS  TMMOUCeP  POSITIONS 


Figure  4.13  Period  10  (Case  A).  Using  all  available  satellites. 
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seafloor  benchmark  experiment 

PHASE  (I 
16  AUGUST  1986 


zoto  lost  tracking  on  SVll 


start  tracking  on  SV6 


2023 
\  2023 


2021* 

'k  2024 


Satellites  used: 

SV6,  SV8,  SV9,  SVll,  SV13 


2025 
\  2025 


+  MR  FW.CON  TRRNOUCEfl  POSITIONS 
*  GPS  TRUSOUCEP  POSITIONS 


Figure  4.14  Period  10  (Case  B).  Constraining  to  geoidal  height. 


seafloor  benchmark  experiment 
PHASE  II 
16  AUGUST  1966 


SV6 ,  SV8,  SV9,  SV13 
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4.  Programs.  Inputs  and  outputs 

Two  pro^ams  written  for  data  analysis  are  summarized  below: 

Program:  COMPARE  POSITION 

Input:  Control  file  with  period  limits  (Appendix  K) 

File  with  MR  Falcon  positions 
File  with  GPS  positions 

Output:  Listing  with  date,  time  tag,  and  GPS  geographic  positions 

and  differences  in  UTM  coordinates 

Source:  See  Appendix  G 

•  Program:  COMPARE  PLOT 

Input:  Control  file  with  period  limits  (Appendix  K) 

File  with  M  R  Falcon  positions 
File  with  GPS  positions 

Output:  Plots  of  tracks  using  GPS  and  MR  Falcon  data 

Source:  See  Appendix  H 

B.  CONCLUSIONS 

The  results  of  the  data  analysis  using  four  satellites  with  the  broadcast  ephemeris 
indicate  that  the  real-time  positioning  by  GPS  of  a  dynamic  platform,  e.g.  a  ship,  under 
the  best  conditions  is  within  15  m  but  in  most  cases  20  m. 

Using  data  from  three  satellites  and  constraining  the  solution  to  the  geoidal 
height,  positional  accuracies  are  within  20  m  under  the  best  conditions  and  30  m  in 
most  of  the  cases. 

It  wras  found  that  the  GPS  positions  were  always  southeasterly  of  the  Mini- 
Ranger  positions.  Due  to  lack  of  sufficient  data  and  due  to  time  limitations,  correlation 
between  the  geometry  of  the  satellites  and  ship  s  heading  could  not  be  investigated. 

The  accuracy  of  the  positions  is  correlated  with  the  relative  motion  of  the 
satellites;  this  was  found  when  SV11  was  ignored  as  in  Case  C,  which  improved  the 
solution.  If  the  variations  in  azimuth  and  elevation  of  all  satellites  are  compared  for  all 
periods  the  data  was  processed,  SVs  6  and  11  have  variations  larger  than  110  degrees, 
u'hile  SVs  8,  9  and  13  have  variations  in  azimuth  smaller  than  35  degrees  (Table  2).  It 


56 


can  be  seen  that  both  SVs  6  and  1 1  reach  their  culmination  during  the  period  the  data 
is  analyzed.  SV6  reaches  its  culmination  around  20:15  and  SV1 1  around  19:50. 

When  SVllls  approaching  culmination,  the  solution  is  not  improved  when  the 
position  is  constrained  to  the  geoid;  after  culmination  however,  application  of  the 
constraint  improves  the  solution  (Tables  12,  13  and  14). 

The  degradation  in  the  accuracy  of  the  GPS  positions  using  either  SV6  or  SV11 
is  evident  in  Table  15  where  the  differences  found  in  these  situations  are  bigger  than  for 
positions  computed  without  data  from  them.  It  is  evident  that  even  in  Case  C  when 
SV6,  with  the  same  characteristics  as  SVll,  is  used  to  compute  the  position,  the 
differences  become  larger. 

There  are  some  theoretical  studies  (Landau,  1986)  about  the  selection  of  satellite 
configuration  in  order  to  get  the  best  results  for  computed  positions.  However,  with 
the  present  data  it  is  not  possible  to  select  a  best  combination,  since  the  data 
acquisition  was  limited  to  a  small  observational  period  and  a  small  number  of  satellites. 
Currently  a  way  of  measuring  the  effect  of  the  geometry  of  the  satellite  configuration  is 
through  che  analysis  of  the  geometric  dilution  of  precision  (GDOP)  parameters.  These 
parameters  include  the  position  dilution  of  precision  (PDOP),  reflecting  the  dilution  of 
precision  in  three  dimensions,  the  the  horizontal  dilution  of  precision  (HDOP). 
reflecting  the  dilution  of  precision  in  two  dimensions,  the  the  vertical  dilution  of 
precision  (VDOP),  reflecting  the  dilution  of  precision  in  the  vertical  dimension,  and  the 
time  dilution  of  precision  (TDOP),  reflecting  the  dilution  of  precision  in  time  (Milliken 
and  Zoller,  1980).  Thus,  a  low  PDOP  provides  a  good  geometric  configuration.  Table 
16  lists  the  observed  satellites  and  their  GDOP  parameters. 

It  was  found  that  the  position  of  the  antenna  is  important,  and  in  the  case  of  the 
NSWC  TI-4100  receiver  some  error  was  caused  by  shadowing  due  to  the  ship's  mast 
when  the  azimuth  of  the  satellites  was  close  to  the  course  of  the  ship. 

The  positions  used  for  comparison  with  the  GPS  positions  are  themselves 
affected  by  the  noisy  signal  of  the  MR  Falcon  positioning  system.  This  problem  could 
be  reduced  if  more  than  two  LOPs  were  used  in  the  computation  of  positions. 
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TABLE  12 


MR  FALCON  POSITIONS  VERSUS  GPS  POSITIONS. 
OBSERVED  DIFFERENCES  IN  METERS  DURING  PERIOD  3 


Time 

Case  A 

Case  B 

Case  C 

H 

M 

S 

DX 

DY 

DP 

DX 

DY 

DP 

DX 

DY  DP 

19 

25 

0.7 

8 

-18 

19 

11 

-29 

31 

- 

— 

19 

25 

5.7 

8 

-17 

19 

12 

-27 

30 

19 

25 

10.7 

9 

-15 

18 

12 

-27 

29 

19 

25 

15.7 

8 

-17 

19 

12 

-28 

30 

19 

25 

20.7 

9 

-16 

18 

12 

-27 

30 

19 

25 

25.7 

9 

-16 

18 

12 

-27 

30 

19 

25 

30.7 

9 

-15 

18 

12 

-27 

29 

19 

25 

35.7 

9 

-16 

19 

12 

-27 

30 

19 

25 

40.7 

10 

-15 

18 

13 

-26 

29 

-- 

- 

19 

25 

45.7 

10 

-16 

18 

13 

-26 

30 

— 

DX,  DY  differences  in  the  UTM  coordinates 
DP  derived  position  difference  from  DX  and  DY 

Case  A  -  using  four  satellites,  or 

three  satellites  with  geoidal  height  constraint 
Case  B  -  geoidal  height  constraint  in  Case  A 
Case  C  -  no  data  during  this  period 
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TABLE  13 


MR  FALCON  POSITIONS  VERSUS  GPS  POSITIONS. 
OBSERVED  DIFFERENCES  IN  METERS  DURING  PERIOD  6 


H 

Time 

M 

S 

DX 

Case  A 

DY  DP 

DX 

Case  B 
DY 

DP 

DX 

Case  C 

DY  DP 

19 

52 

40.7 

17 

-22 

28 

19 

-22 

29 

16 

-21 

26 

19 

52 

41.7 

16 

-24 

28 

18 

-23 

29 

15 

-22 

27 

19 

52 

42.7 

16 

-24 

29 

18 

-23 

30 

16 

-22 

27 

19 

52 

43.7 

17 

-22 

28 

19 

-23 

30 

16 

-22 

27 

19 

52 

44.7 

15 

-24 

29 

16 

-25 

30 

14 

-24 

28 

19 

52 

45.7 

16 

-21 

27 

18 

-22 

29 

15 

-21 

26 

19 

52 

46.7 

10 

-27 

29 

12 

-28 

30 

8 

-27 

28 

19 

52 

47.7 

14 

-23 

27 

15 

-24 

28 

8 

-23 

25 

19 

52 

48.7 

14 

-23 

27 

15 

-23 

27 

9 

-20 

22 

19 

52 

49.7 

15 

-22 

26 

16 

-21 

26 

10 

-18 

21 

19 

52 

50.7 

17 

-20 

26 

18 

-19 

26 

12 

-16 

20 

19 

52 

51.7 

18 

-20 

27 

19 

-19 

27 

14 

-17 

21 

19 

52 

52.7 

19 

-19 

27 

20 

-19 

27 

15 

-17 

23 

19 

52 

53.7 

19 

-19 

27 

20 

-19 

28 

16 

-18 

24 

19 

52 

54.7 

13 

-27 

30 

13 

-27 

31 

10 

-26 

28 

19 

52 

55.7 

12 

-26 

29 

13 

-26 

29 

9 

-25 

26 

DX.  DY  differences  in  the  U'TM  coordinates 
DP  derived  position  difference  from  DX  and  DY 


Case  A  -  using  four  satellites,  or 

three  satellites  with  gecidal  height  constraint 
Case  B  -  geoidal  height  constraint  in  Case  A 
Case  C  -  ignoring  Svll  but  constraining  with 
geoidal  height 


TABLE  14 


MR  FALCON  POSITIONS  VERSUS  GPS  POSITIONS. 
OBSERVED  DIFFERENCES  IN  METERS  DURING  PERIOD  8 


H 

Time 

M 

S 

DX 

Case  A 

DY  DP 

DX 

Case  B 
DY 

DP 

DX 

Case  C 

DY  DP 

20 

7 

0.7 

11 

-15 

19 

9 

-14 

17 

5 

-13 

14 

20 

7 

20.7 

10 

-17 

20 

8 

-15 

17 

4 

-14 

14 

20 

7 

40.7 

10 

-15 

18 

9 

-14 

17 

4 

-12 

13 

20 

8 

0.7 

9 

-18 

20 

7 

-17 

18 

4 

-16 

16 

20 

8 

20.7 

10 

-17 

20 

8 

-16 

17 

5 

-14 

15 

20 

8 

40.7 

10 

-17 

20 

8 

-16 

17 

4 

-14 

15 

20 

9 

0.7 

11 

-17 

20 

8 

-15 

17 

5 

-14 

15 

20 

9 

20.7 

11 

-17 

20 

8 

-15 

17 

5 

-14 

14 

20 

9 

40.7 

11 

-17 

20 

8 

-15 

17 

5 

-14 

14 

20 

10 

0.7 

11 

-15 

19 

9 

-14 

16 

5 

-12 

13 

20 

10 

20.7 

9 

-17 

19 

8 

-16 

17 

5 

-14 

15 

20 

10 

40.7 

11 

-17 

17 

9 

-14 

16 

5 

-12 

12 

20 

11 

0.7 

11 

-16 

20 

8 

-15 

17 

4 

-13 

13 

20 

11 

20.7 

12 

-16 

20 

9 

-14 

17 

4 

-13 

13 

20 

11 

40.7 

12 

-16 

19 

8 

-14 

16 

4 

-12 

13 

20 

12 

0.7 

12 

-15 

19 

9 

-13 

16 

4 

-11 

12 

DX,  DY  differences  in  the  UTM  coordinates 
DP  derived  position  difference  from  DX  and  DY 

Case  A  -  using  four  satellites,  or 

three  satellites  with  geoidal  height  constraint 
Case  B  -  geoidal  height  constraint  in  Case  A 
Case  C  -  ignoring.  S V 1 1  but  constraining  with 
geoidal  height 
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TABLE  15 


MR  FALCON  POSITIONS  VERSUS  GPS  POSITIONS. 
OBSERVED  DIFFERENCES  IN  METERS  DURING  PERIOD  10 


Time  Case  A 

H  M  S  DX  DY  DP 


H 

M 

S 

DX 

DY 

DP 

20 

20 

7.7 

14 

-15 

20 

20 

20 

8.7 

14 

-15 

20 

20 

20 

9.7 

11 

-13 

17 

20 

20 

10.7 

9 

-13 

15 

20 

20 

11.7 

8 

-12 

14 

20 

20 

12.7 

7 

-12 

13 

20 

20 

44.7 

4 

-10 

10 

20 

20 

46.7 

4 

-10 

10 

20 

20 

47.7 

3 

-11 

11 

20 

21 

35.7 

3 

-10 

11 

20 

21 

36.7 

3 

-10 

11 

20 

21 

37.7 

8 

-13 

15 

20 

21 

38.7 

10 

-14 

18 

20 

21 

39.7 

12 

-15 

19 

20 

21 

40.7 

13 

-16 

21 

DX 

Case  B 
DY 

DP 

DX 

Case  C 
DY  DP 

9 

-13 

15 

4 

-10 

11 

9 

-13 

15 

4 

-10 

11 

8 

-12 

15 

4 

-10 

11 

7 

-13 

14 

3 

-11 

11 

7 

-12 

14 

3 

-10 

11 

6 

-12 

13 

j 

-10 

11 

4 

-10 

10 

4 

-10 

10 

4 

-10 

10 

4 

-10 

10 

3 

-11 

11 

3 

-11 

11 

3 

-10 

11 

3 

-10 

10 

3 

-10 

11 

3 

-10 

11 

8 

-13 

15 

8 

-13 

15 

10 

-14 

18 

10 

-14 

18 

12 

-15 

20 

12 

-15 

20 

13 

-16 

20 

13 

-16 

21 

DX,  DY  differences  in  the  UTM  coordinates 
DP  derived  position  difference  from  DX  and  DY 

Case  A  -  using  four  satellites,  or 

three  satellites  with  seoidal  height  constraint 
Case  B  •  geoidal  height  constraint  in  Case  A 
Case  C  -  ignoring,  Sv  1 1  but  constraining  with 
geoidal  height 


V.  RECOMMENDATIONS 


GPS  data  from  other  days  of  the  Phase  II  cruise  should  be  used  with  the 
modified  programs  discussed  here  to  further  check  our  conclusions. 

Tests  should  be  made  in  order  to  find  the  reason  for  offsets  in  a  southeasterly 
direction  and  evaluate  their  correlation,  if  any,  with  the  geometry  of  the  satellites  used. 

To  test  the  conclusion  that  the  accuracy  depends  on  the  large  variations  of  the 
azimuth  of  the  satellites.  SVs  other  than  SV6  and  S  VI 1  should  be  used. 

The  relation  between  satellites'  culminations  during  observations  and  the 
accuracy  of  the  positions  should  be  further  investigated  to  provide  a  criterion  for 
satellite  selection. 

A  better  antenna  site  on  the  ship  should  be  found  to  avoid  shadowing  and  to 
minimize  the  effects  of  pitch  and  roll.  Antennas  of  two  similar  systems  should  be 
installed  at  both  high  and  low  elevations  to  determine  the  effects  of  pitch  and  roll. 

The  Kalman  filter  should  be  improved  to  avoid  the  oscillation  when  data  are 
noisy.  Other  algorithms  should  be  implemented  and  their  results  compared  to  the 
results  of  the  KALMN2  program. 

To  provide  a  real-time  position  computation  the  KALMN2  program  or  a  similar 
program  should  be  installed  in  a  transportable  computer.  A  system  should  be  designed 
to  have  one  computer  converting  and  logging  the  data,  another  processing  and 
displaying  results,  saving  them  or  sending  them  to  still  another  computer  or  data 
recording  system,  and  there  should  be  a  master  system  to  control  the  synchronization 
of  time  recorded  with  all  data  and  to  control  the  flow  of  information. 

Processing  of  the  same  data  should  be  done  using  the  precise  ephemens  to  see 
how  the  solution  is  improved. 

Processing  in  differential  mode  should  be  done  to  see  how  much  the  solution  is 
improved. 

Comparison  of  data  from  other  receivers  should  be  done  to  check  which  type  ct 
receiver  provides  higher  accuracy. 

Next  an  experiment  should  be  designed  to  obtain  point  or  differential  positions  in 
real  time  for  moving  platforms  (survey  ship'  with  little  or  no  shore  support. 


To  allow  belter  evaluation  of  the  OPS  data,  positions  computed  w.th 
two  lines  of  position  should  be  used  as  reference. 
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APPENDIX  A 

PROGRAM  FALCON.  SOURCE  LISTING 


//EZEQUIEL  JOB  <  0812 ,9999 1 ,  •  EZEQUIEL'  ,CLASS*C 
//•MAIN  ORG=NPGVW1.0812P, LINES* <  99)  .CAROS* t  99) 

//•FORMAT  PR .DDNAME sGO . FT06F001 > 

//•FORMS*S6Pl 

//  EXEC  FORTVCLG 

//FORT .SYSIN  DD  * 

C 

C  AUTHOR :  AUGUSTO  EZE OUZEL 

C  DATE:  APRIL  10,  1987 

C 

C  DESCRIPTION: 

C  THE  PROGRAM  TAXES  THE  RAM  RANGE -RANGE  DATA  FROM  AN  INPUT  FILE, 

C  COMPUTS  THE  POSITIONS  AND  ANO  THE  COURSE  BETWEEN  TWO  ANTEMA 

C  POSITIONS  ANO  PRINTS  THE  UTM  COORDINATES,  THE  COMPUTED  COURSE 
C  ANO  THE  GPS  TIME  OF  EACH  OBSERVATION. 

C  THE  PROGRAM  RUNS  IN  MVS 

C  I/O  SPECIFICATIONS:  SEE  END  OF  THIS  JOB. 

C 

C  THIS  PROGRAM  IS  RESTRICT  AS  IS  TO  THE  SEAFLOOR  BENCHMARK  EXPERIMENT 

C  FOR  OTHER  USES  CO  NOT  FORGET  TO  ENTER  THE  COORDINATES  OF  THE 

C  STATIONS  ANO  OTHER  PARAMETERS. 

C 

C  COORDINATES  ARE  IN  MGS 72  DATUM. 

C 

C  THE  PROGRAM  HILL  RUN  WITH  ANY  AMOUNT  OF  DATA  ONLY  LIMITED  TO  DISK 
C  SPACE 
C 

C  ANY  BLANK  LINES  HILL  TERMINATE  THE  PROGRAM  IN  ERROR 

C 

C 

C 

C 

c 

REAL  *8  X<  2  )  >Y<  2  )  ,XPOS,YPOS,XOLD,YOLD 

INTEGER  LADEG1 , LAOEG2 , LAMIN1 , LAMIN2 , LODEG1 , LOOEG2 , LOMINI , LOMIN2 
INTEGER  MONTH .YEAR .DAY .HOUR ,MIN , LADEG . LAMIN , LODEG , LOMIN 
REAL  LASEC1 , LASEC2 , LOSEC1 , LOSEC2 , LASEG , LOSEG .SECS 
C 

C  STATION  1  FERRIER 

C 

LA0EG1*36 
LAMIN1=33 
LASEC1*53. 748 
LQ0ES1*121 
LOMINI *53 
10SEC1*52. 959 
C 

C  STATION  2  DOME  ECC 

C 

LA0E62*36 
LAMIN2*18 
LASEC2*20. 95 
LOOEG2*121 
L0MIN2*S4 
LOSEC2*00.50 
C 

C  COMPUTE  X-Y  FOR  CONTROL  STATIONS 
C 

CALL  GPUTMI LA0EG1 .LAMIN1 .LASEC1 .L0DEG1 .LOMINI .LOSEC1 ,Y(  1 )  ,X1 1  )  ) 
CALL  GPUTMf  LADEG2 .LAMIN2 .LASEC2 .LODEG2 . LOMIN2 , LOSEC2 , Y I  2  )  ,XI 2  )  ) 
C 
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C  ENO  OF  INITIAL  OATA 
C 

C 

C  COMPUTS  THE  FIRST  POSITION  TO  HAVE  A  STARTING  POINT  FOR 
C  EVALUATION  OF  THE  COURSE 

C 

READ! 50 >  •  ,£N0*100)  MONTH , DAY .YEAR .HOUR , MIN, SECS, 

#  KOOE 1 , RANG I . ISTRE 1 ,KOOE 2 .RANGE , I STRE 2 

CALL  RRXY (X,Y, RANG 1. RANGE, XOLD.YOLD) 

C 

C  READ  IN  OATA  FROM  INPUT  FILE 
C 

50  CONTINUE 

READ! 50,  «  , ENO *100  )  MONTH, DAY .YEAR. HOUR, MIN, SECS, 

»  K00E1.RANG1, ISTRE1 ,K00£2 , RANGE, ISTREE 
C 

C  COMPUTE  X-Y  AND  GP  OF  RR  POSITION 
C 

CALL  RRXYIX.Y ,RANG1, RANGE, XPOS, YPOS  ) 

C 

C  COMPUTES  THE  COURSE 
C 

CALL  HEAO! XOLO , YOLD ,XPOS , YPOS , COURSE , IERR ) 

C 

C  IF  LANE  JUMPS  DOES  NOT  SAVE  THE  POSITION 

C 

C  IF!  IERR.  LT.  0  1  WRITE!  6,1  HERR 

Cl  FORMAT!  16) 

IF! IERR. LT . 0  )  GOTO  50 
C 

C  COMPUTES  THE  GEOGRAPHIC  COORDINATES 
C 

CALL  UTMGPI  YPOS ,XPOS , LAOEG . LAMIN, LASEG , LOOEG , LOMIN , LOSEG ) 

C 

C  OUTPUT  THE  RESULTS 
C 

WRITE! 6.60  1  MONTH, DAY  .YEAR, HOUR, MIN, SECS, LADEG, LAMIN, LASEG 

#  , LOOEG, LOMIN, LOSEG, COURSE 

60  FORMAT! 2! IX, 12  ) , IX ,14 ,3! IX ,15 , IX ,12 ,1X ,F7 . 3  ) , IX , F5 . 0  ) 

C 

WRITE  I  7  )  MONTH , DAY .YEAR .HOUR , MIN , SECS , LAOEG , LAMIN , LASEG 

#  , LOOEG, LOMIN, LOSEG 
C 

C  SAVES  THE  COURSE  IN  A  SEPARATE  FILE 

C 

WRITE! 8)  MONTH, DAY, YEAR, HOUR .MIN, SECS, COURSE 
C 

C  SAVES  POSITION  FOR  COURSE  COMPUTATION 

C 

XOLO=XPOS 
YQLO*YPOS 
GOTO  50 
C 

C  ENO  OF  PROGRAM 

C 

100  CONTINUE 
STOP 
ENO 
C 

. . . 

SUBROUTINE  RRXY! XR ,YR ,R1 ,R2 .XCO.YCO ) 

C 

IMPLICIT  REALMS  (A-H.0-21 
DIMENSION  XRI  2  )  ,YRI  2  ) 

A  *  OSqRT! (XR! 2  )-XR< 1 1  l*»2  ♦  (  YRI  2 ) — YR  (  1  I )*»2 ) 

CR  *  I  XR!  2  )-XRI 1 ) )  /  A 
SR  *  (  YR!  2  )-YR(  1 )  )  /  A 
XP  »  (R1«R1  -  R2»R2  ♦  A»A 1  /  (2»A) 

ARG  *  R1*R1  -  XP*XP 
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YP  a  DSGRT 1 ARG  ) 

XCO  =  XP*CR  ♦  YP*SR  ♦  XR(  1 1 
YCO  *  XP*SR  -  YP*CR  +  YRI 1 ) 
RETURN 
ENO 


SUBROUTINE  GPUTMl LADEG .LAMIN.LASEG.LQOEG.LOMIN.LOSEG .NORTH .EAST 1 


DOUBLE  PRECISION  A ,R ,N,AP ,BP .CP ,DP ,EP ,S,R1 .ESG.ESGP .RH.RP ,KO 
COUBLE  PRECISION  R2 ,R3 ,R4,RS .P ,P2 ,P3 .P4.P5 ,P6 >A6 .BS.SINSEC 
DOUBLE  PRECISION  PHI > OLAN. NORTH .EAST .B.PHIMIN.PI .LON 
REAL  LASEG.LOSEG 


THIS  SUBROUTINE  COMPUTS  THE  UTM  COORDINATES  OF  SP  IN  WGS  72 
IN  ZONE  10  CENTRAL  MERIOIAN  125  00  00  H 


C 


C 

C 

c 


c 

c 

c 


c 


c 


c 


c 


c 


CMs-123.000 

PHI sOFLOAT ( LAOEG  )  +DFLCAT ( LANIN  )/60 . DO+OBLE I LASEG 1/5600 . 00 
LONaQFLOATl LOOEG  l+DFLOATI LOMIN  1/60 . 00+0BLEI LOSEG  1/3600. DO 
LON=-LON 

OLAIIst  LON-CM  1*3600 .  ODO 

A=6378135 . 000 
R-298. 2600 

K0=0 . 999600 

B  =  A*<  R-l . 00  1/R 

N  =  IA-B1/1A+B1 

A P  *  A*<  (  l.DO-N) *5.00/4. 00*< N**2-N**3  U81.  DO/64. 00*< N»*4-N**5 11 
BP  =  3. 00/2. 00*A*<  (N-N**2 1*7. 00/8. 00*1 N**3-N**4  1*55. 00/64. D0*N»*5  1 
CP  =  15. 00/16. 00*A*<N**2-N**3*S. 00/4. 00*(N**4-N**5> 1 
DP  a  35.00/48. DO*A*< N**3-N**4*ll . 00/16 . DQ*N**5 1 
EP  a  315 . 00/512 . 00*A*<  N**4-N**5  1 
PHIMIN  =  PHI*60. 00*2. 908882086660-4 

PI=OARCOS( -1.00  1 

PHI=PHI/180.D0*PI 

S  =  AP*PHIMIN-BP*OSIN(2.00*PHI  l«CP«0SIN( 4 . DO*PHI 1 
4  -  DP*DSIN(6.00*PHI  )+EP*0SIN(8.00»PHI  1 

R1  ~  KO*S 

SINSEC  *  ( 1.00/3600. DO  1/180. D0*PI 
SINSEC=OSIN(  SINSEC I 

ESQ  *  I  A**2-B**2  )/A**2 
ESOP  »  ESQ/I  1 . DO-ESQ  1 

RM  *  A*! 1 . DO-ESQ l/( OSQRT ( 1 . 00-ESQ*OSIN( PHI )**2 1 )**3 
RP  «  RH*( 1.00*ESQP*0COS< PHI  )**2  1 

R2  *  RP*OSIN( PHI )*OCOS( PHI  l*SINSEC**2/2 . D0*K0*1 . 008 
R3  *  SIN5EC**4*RP*DSIN( PHI  )*DCOS( PHI  1**3/24.00*1 5. DO-DTANI PHI  )**2 
4  ♦  9. DO*€SQP*0CCSI PHI  l**2*4. DO*ESGP*ESQP*OCOS( PHI  )**4  1*K0*1 .016 
R4  ■  RP*OCOSIPHI1*SINSEC*KO*1.04 

R5  *  SINSEC**3»RP*OCOS( PHI  l**3/6 .00*1 1. OO-OTANI PHI )**2 
4  ♦  ESQP*OCOSI PHI  )**2 >*KO*1.012 

P  *  . 0001DO*OLAM 
P2  a  P**2 
P3  a  p**5 
P4  a  P**4 
P5  a  P**5 
P6  a  p»»6 

A6  *  P6»SINSEC**6«RP*OSIN(PHI1»OCOSIPHI  )**5/720.00 


■  rV*  — V  -  S  v  1  h'  1' ra"T 16'  |J[ 


oooono  o  on  oooon  no  o  o 


unuMui  pjiu  u  * j  *jtrj»  *-V  *v*v  »vi  u  * 


*  *  (61. 00-58. 00*DTANIPHI)**2+0TAN(PHI)**4 

$  ♦  270.00*ES9P*DCOS<PHI  )**2-330 .00*ES9P*0SIN< PHI  >**2  )*K0*1.D24 

85  *  P5*SINSEC**5*RP*OCOS< PHI  )**5/120.D0*f 5.D0-18.D0*DTAN< PHI  )**2 
$  ♦  0TAN( PHI  )**4+14.DO*ES9P*OCOS< PHI  1**2 

*  -  58.00*ES9P*DSIN(PHI)**2)*KO*1.D20 

N0RTH*R1+R2*P2+R3»P4+A6 

EAST* I R4*P+R5*P3+B5 1+500000 . 00 

RETURN 

END 


C 


SUBROUTINE  HEAD! XOLO ,YOL0 ,XPOS , YPOS .COURSE . IERR  ) 


C 


REAL*8  XOLO, YOLO ,XPOS,YPOS 
IERR*0 

DX=XPOS-XOLD 
OY=YPOS-YOLO 
RA0IUS*SQRT( 0X**2+DY**2  ) 

IF( RADIUS. GT. 20.  )  IERR=-1 
IF( RADIUS. ST. 0.  )  GO  TO  10 
COURSE *0.0 
RETURN 
10  CONTINUE 

20  COURSE*  ARCOS1 DY/RADIUS  ) 

COURSE =C0URS£*180./S. 141592654 
IF(OX.LT.O-)  COURSE*360. -COURSE 
RETURN 
ENO 


SUBROUTINE  UTMGP ( NORTH , E AST , LAOEG , L AMIN , LASEG , LODEG , LOMIN , LOSES  > 


DOUBLE  PRECISION  A,R,N,B.  SP,CP,DP,EP,S,R1,ESO,ES9P,RM,RP,KO 
DOUBLE  PRECISION  R7.R8 ,R9 ,ES ,9,92 ,93 ,94 ,95 ,96 ,D6 , AP .SINSEC 
00U8LE  PRECISION  RIO ,DLAM,NORTH .EAST ,RPHI ,DPHI  .PHIMIN 
DOUBLE  PRECISION  EPRIME , DELTA, 06NUM, CM, PHI , LON, PI 
REAL  LASEG, LOSEG 

THIS  SUBROUTINE  COMPUTS  THE  UTM  COORDINATES  OF  GP  IN  WGS  72 
IN  ZONE  10  CENTRAL  MERIDIAN  125  00  00  H 

CM=-123 . 000 

A=6378135.0DO 
R=298. 2600 

K0=0 . 999600 

PI*OARCOS( -1.00) 

8  *  A*( R-l . DO  )/R 

N  *  (  A-B )/( A+8  ) 

AP  *  A*(  ( 1 .  DO-N  J+5.D0/4. 00*(  N**2-N**3  )+81 . 00/64. 00*( N**4-N**5  )  ) 

BP  *  3.D0/2 .D0*A*< ( N-N**2  f+7. 00/8. D0*t  N**3 -N**4  ) +55 . DO/64 . D0*N**5 ) 
CP  =  15. 00/16. 00*A*(N**2-N**3+3. DO/4, 00*1 N**4-N**5 ) ) 

OP  a  35. 00/48. D0*A*(N**3-N**4+11. DO/16. D0*N**5 > 

EP  *  315.00/512 . 00*A*l  N**4-N*«5 ) 

FIRST  APROXIMATION  OF  PHI 

PHI*NORTH/30 . 8D0/3600 . DO 

COMPUT  TRUE  MERIDIONAL  OIST  ANO  APROXIMATE  PHI 

DO  100  1*1,8 

PHIMIN*PHI*60 . 00*2 . 908882086660-4 
RPHI*PHI/180 . 00*PI 
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S  =  AP*PHIMIN-BP*OSIN(2.00«RPHX)+CP*OSINI4.00*RPHI) 
$  -OP*OSIN<  6.D0*RPHI )+EP*OSIN(8.00*RPHI ) 

R1  =  KO*S 
OELTA=NORTH-RI 
PHIX0ELTA/30. 8/3600. I+PHI 
100  CONTINUE 

RPHI=PHI/180.D0*PI 


SINSEC  =  ( 1.00/360000 1/180. 00*PI 
S1NSEC  =  DSINI  SINSEC) 

C 

ESQ  *  (  A**2-B**2  )/A**2 
ESQP  *  ES<3/(  1. DO-ESQ) 

RM  *  A*<1.  OO-ESQ  )/<DSQRT(1.00-ESQ*0SIN(RPHI)**2))**S 
C 

RP  =  RM*<1.0O+ESQP*OCOS(RPHI)**2) 

R7*< 0TAN1 RPHI  )/( 2 . D0*RP**2*SINSEC  ) )*( 1 . DO+ESQP»OCOS( RPHI I**  2 I 

*  *1.012/K0**2 

R8=< DTANI RPHI )/( 24 . D0*RP**4*SINSEC  ) )*( 5 . DO+3 . 00*0TAN<  RPHI )**2 

*  +6 . DO*ESQP*DCOS! RPHI )**2-6 . O0*ESGP*0SIN< RPHI  )**2-3 . D0*ESQP**2 

*  *DCOS(  RPHI  )**4-9.00*ESQP**2*0COS<  RPHI  >**2*0SIN( RPHI  1**2  ) 

*  *1.D24/K0**4 

R9=l . OO/OCOSl RPHI  )/( RP*SINSEC )»1 . D6/K0 

R10=l . OO/OCOSl RPHI  )/( 6 . D0*RP**3*SINScC  )*( 1 . 00*2 . 00*0TAN< RPHI  )**2 
6  +ESQP*DCOS(  RPHI )**2 )*1.018/K0**3 

C 

EPRIME=EAST -500000 . 00 
Q  s  ,00000100«EPRIHE 
Q2  *  <3**2 
C33  =  Q**S 
04  *  0**4 
Q5  a  Q**5 
06  =  0**6 

06NUM=Q6»OTAN( RPHI ) 

06  =  1  06NUM )/(  720 . 00*RP**6*SINSEC 1*1 61 . 00+90. DO*DTANI RPHI )**2 

*  +45 . DO*OTAN<  RPHI  1**4+107 . DO*ESOP*OCOSl RPHI  )**2 

S  -162 . 00*ESQP*0SIN(  RPHI )**2-45. DO*ESQP*OTANI RPHI  )»*2 

*  *OSINl RPHI  )**2  )*1 . D36/X0**6 

ES*< Q5*l . OO/OCOSl RPHI  ) )/( 120 . 00*RP**5*SINSEC 1*1 S . 00+28.00 

*  *OTAN( RPHI  )**2+24 . 00*0TAN( RPHI  )**4+6 . 00*ESQP*DCOS( RPHI )**2 

*  +8 . 00*ESQP*0SIN( RPHI  )**2  )*1 . 030/KO**5 
C 

DPHI=( -R7*Q2+R8*Q4-06 1/3600 . 00 
OLAM=( R9*Q-R10*Q3+ES  1/3600.00 
PHI=PHI+DPHI 
LON=CM+DLAM 

CALL  OMSIPHI.LAOEG.LAMIN.LASEG) 

CALL  OMSILON,LOOEG,LOMIN,LOSEG) 

RETURN 

ENO 

C 


SUBROUTINE  DMSI DEC , LDEG .MIN ,SEC  ) 


XNUH»OABS< DEC  1 
LOE6*OINT ( XNUH ) 

XMIN*I XNUM-OFLOATI LOEG ) 1*60 . 00 
MIN*OINT ( XMIN  I 

XNUM«IXMlN-OFLOAT( MIN)  1*60.00 
SEC=SNGLIXNLM) 

IF  I LOEG. GE. 360)  LOEG*L0EG-36O 

RETURN 

ENO 

/* 

//GO.FT06F001  DO  SYSOUT=* 


//GO. FT07F001  00  OSN*MSS .S0812. FALCON .ANTENNA . POS ,DISP*OLD 
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APPENDIX  B 

PROGRAM  COURSE  SMOOTH.  SOURCE  LISTING 


//40BSM0TH  JOB  I  0812,9999 ) , 1 EZEQUIEL ‘ .CLASSIC 
//•MAIN  0RG*NPGVH1 . 0812P , LINES* (  99 ) .CARDS*! 99  ) 

//•FORMAT  PR  »00NAME  *GO . FT06F001 , 

//*FORMS=S£Pl 

//  EXEC  FORTVCLG 

//FORT . SYSIN  00  • 

C 

C  AUTHOR:  AUGUSTO  EZEOUIEL 
C  DATE:  JANUARY  05,  1986 
C 

C  DESCRIPTION: 

C  THE  PROGRAM  MAKES  THE  RUNNING  AVERAGE  OF  THE  COURSE  COMPUTED 

C  BY  THE  PROGRAM  <  COMP  >  IN  ORDER  TO  TAKE  OUT  THE  EFECTS  OF  THE 

C  PITCH  ANO  ROLL  EFFECTS.  THE  FIRST  ANO  LAST  TWO  POSITIONS  ARE 
C  REJECTED  AS  NO  RUNNING  AVERAGE  IS  POSSIBLE  FOR  THEM. 

C  THE  PROGRAM  RUNS  IN  MVS 

C  I/O  SPECIFICATIONS:  SEE  ENO  OF  THIS  JOB. 

C 

C  THIS  PROGRAM  IS  RESTRICT  AS  IS  TO  THE  SEAFLOOR  BENCHMARK  EXPERIMENT 

C 

C  THE  PROGRAM  WILL  RUN  WITH  ANY  AMOUNT  OF  DATA  ONLY  LIMITED  TO  DISK 
C  SPACE 
C 

C  ANY  BLANK  LINES  WILL  TERMINATE  THE  PROGRAM  IN  ERROR 

C 

C 

C 

C 

C 

REAL  COURSE ( 5  )  >SECSl 5 ) 

INTEGER  MONTH!  S ) .YEAR! 5 )  ,DAYI 5 ) .HOURI 5  )  .MINI  5 )  .FLAG 
LOGICAL  SAVE 
C 

C  INITIALIZATION  OF  FLAGS 

C 

SAVE=  .FALSE. 

FLAG=1 

C 

C  READS  THE  FIRST  FIVE  COURSES 

C 

QO  100  4*1,5 

READI  51  >EN0*1100  )  MONTH!  J  )  ,OAY(J  )  .YEARU  )  .HOUR!  J  ) , 

*  MINI  J).SECS(J), COURSE!  J  I 

150  FORMAT! 21 IX, 12), IX, 14. IX, 12, IX, 12, IX, F4. I, IX, F07. 2) 

C  WRITE!  6 , 150  )  MONTHtJl.OAY!  J),  YEAR!  J), HOUR!  J), 

C  $  MIN!  4), SECS!  J), COURSE!  J) 

100  CONTINUE 
C 

C  AVERAGE  THE  COURSE 

C 

105  AVG*!  COURSE 1 1 ) ‘COURSE t 2  ) ‘COURSE ( 5  )  ‘COURSE ( 4  ) ‘COURSE I S  )  )/5 . 0 
C 

C  SAVES  THE  AVERAGE  IN  THE  MIOOLE  POSITION 
C 

COURSE! 5 )*AVG 
C 

C  DOES  NO  SAVE  THE  FIRST  TWO  POSITONS  OF  THE  FILE  AS  NO  RULING 
C  AVERAGE  EXISTS  FOR  THEM 

C 

IF  (SAVE)  GOTO  50 
C 


FLAG*FLAG+l 

I F ( FLAG. EG. 3  )  SAVE*. TRUE. 

GOTO  70 
C 

C  WRITES  THE  DATA  WITH  THE  AVERAGED  COURSE  INTO  THE  FILE 
C 

50  WRITE 1 7  I  MONTH! 1 )  »OAY(  1 ) , YEAR 1 1 )  .HOUR! 1 )  >MXN( 1 )  .SECS! 1 )  (COURSE I  1  ) 
C 

WRITE! 6 , ISO  1  MONTH! 1 1 ,OAY( 1 1 .YEAR! I )  .HOUR! I ) , 
t  MINUI.SECSU  »,  COURSE!  1) 

C 

C  MOVES  FORWARD  THE  DATA 
C 

70  CONTINUE 

00  75  J*2,5 
MONTH! J-I >*M0NTHI  J ) 

OAY! J-l  l*OAY( J  I 
YEAR!  J-1)*YEAR(  J) 

HOUR!  J-I  )*HOUR(J) 

MIN! J-I  )*MIN( J  ) 

SECSU-1  )=SECSU ) 

COURSE ( J - 1 )  ^COURSE ( J ) 

75  CONTINUE 
C 

C  READS  ONE  MORE  COURSE 

C 

READ! 51 >ENO*1000 )  MONTH! 5 ) .DAY! 5 ) . YEAR! 5  )  .HOUR! 5 ) . 

$  MINI 5), SECS! 5) .COURSE!  5) 

C 

C  GO  BACK  TO  COMPUT  A  NEW  AVERAGE 
GO  TO  105 

C 

C  ENO  OF  FILE  FOUND  DURING  THE  INPUT  OF  A  FITH  COURSE 
C  SO  SAVES  THE  REMAINING  TWO  AVERAGES 

1000  CONTINUE 

OO  1001  J*1,Z 

WRITE! 7,150!  MONTH! J ) ,DAY<  J I .YEAR! J  )  .HOUR! J ) , 

*  MIN ( J  I. SECS! 4  l, COURSE!  J) 

WRITE! 6,150)  MONTH! J  I  .DAY! J  )  .YEAR! J  )  .HOUR! J  I , 

<  MINI J  )  .SECS! J  I . COURSE U ) 

1001  CONTINUE 

STOP 

C 

C  NO  OUTPUT  IS  MADE  IF  AN  END  OF  FILE  HAS  FOUND  DURING  THE  READING 

C  OF  THE  FIRST  FIVE  COURSES 

1100  CONTINUE 

WRITE!  6,1401 

140  FORMAT  ( '  ERROR.  ENO  OF  FILE  OURING  THE  FIRST  SET  OF  COURSES' ) 
STOP 
ENO 

/* 

//SO. FT06F001  00  SYSOUT*» 

//SO.  FT07F001  00  OSN*MSS.SO«12. FALCON. COURSE. SMOOTH. DISP»SHR 
//GO . FT51F001  00  DSN*MSS.S0512. FALCON .COURSE . ROUGH ,DISP*SHR 

// 
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APPENDIX  C 

PROGRAM  TRANSDLC  FALCON.  SOURCE  LISTING 


//J08TRANS  JOB  ( 0812, 9999 h 'EZEQUIEL1 .CLASSIC 

//•MAIN  0RG*NPGVM1.0812P,LINES=( 99), CARDS*! 99) 

//•FORMAT  PR  >OONAM£*GO . FT06F001 • 

//•FORMS *SEP1 

//  EXEC  FORTVCLG 

//FORT.SYSIN  00  • 

C 

c 

C  AUTHOR:  AUGUSTO  EZEQUIEL 

C  DATE:  MARCH  04,  1987 

C 

C  DESCRIPTION: 

C 

C  THIS  PROGRAMS  TAKES  EACH  POSITION  OF  THE  ANTEMA 

C  THE  COMPUTED  COURSE  I  FILTERED  OR  NOT  )  .THE  PITCH 

C  ANO  ROLL  OATA,  ANO  COMPUTES  THE  POSITIONS  OF  THE 

C  TRANSDUCER,  APLYING  THE  SEVEN  PARAMETER  TRANSFORMATION 

C  TO  THE  OFFSETS  OF  THE  TRANSDUCER  IN  RELATION  TO  THE  ANTENNA 

C  CONSIDERED  THE  CENTER  OF  COORDINATE  SYSTEM  ANO  THEN  ADDING 

C  THE  CORRECTIONS  TO  THE  COORDINATES. 

C 

C  THE  PROGRAM  RUNS  IN  MVS 

C 

C  I/O  SPECIFICATIONS:  SEE  END  OF  THIS  JOB. 

C 

C  THIS  PROGRAM  IS  RESTRICT  AS  IS  TO  THE  SEAFLOOR  BENCHMARK 
C  EXPERIMENT 

C 

C  THE  PROGRAM  MILL  RUN  HITH  ANY  AMOUNT  OF  DATA  ONLY  LIMITED  TO  DISK 
C  SPACE 

C 

C  ANY  BLANK  LINES  MILL  TERMINATE  THE  PROGRAM  IN  ERROR 

C 

c 

c 

DOUBLE  PRECISION  XPOS ,YPOS ,TIME1 , TIME! .TIMES ,SECD ,SECH .RATE 
DOUBLE  PRECISION  TIME4,  TIMES 

REAL  COURSE .PITCH!  2  )  .ROLL! 2  1  .SEClt  2  )  ,SEC2 .PIT ,ROL .SEC j! 2 ) 

REAL  OX.DY .DZ.OFFX.OFFY .OFFZ  .HEAD! 2  ) , LASEG , LOSEG 
INTEGER  MONTHl!  2  )  .0AY11  2  )  .YEARIf  2  I  .HOURK  2  )  ,MIN1(  2  ) 

INTEGER  MONTH2  .DAY 2  , YEAR2  .HOUR 2  ,MIN2 

INTEGER  MONTHS! 2  1 ,0AY5( 2  !  .YEARS! 2 ) .HOURS! 2  ),MIN3< 2 ) 

INTEGER  LAOEG, LAMIN, LOOEG.LOMIN 
C 

C  INITIALIZATION  OF  CONSTANTS 

C 

SECH«3600.D0 
SEC0*24 . D0*S£CH 

C 

0FFX«-5.345 

OFFY»-9.453 

OFFZ*-16.152 

C 

C  READS  TWO  SETS  OF  PITCH  ANO  ROLL  OATA 
C 

REA0<51,».END=100  )  MONTHl!  1  I.OAYl!  1 1.YEARl!  1  )  .HOURl!  1 )  .MINI!  1 ), 
9  SEC1I  1  !  .PITCH!  1  )  .ROLL!  I  ) 

TIM£1=DFL0AT(  DAY1!  1 )  )•SECO♦OFLOATf  HOURK  1 )  )*SECH 
*  ♦  D FLOAT!  MINK  1  I  1*60  .  00*DBlE  l  5ECK  1  )  ) 

C 

READ! 51 ,*  >END*100  )  MONTHlI 2  )  .DAYlf 2  ) , YE ARl! 2  )  .H0UR1! 2  )  .MINI! 2  ), 


n  o  o  n  o  o  o  o  o  ooo  o  o  o  non  o  ono  o  o  o  n 


*  SECK2),PITCH(2).ROLL(2) 

TIME2*OFLOATl  0AY1!  2 1 1*S£CO*OFLOAT(  HOURl!  2  )  l*SECH 

$  ♦  DFLOAT! MINI!  2  )  1*60 . DO+DBLE!  SECK  2  )  I 

READS  THO  SETS  OF  COURSE 

READ) 52  >END*100  )  MONTHS! 1 !  ,0AY5f 1 ) , YEARS! I )  ,HOUR3( 1 )  ,MIN3< 1 ) > 

$  SECS(l).HEAOU) 

TIME4*0FL0AT<  OAYS! 1 ) >*SECO*DFLOATI HOURS! 1 )  l*SECH 
$  ♦  0 FLOAT! MXN3! 1 )  1*60 ■ DO+OBLE! SECS! X )  ) 

READ! 52  >END'100  )  MONTHS! 2  )  .DAYS! 2  )  t YEARS! 2  )  .HOURS! 2  )  ,MIN5< 2  ’  > 

*  SECS!  2  I. HEAD!  2  I 
TIME5*OFLOAT( OAYS! 2  I  )*SECO*OFLOAT( HOURS! 2  I  )*SECH 

*  ♦  DFLOAT ! MXN3I 2  I  1*60 . DO+DBLE I  SECS! 2  I  ) 

READS  ONE  POSITION  OF  THE  SHIP 

20  READ! 50,£ND*100  I  MONTH2,OAY2,YEAR2,HOUR2,MIN2,SEC2, 

*  LA0E6 >  LAMIN , LASEG >  LOOE6  >  LOMXN  .LOSES 

TIM£S*OFLOAT(  DAY2  >»SECO*OFLOATI  HOUR 2  )»SECH 

*  ♦  DFLOAT! MIN2  )*60 . DO+OBLE l SEC2  ) 

POSITION  IN  TIME  BETWEEN  THE  TWO  SETS  OF  PITCH  AND  ROLL  DATA 
AND  THE  TWO  SETS  OF  COURSE  DATA 
25  CONTINUE 

POSITION  TO  EARLY  FOR  PITCH  ANO  ROLL  DATA 
IFfTIMES.GE. TIMED  SO  TO  40 

WRITE  I  6 . SO  )  M0NTH2 , 0 AY2 , YE AR2 , HOUR 2 , MINE .SEC 2 
SO  FORMAT!  '  PSN  3  1 ,12 , IX , 12 , IX ,14 ,1X,I2  , IX ,12 , IX ,F4 . 1 , 

S  •  REJECTED.  NO  PITCH  ANO  ROLL  DATA'  ) 

GOTO  2u 

POSITION  TO  EARLY  FOR  PITCH  ANO  ROLL  DATA 
40  CONTINUE 

IF1TIME3.GE.TIME4)  GO  TO  60 

WRITE! 6,50)  MONTH2.0AY2 , YEAR2 ,H0UR2 ,MIN2 ,SEC2 
SO  FORMAT!’  PSN  3  ' ,12 , IX ,12 , IX ,14, IX, 12 , IX , 12 , IX ,F4 . 1 , 

*  '  REJECTED.  NO  COURSE  DATA’) 

GOTO  20 

POSITION  LATER  THAN  THE  SECOND  SET  OF  PITCH  AND  ROLL  DATA 
60  CONTINUE 

IF! TIMES. LE.TIME2)  GOTO  70 

MOVES  THE  PITCH  ANO  ROLL  DATA  TO  THE  FIRST  SET 

MONTH!! 1  )*M0NTH1< 2  ) 

DAY!!  X  )*OAYK  2 ) 

YEAR1<  X  )*YEAR1<  Z ) 

HOURl!  X  I*  HOURl!  2  ) 

MINK  X  )*MIN1(  2  ) 

SECK  2  >*SEC1!  2  ) 

PITCH! X  )*PITCHI  2  ) 

ROLL! 1 >*ROLL< 2  ) 

TIME1*TIME2 

READS  A  SECOND  SET  OF  PITCH  AND  ROLL  DATA 

RE  AO!  5X  ,*,END»100  )  MONTHl!  2  I.OAYll  2  I.YEARK  2  1  .HOURl!  2  I, MINX!  2  ), 

*  SECl! 2  )  .PITCH! 2  I  .ROLL! 2  ) 

TIME2*0FL0AT(  OAYlt  2  )  )»SECO*OFLOAT!  HOURl!  2  I  l«$ECH 

t  ♦  DFLOAT!  MINK  2)  l»60 . 0O«OBLE l  SECK  2  ) ) 

GOTO  25 


"4 


non  n  non  on  non  non  non  o  n  n  n  non  non  non  non  non 


POSITION  LATER  THAN  THE  SECOND  SET  OF  COURSE  DATA 
70  CONTINUE 

IF( TINES. LE. TIMES)  GOTO  30 

MOVES  THE  COURSE  DATA  TO  THE  FIRST  SET 

MONTHS) I )*M0NTH5(  2  ) 

OAYS)  1)=DAY3<  2) 

YEARS) I  !*YEAR3( 2  ) 

HOURS) 1  )=H0URS>  2  ) 

MTNS) 1  )=MIN3( 2  ( 

SEC3I  2  )*SEC3) 2  ) 

HEAO) 1 )=HEA0( 2  ) 

TIME4*TIME5 

REAOS  A  SECOND  SET  OF  PITCH  AND  ROLL  DATA 

READ) 52 ,END*100  )  MONTHS) 2 ) ,0AY3( 2  )  .YEARS) 2  )  .HOURS) 2  )  .MIN3I 2  ) , 

*  SECS)  2), HEAO)  2) 

TIME5«DFL0AT( DAYS) 2 ) l«SECD*OFLOAT) HOURS) 2  )  >*SECH 

*  *  DFLOAT(MINS) 2  )  1*60 . 00+0BLE1 SECS) 2 ) ) 

GOTO  25 

POSITION  HITHIN  THE  TWO  SETS  IN  TIME 

SO  RATE*) TIMES-TIME1 )/<  TIME2-TIME1 ) 

PIT*PITCH(  1 )♦!  PITCH) 2 l-PITCH) 1 )  )»RATE 
ROL*ROLL<  X )♦(  ROLL) 2  J-ROLL) 1 )  1FRATE 
RATE*) TIMES-TIME4 )/) TIME5-TIME4  ) 

COURSE =HE AO ( X  )♦) HEAD) 2  ) -HEAD 1 1  )  )*RATE 

COMPUTS  THE  CORRECTIONS  AND  THE  POSITIONS  OF  THE  TRANSDUCER 

CALL  PTHI OFFX.OFFY ,OFFZ. COURSE .PIT ,ROL ,DX ,0Y ,DZ  ) 

COMPUTES  THE  UTM  COORDINATES 

CALL  GPUTMI LADEG . LAMIN , LASEG , LODEG . LOMIN , LOSEG , YPOS ,XPOS  ) 

XPOS*XPOS+OX 

YPOS=YPOS*OY 

CONVERTS  BACK  TO  GP 

CALL  UTMGP ( YPOS ,XPOS . LAOEG , LAMIN , LASEG , LODEG . LOMIN , LOSEG ) 
OUTPUTS  THE  RESULT 

WRITE) 7)  MONTH2 ,DAY2 .YEAR2 .HOUR 2 ,MIN2 ,SEC2 , 

*  LADEG . LAMIN , LASEG , LODEG . LOMIN , LOSEG 

WRITE) 6,90)  M0NTH2 >DAY2 .YEAR2 .H0UR2 ,MIN2 ,SEC2 , 
a  LADEG , LAMIN . LASEG , LODEG . LOMIN , LOSEG 

90  FORMAT) 21 XX, 12  )  ,XX,I4 ,1* ,12 , IX ,12 , IX ,F4 . 1 , 2) XX, IS ,XX,I2 ,1X,F7.4 ) ) 

REAOS  ONE  MORE  POSITION 
GOTO  20 

END  OF  PROGRAM 

100  CONTINUE 
STOP 
END 

SUBROUTINE  PTH) XOFF ,Y0FF ,Z0FF .COURSE .PITCH .ROLL , OX ,DY ,DZ ) 

THIS  SUBROUTINE  COMPUTES  THE  CORRECTIONS  TO  THE 
COORDINATES  DUE  TO  THE  OFSET  OF  ANTENA  IN  RELATION 
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c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


TO  THE  TRANSDUCER 

PI*ARCOSI-1.0) 

ROU.=ROU.*PI/180. 

PITCH»PITCH«PI/180 . 

COURSE «C0URSE»PI/180 . 

"SS  StL°"is”“l  “t*  t°  mTC"  ™E  »«« ™* 

ROLLs-ROLl 

COSPHJs  COS( ROLL ) 

SINPHI  =  SINI ROLL ) 

COSKs  COS  (COURSE) 

SINKs  SINI COURSE) 

COSWs  COS( PITCH  ) 

SINhs  SINi PITCH) 


OX 


■  XOFF  a  COSPHI  a  COSK  ♦ 

YOFF  a  (  COSH  a  SINK  ♦  SIlsH 


ZOFF  a  (  siW  a  SINK  -  COSH  a  SINPHI 


»  SINPHI  a  COSK  )  ♦ 


*  COSK  ) 


OY  =  XOFF  a  COSPHI  a  SINK  a  („i.  )  * 

YOFF  »  I  COSH  a  COSK  -  SIFK  a  SINPHI 
ZOFF  a  (  Sim  a  COSK  ♦  COSH  a  SINPHI 

OZ 


a  SIMC 
•  SINK 


)  ♦ 


XOFF  a  SINPHI  - 
YOFF  a  (  SINH  a  COSPHI  ) 
ZOFF  a  (  COSH  a  COSPHI  ) 
RETURN 
END 


^ROUTINE  UTM6PI NORTH  ,E  AST,  LAOEG>LAHIN*LASE<5»LOOES»LCMIN,  LOSE® ) 
SK  PRECISION  ^,"’f:?P>S-Rl'ES«I.ESqP,RN,RP>KO 

SSS  mIm!X0M 

s  srsr-  or  w  «  ■*“  * 

CMs-123.000 

A*6378135. 000 
R*298.26D0 

K0=0 . 999600 

PI*OARCOS( ”1.00 ) 

B  *  Aa(  R*1.00  )/R 


N  « 

AP 

BP 

CP 

OP 

EP 


(A-B)/(  A  *8) 

!  t#ii^:O?'NUS'DO/^-D0*>N*a2-N*a3  )»81. 00/64  00a( Naa^.Maae  i  > 

•  W  OO/iroOaii^^^rooTnn^5*^  ’  *55 ' 00/64  '  00*N<H*5 1 

31S ,  00/512 . 00a*a(  Naa<t-N**5  | 


FIRST  APROXI NATION  OF  PHI 
PHIsnoRTH/30 . 800/3600 . 00 
COMPUT  TRUE  MERIDIONAL  DIST  ANO  APROXIMATE  PHI 
00  100  1*1,8 


* 


a*»"*  9‘  »‘9,%  .9.%  A  A  AA  A 


oooooo  o  n  noon 
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PHIMIN*PHI»60 . 00*2 . 906662086660 -4 
RPHI SPHI/160 . 00*PX 

S  *  AP*PHIMIN-BP*0SINI 2 . D0*RPHI  )  *CP*0SINI 4 . D0*RPHI I 
*  -DP*0SINI 6 . 00*RPHI I *EP*0SIN( 6 . D0*RPHI ) 

R1  *  K0*S 
0ELTA*N0RTH-R1 
PHI*<0ELTA/30. 8/3600.  )*PHI 
100  CONTINUE 

RPHI =PHI/180 . 00*PI 


SINSEC  2  I  1.00/360000 1/180. 00»PI 
SINSEC  s  0SINI SINSEC) 

ESQ  *  <A**2-B*»2)/A**2 
ESQP  *  ESQ/I  1.00-ESQ  1 

RM  a  A*( 1. DO-ESQ  )/( DSQRT( l.DO-ESQ*OSIN( RPHI  )»*2  )  )**3 
RP  =  RH*I  1 . DO*ESQP*OCOS( RPHI  )**2  ) 

R7*(  0TAN(  RPHI  I/I  2 . 00*RP**2*SINSEC  11*11. 00 ♦ES«P*OCOS(  RPHI  )**2  I 
*  *1.012/KO**2 

R8=IDTANI RPHI  I/I  24. 00*RP**4*SINSEC  I  1*1 S. 00*3. D0*DTAN( RPHI l**2 
4  +6 . 00*£SQP*CCOSI  RPHI )**2-6 . D0*ESQP*0SINf  RPHI  l»*2-3 . 00*ESQP**2 

$  »0COSI RPHI  )**4-9.D0»ESQP**2*0C0SI RPHI  )**2*0SIN( RPHI  l**2 I 
4  *1 . 024/K0**4 

R9*l .  OO/OCOSI  RPHI  l/l  RP*SINSEC  1*1 . 06/KO 

R10  =  l . DO/DCOSl RPHI  1/16. 00*RP**3*SINSEC  1*1 1 . 00*2 . 00*0TAN( RPHI  )**2 
4  *ESQP*OCOSl  RPHI  )**2  1*1 . D16/K0**3 

EPRIHE *E AST -S00000 . DO 
Q  *  ,000001D0*EPRIME 
Q2  *  0**2 
Q3  *  0**3 
04  *  Q**4 

Q5  s  Q**S 
Q6  *  Q**6 

D6NUM*Q6*DTANI RPHI  I 

06*1 DbNUM  I/I  720 . 00*RP**6«SINSEC >*( 61. 00*40 . D0»DTANI RPHI  )**2 
4  *45  - 00*0TAN( RPHI  1**4*107. 00*ESQP*OCOSI RPHI  l»*2 

4  -1*2. 00»ESQP*OSINI RPHI l**2-45. 00*CSQP*0TAN( RPHI )**2 

4  »OSINt  RPHI  i**2  1*1 . 036/K0**6 
ES  =  i  05*1. OO/OCOSI RPHI  I  I/I  120 . 00*RP**5*SINSEC  1*1 5 . 00*28 . DO 
4  *0T AN! RPHI  )**2 *24 . 00*OTANI RPHI  l**4*6 . 00*ESQP*0COSl RPHI  )**2 

4  *B  00»ESOP*OSIN( RPHI 1**2  1*1 . D30/X0**5 

OPHI *1  -R  7*02  *R8*O4-06  1/3600.00 
OLAM*t  R«*0-R  10*03 *E5  1/3600 .  DO 
PHI *PHI *OPHI 
LON«CM*OLAH 

CALL  DMSI PHI, LAOEG, LAMIN, LASEG I 
CALL  DMSILON.LOOEG.LOMIN.LCSEGl 
RETURN 
END 

SUBROUTINE  GPUTM1 LAOEG .LAMIN .LASEG ,LOOEG , LOHIN , LOSEG .NORTH .EAST  I 

DOUBLE  PRECISION  A ,R ,N,AP ,BP ,CP ,0P ,EP ,S ,R1 , ESQ, ESQP .RM.RP ,K0 
DOUBLE  PRECISION  R2 ,R3 ,R4 ,R5 ,P  .P2  ,P3  ,P4.P5  ,P6  .A6  ,B5  .SINSEC 
DOUBLE  PRECISION  PHI ,OLAM  .NORTH , EAST  ,B .PHIMIN, PI , LON 
REAL  LASEG. LOSEG 

THIS  SUBROUTINE  COMPUTS  THE  UTI4  COORDINATES  OP  GP  IN  HCS  72 
IN  ZONE  10  CENTRAL  MERIDIAN  123  00  00  H 

CM.-123.000 

PHI.OFLOATI LAOEG  )*OPLOAT( LAMIN  1/60 . DO+DBLE I LASEG  1/3600 . 00 
l0N«0FLQATI LOOEG  l*OFLOATt  LOMIN  1/60 . DO*DBLEI LOSEG  1/3600 . 00 
LON.-LON 

OLAM.I LON-CM  1*3600 . 000 


>  '  *  .*  ,» 


vniiiiiiui 
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A*63781S5.0D0 
R*298. 2400 
C 

K0»0. 999600 
C 

8  *  A*<R-1.D0)/R 
C 

N  *  (  A-8 >/( A*B  ) 

AP  *  A»(  ( l.D0-N>*5. 00/4. 00*(  N**2-N**3 >*81. 00/64. 00*( N**4-N**5 > > 

BP  s  S. 00/2. 00*A#(  <N-N**2  1*7. 00/8. 0Q*t N**3-N**4  1*55. 00/64. 00*N**5 > 
CP  *  15. 00/16. D0*A#(N**2-N**S*S. 00/4. D0#( N**4-N**5 >) 

OP  *  15. 00/48. 00*A*<N**3-N**4*11. 00/16. 00#N**5 ) 

EP  =  315. 00/512. 00#A«< N**4-N**5 > 

PMIMIN  >  PHI*60. 00*2.408882086660-4 
C 

PI=OARCOS( -1.00 ) 

C 

PHI spHI/180 . D0*PI 
C 

S  3  AP*PHIMIN-BP*OSIN( 2 . 00*PHI  1  *CP*DSIN( 4 . 00*PHI  l 
t  -  DP*DSIN(  6 . D0*PHI J *EP*0SIN( 8 ■ D0*PHI  ) 

R1  >  KO*S 
C 

SINSEC  3  ( 1.00/3600.00  1/180. Q0*PI 
SINSEC «0SIN<  SINSEC I 
c 

ESQ  ■  ( A 6*2-86*2 )/A**2 
ESQP  >  ESQ/I  1.00-ESQ) 

RK  3  a*( 1 . OO-ESQ  )/( DSQRTI 1 . DO-ESQ*OSINI PHI )**2 ) )**S 
C 

RP  s  RM*t l.DO*ESQP*DCOSIPHI  >**2 ) 

R2  3  RP»OSIN< PHI  )*OCOSI PHI  )*SINSEC*«2/2 . DO*KO*1 . 008 
R3  *  SINSEC**4*RP*0SIN( PHI  )*OCOS( PHI  )**S/24 . D0*( 5 . DO-DT AN( PHI  )**2 
t  *  9. 00*ESQP*0COS( PHI  >#*2*4 . 00*£SQP*ESQP*DCOS( PHI  )**4  l«KO*l .016 
R4  3  RP«DCCSI PHI  )*SINSEC*K0*1.Q4 

R5  3  SXNSEC»*3*RP*OCOSI  PHI l**S/6 . 00*(  1 . 00-DTANI  PHI )*»2 

*  ♦  ESQP*OCOS( PHI  )**2  l*KO*1.012 
C 

P  3  .000100H0LAM 
P2  *  p**2 
P3  *  P**3 
P 4  *  P**4 
P5  a  p**5 
P6  «  P**6 
C 

A6  3  P6*SINSEC»*6»RP*0SIN( PHI  )*DCOS( PHI >**5/720.00 

*  »  I  61 . DO-58. D0*0TAN( PHI )**2*0TANI PHI )**4 

i  ♦  270 . 00*ESQP*OCOS( PHI  )**2-330 . 00*ESQP*0SIN( PHI )**2 )«K0*1 . 024 
85  *  P5*SINSEC**5»RP*0C0S(PHI >*#5/120. 00*15. 00-18. 00*OTAN1PHI  1**2 

*  *  0TAN1PHI  >#*4*14. 00*ESQP*DC0S(  PHI )**2 

*  -  58.00*ESQP*0SIN(PHI l**2 >**0*1.020 
NORTH *R1*R2*P2*R3*P4*A6 

EASTs( R4*P *RS»P3 *85  )  *500000 . DO 
RETURN 
ENO 
C 

SUBROUTINE  DHSI DEC  »LDEG .MIN , SEC > 

C 

DOUBLE  PRECISION  DEC.XNUM.XMIN 

C 

XNUM«OABSI DEC ) 

IOE6«OINT <  XNUM ) 

XMIN»<  XNUM-OFLOATI  LOEQ  I  1*60 . 00 
MIN30INT(XMIN) 

XNUM=(XMIN-DFLOATlMIN>  >*60.D0 
SECsSNGLI XNUM ) 

IF  (10EG.GE.36O)  LDEG*LDEG-560 

RETURN 

ENO 

/* 
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APPENDIX  D 

PROGRAM  PLOT.  SOURCE  LISTING 


//EZEQUIEL  JOB  (  0812  ,M« ) ,  ‘  EZEQUIEL  ‘  ,CLASS*C 
//  EXEC  FRTVCLGP 

//FORT . SYSIN  00  * 

C 

C  PROGRAM  PLOT 

C 

C  RUNS  IN  FORTRAN  VS 

C 

C  AUTHOR:  AUGUSTO  EZEQUIEL 

C 

C  DATE  :  26  JANUARY  1987 

C 

C  THIS  PROGRAM  MAKES  A  PLOT  OF  THE  POSITIONS  OF  THE  SHIP 

C 

C 

DOUBLE  PRECISION  XPOS,YFOS,XLEFT,YLEFT,TIME1,TIME2,SECD,SECH 
REAL  SECS ,XPLT , YPLT ,8LXH ,BLYH .SCALE .VALUE . LASEG ,LOSEG 
INTEGER  MONTH, YEAR, DAY. HOUR, MIN, IPEN.LAOEG, LANIN, LOOEG.LOMIN 

PLOTTER  INITIALIZATION 

CALL  PLOTS  (0,0,0) 

CALL  PLOT  (2. ,2. ,-3) 

IPEN*3 

DIMENSIONS  OF  SHEET,  LEFT  CORNER  AND  SCALES 

BLXH=36. 

BLYH=49. 

XLEFT =564800 . DO 
YLEFT  =4037000 . 00 
SCALE *1 . /10000 . 

INITIALIZATION  OF  CONSTANTS 

TIME 1*0. DO 
SECH=3600 . 00 
SEC0=24.D0«SECH 

GRIO 


CALL  PLOT  (0.1, 0.1, 3) 

CALL  PLOT  (0.1 ,BLYH*0 .1,2) 

CALL  PLOT  ( 8LXK-Q . l,BLYH-0 .1,2) 

CALL  PLOT  ( BLXH-0 . 1 ,0. 1,2 ) 

CALL  PLOT  10.1, 0.1, 2) 

TITLE 

CALL  TITLE 
SCALE 

CALL  METER)  SCALE) 

READ  THE  DATA 
10  CONTINUE 

READ*  51 ,EN0:50  (MONTH ,0AY .YEAR .HOUR , 

*  MIN, SECS, LAOEG.LAMIN. LASEG, LODES, LOMIN.LOSEG 


XI) 


u  u  o  oou  o  o  u  oou  o  o  o  o  u  o  o  o  o  u  u  o  oou 


TfCYTT 


rv*  <r,T.  ■  ,  * ' ' 


COMPUTES  THE  UTM  COORDINATES 

CALL  GPUTM!  LAQEG , LAMIN , LASEG , LODEG .LOMIN , LOSEG , YPOS >XPOS ) 
COMPUTES  THE  TIME  IN  SECS 

TIM£2=SEC0*FL0AT(  DAY  )+S£CH*FLOAT(  HOUR  )+60 .  D0*FL0AT( MIN  USECS 

COMPUTES  THE  PLOTTER  COORDINATES 

XPLTMXPOS-XLEFT  )*SCALE*100. 

YPLT  =  ( YPOS-YLEFT )*SCAL£*100. 

TESTE  IF  INSIDE  AREA 

IFIXPLT.LT.O.  )  GOTO  50 
IF< YPLT.LT.O. I  GOTO  50 
IFIXPLT.GT.BLXH)  GOTO  50 
IFIYPLT.GT.9LYH I  GOTO  50 

TESTS  IF  POSITIONS  ARE  AWAY  MORE  THEN  10  SECS  IN  TIME 

IFITIME2-TIME1.GT.10.00)  IP£N=5 
TIM£1=TIME2 

PLOTS  THE  POSITION 

CALL  PLOT(XPLT,YPLT,IPEN) 

XPEN=2 

PLOTS  THE  TIME  EVERY  OS  MINUTES 

IF(MOO<MIN,05I.GT.O.OR.IFIXISECS>.GT.O)  GO  TO  10 
VALUE »FL0AT( HOUR  l»100*FLOAT! MIN ) 

CALL  NUMBER! XPLT*0.  IS,' YPLT*0. 15,0.25,  VALUE,  0.  ,-H 
CALL  SYMBOL! XPLT.YPLT, 0.15,5,0.  ,-l) 

GOTO  10 

PEN  UP  WHILE  THE  POSITIONS  ARE  OUT  OF  THE  SHEET 

50  IPENs? 

GOTO  10 

END  OF  PLOT 


50  CALL  PLOT!  0. ,0. ,*999) 

STOP 

END 

C 

SUBROUTINE  METER! SCALE) 

C 

X0*24.0 

Y0=2.0 

CALL  PLOT(XO,YO,-5) 

C 

DO  10  J*l,12 
XP«  FLOAT! J  )-1.0 
CALL  PLOT! XP, 0.0, 5) 

CALL  PLOT! XP ,0. 2S,2  ) 

10  CONTINUE 
C 

DO  20  JS1 , 9 
XP  =  FL0AT(  J)*0.1 
CALL  PLOT! XP  ,0.0.5) 

CALL  PLOT! XP ,0 . 2 1 2  ) 

20  CONTINUE 
C 

CALI  PLOT ( 0 . ,0.0,5 ) 

CALL  PLOT! 11. ,0. ,2 ( 

SI 


f 

I 

I 

I 


r>  o  r>  o  o  o  n  noon  non  non  o  o  oo 


CALL  PLOTtO.  ,0.2,5) 

CALL  PLOT! 11. > .2,2  ) 

C 

VALUC*1. /SCALE/100. 

CALL  NUNBER!  -0.25,0. 27, 0.25, VALUE, 0. 0,-1) 

CALL  NUPBER!  0.95 ,0.27 ,0.25, 0.0,0. 0,-1 1 
VALUE  *5 . 0/SCALE/100 . 

CALL  NUMBERI5. 775. 0.27, 0.25, VALUE, 0. 0,-1  ) 

VALUE=10 .0/SCALE/100 . 

CALL  NUMBER! 10 . 7 ,0 . 27 ,0 . 25 , VALUE ,0 . 0 , -1 ) 

CALL  SYMBOL! 5 .5,-0. 50 , . 25.6HMETERS ,0 . 0,6) 

C 

CALL  SYMBOL! 5. 5, 24. 0,2. 0,62,0. 0,-1 ) 

CALL  SYMBOL! 5. 7, 24. 5, 0.4, 85,0. 0,-1) 

X0=-X0 

YO=-YO 

CALL  PLOT (  XO ,  YO  >-3  1 
C 

RETURN 

EN0 

C 

SUBROUTINE  TITLE 
C 

XO=26.5 
Y0S43 . 05 

CALL  PLOT! XO,YO, -3  ) 

C 

CALL  SYMBOL! 00. 0,3. 5,0. 25, 29H  NAVAL  POSTGRADUATE  SCHOOL  >0.,29) 
C 

CALL  SYMBOL! 00.0,3.0,0. 25.29HSEAFLOOR  BENCHMARK  EXPERIMENT, 0. ,29) 


CALL  SYMBOL!  00 . 0 ,2 .5,0. 25>29H  PHASE  II  ,0.,29) 

CALL  SYMBOL! 00 .0 ,2 .0,0 . 25.29H  R/V  POINT  SUR  ,0.,29) 

CALL  SYMBOL! 00. 0,1. 5,0. 25, 29H  15/17  AUGUST  1986  ,0.,29) 

»  CALL  SYMBOL! 00 .0 ,1 . 0 ,0 . 25 ,29H  GESAR  SOLUTIONS  ,0.,29) 

CALL  SYMBOL! 00 . 0 , 1 . 0 ,0 . 25 , 29H  TI4100  ,0.,29) 

»  CALL  SYNBOL!  00 . 0  >  1 . 0 ,0 . 25 » 29H  MINI  RANGER  FALCON  ,0.,29) 

CALL  SYMBOL! 00 . 0 ,0 . 5 ,0 . 25 , 29H  ANTENNA  POSITIONS  ,0.,29) 

»  CALL  SYMBOL! 00. 0 ,0 .5,0. 25, 29H  TRANSDUCER  POSITIONS  ,0.,29) 

CALL  SYMBOL!  00. 0,0. 0,0. 25, 29H  LCDR  AUGUSTO  EZEWJIEL  ,0.,29) 


CALL  PLOT! -0.2, -0.2, 3) 
CALL  PLOT! -0.2,3.95,2) 
CALL  PLOT! 7.5,3.95,2 ) 
CALL  PLOT! 7. 5,-0. 2 ,2 ) 
CALL  PLOT! -0.2, -0.2, 2 ) 

XO*-XO 

YO*-YO 

CALL  PLOT! XO,YO,-3 ) 

RETURN 

ENO 


SUBROUTINE  GPUTM!  LAOEG , LAMIN , LASEG , LOOEG , LOMIN , LOSEG .NORTH , EAST  ) 
OOU8LE  PRECISION  A,R,N,AP,BP,CP,DP,EP,S,R1  ,ESG,£SGP  ,RM,RP ,KO 


UMUKMt 
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c 
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c 

c 

c 
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c 
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c 

c 


c 


c 


c 


/* 


DOUBLE  PRECISION  R2 ,R3 ,R4,R5,P,P2 ,P3 ,P4,P5>P6 ,A6 ,85, SINSEC 
DOUBLE  PRECISION  PHI ,0LAM, NORTH .EAST .B.PHIMIN, PI , LON 
REAL  LASEG.LOSEG 

THIS  SUBROUTINE  COMPUTS  THE  LTTM  COORDINATES  OF  GP  IN  MGS  72 
IN  ZONE  10  CENTRAL  MERIDIAN  123  00  00  W 

CMs-123.000 

PHI=0FLOAT( LADEG  l+DFLOATI LAMIN  )/60 . D0+0BLE I LASEG  1/3600 . DO 
LON=0FLOAT( LOOEG  l+OFLOATI LOMIN  1/60 . 00+0BLE ( LOSEG  1/3600 . 00 
LON=-LON 

0LAM=( LON-CM  1*3600 . 000 

A=6378135.000 
Rs298. 2600 

K0*0. 999600 

S  =  A*<  R-l .  DO  )/R 

N  =  ( A-B  )/( A*8  ) 

AP  =  A*l  (  1.00-N)  +5. 00/4.00*1  N**2-N**3  1+81. 00/64. 00*C N**4-N**5 1  ) 

BP  =  3. 00/2. D0*A*( <N-N**2  1*7. 00/8. 00*( N**3-N**4  1*55. 00/64. 00*N**5 1 
CP  =  IS . 00/16 . D0*A*( N**2-N**3+3 . 00/4 .  D0*( N**4-N**5 1 ) 

OP  =  35. 00/48. D0*A*(N**3-N**4*11.  DO/16.  00*N*#5) 

EP  =  315. 00/512. 00*A*(N**4-N**S  1 
PHIMIN  =  PHI*60. 00*2. 908882086660-4 

PI =DARCOS( -1 . 00  ) 

PHI=PHI/180.00*PI 

S  =  AP*PHIMIN-BP*OSIN< 2.00*PHI  !*CP*DSINI 4.00*PHI ) 

$  -  0P*0SIN(6.00*PHI )*EP*0SINI8.D0*PHI ) 

R1  =  KO*S 


SINSEC  =  <  1. 00/3600. 00 1/180. 00*PI 
SINSEC-OSINf  SINSEC ) 

ESQ  =  (A*#2-B**2)/A**2 
ESQP  *  ESQ/I  1.00-ESQ) 

RM  =  A#<1.00-ESQ)/<OSQRT<1.00-ESQ*DSIN(PHI)**2))**3 
RP  =  RM*<1.0Q*ESQP*0COS<  PHI  1**2  1 

R2  =  RP*0SIN(  PHI  )*DCOSI  PHI  )*SINSEC**2/2 .  D0*K0*1 . 008 
R3  =  SINSEC**4*RP*0SIN< PHI )*DCOSI PHI  )**3/24 . D0*( 5 . DO-OTANf PHI )**2 
i  ♦  9.0O*ESQP*DCOS(PHI)**2*4.00*ESQP»ESQP*0C0S(PHIl**4)*KO*l.D16 
R4  *  RP*OCOS( PHI )*3INSEC*K0*1 . D4 

RS  -  SINSEC**3*RP*OCOS< PHI  )**3/6.D0*t l.DO-DTANI PHI )**2 

*  ♦  ESQP»OCOS(PHI  )**2)*KO*1.012 

P  =  . 0001DO*DLAM 
PZ  -  P**2 
P3  =  P**3 
P4  *  P**4 
PS  s  p**5 
P6  *  P**6 

A6  *  P6*SINSEC**6*RP*0SIN(PHI  )*DC0SI PHI  )**5/720.D0 

*  *  (61. 00-58. 00*0TAN(PHI  )**2*0TANf PHI  >**4 

S  *  270.00*ES8P*DCOSt PHI  1**2-330. DO*ESQP*DSIN< PHI  )**2  )«KO*l.D24 
85  *  P5*SINSEC**5*RP*OCOS(  PHI  1**5/120. 00*15.00-18. D0*0TAN(PHI  )**2 
$  ♦  OTANI PHI  )**4*14 . DO*ESQP*OCOSl PHI  i**2 

*  -  58 . 00*ESQP*OSIN( PHI  1**2  l*KQ*l . 020 
NORTH =R1*R2*P2*R5*P4*A6 

EAST*! R4*P+R5*P3*B5 1*500000 . 00 

RETURN 

ENO 


//GO. PLOTPARM  00  * 


S3 
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’YMAXSS1-  'SCALE  =  1-  >^2.540  4END 
//GO. FT51F001  00  DSN=MSS.S0812 . GPS. POS. ANTENNA  >0ISP=SHR 
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APPENDIX  E 

PROGRAM  CVF1CA.  SOURCE  LISTING 


//JOBCVFI  JOB  <081Z,9999I,‘EZEQUIEL‘ ,CLASS*C 

//•MAIN  0RG=NPGVM1 . 0812P ,LINES  =  < 99 )  ,CAR0S=< 99 ) 

//•FORMAT  PR,00NAME=GO.FT06F001, 

//•FORMS =SEP1 

//  EXEC  FORTVCLG 

//FORT . SYSIN  00  * 

PROGRAM  CVFICA 

C  AUTHOR:  AUGUSTO  EZE OUZEL 

C  DATE:  MARCH  17,  1987 

C 

C  DESCRIPTION: 

C  THIS  PROGRAM  CONVERTS  THE  FICA  FILES  INTO  THE 

C  INPUT  FORMAT  OF  THE  KALMN  PROGRAM  ANO  GIVES 

C  GENERAL  INFORMATION  ABOUT  THE  DATA 
C 

DOUBLE  PRECISION  FPII 500  )  ,TP ,PB  ,H0  ,RL10FF,RL20FF 
DOUBLE  PRECISION  01 .02 ,PI ,C ,V0S,K1 ,K31 ,K32 ,TLDLL ,TLPLL 
DOUBLE  PRECISION  POBPC  ,P08C  ,CN1(  4  1  ,CN2l  4  )  ,CR1<  4  )  ,CR2<  4  ) 

DOUBLE  PRECISION  OOPli 4 ) ,00P2( 4  )  ,BLC ,BPDC ,SGRlt 4 > ,SGR2t 4 ) 

DOUBLE  PRECISION  SGD1( 4  J ,SG02( 4 ) ,K2,K4,TT ,M,AG,ESO 
DOUBLE  PRECISION  TOCS , AS( 3  t ,A0C ,CRSS, DNS >MOS,CUCS,ES, CUSS 
DOUBLE  PRECISION  SOAS ,T0ES .CICS.OMEGS ,CISS ,IOS ,CRCS ,WS 
DOUBLE  PRECISION  OLAT ,OLOG ,DHT 
DOUBLE  PRECISION  OMEDS ,A0E ,IOOTS,X,V ,Z,R0BL 
REAL  LASEC , LOSEC ,SEC 

INTEGER  INTGI 500  J .BLOCK , TYPE ,NCI ,NII ,NFI , LAMIN , LOMIN , LAOEG , LOOEG 
INTEGER  TRACER, MOVELI  4  )  ,ISAT<  4  )  ,IN0EX,SAT ,IKKNOS,IEDATS,NSAT 
INTEGER  10, STOAY, MONTH, YEAR, HOUR, MIN, DAY 
INTEGER  STATI4) 

CHARACTER»8  CHI I  500  )  .TITLE ( 10  I 
CHARACTER«5  BLK 
LOGICAL  SAVE, FIRST 

SETS  DEFAULT  METEO  DATA 

OATA  TP.PB.HO  /IS . , 980 . ,75 . / 

SETS  THE  PARAMETERS  OF  WGS72 

OATA  AG,  ROBL  /«378. 135000 ,298. 26D0/ 

FLAG  TO  SAVE  DATA  ONLY  WHEN  THERE  IS  NAVIGATION  OATA 
FOR  ALL  TRACERS 

SAVE*. FALSE. 

FLAG  TO  SAVE  INFORMATION  DATA  AS  TYPE  1  ONLY  ON  FIRST 
CASSETTE 

FIRST*. TRUE. 

ECCENTRICITY  SQUARED  OF  REFERENCE  ELLIPSOID 
ESQ*(  2.00-(  l.OO/ROBLl  i/ROBL 


<31  =  154 . 00 
<32=120,00 

COMPUTS  THE  VALUE  OF  PI  IN  DOUBLE  PRECISION 


35 


non  o  o  n  o  non  o  n  o  on  n  non  noonoonno 


C 


Pl=OARCQS! -1.00) 


SPEED  OF  LIGHT  (KM/SEC) 

0=299792.45600 
NOMINAL  SATE LITE  FREGUENCY 
VOS=10 . 2306 

CALCULATION  Of  CONSTANTS  TO  BE  USED  IN  THE  SGIMA  COMPUTATIONS 
Kl  =  ( C/VOS  l**2 

K31=!C/( 2.00*PI*VOS»91> )**2 
K32=!C/!2.00*PI#VOS»Q2)  I»*2 


RL10FF«-6000.00 
RL20FF=7600 . 00 


M=4.D0 
TL0LL= . 700 
TLPLL= . 7DO 
PDBPC=100 . 00 
PDBC=100 . DO 

SETS  THE  FLAGS  FOR  NOT  HAVING  THE  EPHEMIS  DATA 

00  3  1=1,4 

3  STATU  1*0 

PRINTS  OUT  THE  IDENTIFICATION  OF  THE  PROGRAM 


READS  THE  TITLE  FROM  THE  FILE  AND  CONTROL  INFORMATION 

READ! 31,5,£NO*2000  ITITLE 
WRITE  I  6,5)  TITLE 
READ!  31 ,5,EN0=2000 )TITLE 
WRITE) 6 ,5 ,END=2000  )  TITLE 
5  FORMAT! 10A8) 

REA0I31,»,EN0«2000)  MONTH ,ST0AY .YEAR 
WRITE! 6,6)  MONTH, STOAY, YEAR 
6  FORMAT!///, 

9  '  OATE  OF  STARTING  GPS  WEEK  !  MONTH  DAY  YEAR)  ',13, 13, IS) 

INPUTS  THE  DATA  FROM  THE  COWANO  FILE 

READ! 31 ,» >END*2000 )  NSAT ,  K2,  K4 
READ! 31,»,END*2000 )  X,Y,Z 
WRITE ! 6 , 1 )  NSAT ,K2 ,K4 ,X ,Y ,2 
1  FORMAT  I '  NUMBER  OF  TRACERS  ‘,120,/, 

9  '  RANGE  SIGMA  BIAS  FACTOR  ’ ,E20.8,/, 


$6 


ft  '  DOPPLER  SIGMA  BIAS  FACTOR  ',£20. 8,/, 
9  '  ESTIMATE  RECEIVER  POSITION',/, 

9  '  X  COORD' ,F20.6,/, 
9  '  Y  COORD' ,F20.b,/, 
ft  '  2  COORD' >F20.6,/) 


C 


c 

C  READS  THE  CONTROL  BLOCK  INFORMATION 

C 

10  READ! SO >20  ,END=2000  )  BLK .BLOCK ,NFI , Nil , NCI 
20  FORMAT!  AS  ,415  I 

C 

C  READS  THE  FLOATING  PART  DATA  IF  ANY 
C 

IF  (NFI.EQ.O)  GOTO  35 

READ! 50 ,30 ,END*2000  )  ( FPIII) ,1*1  ,NFI ) 

30  FORMAT! 4020. 16) 

C 

C  READS  THE  INTEGER  PART  IF  ANY 

C 

35  CONTINUE 

IF  INII.EQ.O)  GOTO  45 
READ! 50 ,40 ,EN0*2000  !  ( INTG!  I ) ,1*1 ,NII ) 
40  FORMAT! 6112) 

C 

C  REAOS  CHARACTER  PART  DATA  IF  ANY 

C 

45  CONTINUE 

IF  tNCl.EQ.O)  GOTO  55 

READ! 50 ,50 ,ENO=2000  )  (CHI! I  ), 1*1, NCI > 

50  FORMAT! 10! A8 ) ) 

55  CONTINUE 

IF  (BLOCK. NE. 6  I  GOTO  100 


C 

C  USER  EPOCH  TIME  OF  PSEUDO  RANGE 

C 

TT *FPI ( 3  ) 

C 

C  COMPUTES  THE  OATE 

C 

CALL  SUB! TT, DAY, HOUR, MIN, SEC ) 

DAY =OAY +STDAY 
C 

C  IF  NO  EPHMERIS  DATA  AVOIDS  COMPUTATIONS 
C 

IF! .NOT. SAVE  )  GOTO  10 
C 

C  LI  CARRIER  SIGNAL  TO  NOISE 

C 

CNK  1  )*FPI(  4  ) 

CNK  2  )*FPI(  5  ) 

CNK  3  )*FPIt  6  ) 

CN1I4  )*FPI(  7 ) 

C 

C  L2  CARRIER  SIGNAL  TO  NOISE 

C 

CN2! 1 )*FPI(8 ) 

CN2I 2  )=FPIl 9  ) 

CN2!  3  )  =  FPI<  10  I 
CN2<  4  )=FPI( 11 ) 

C 

C  LI  PSEUOO  RANGE  <  KM  IN  FICA  FILES  CONVERTED  TO  SECONDS  FOR  KALFtJ2 ) 
C 

CR1!  I  )*FPI(  12  l/C 

S7 


CRH  2  )*FPI(  13  )/C 
CRll  5 )=FPI( 14  )/C 
CR1!4)*FPI(15)/C 
C 

C  L2  PSEUOO  RANGE  t KM  IN  FICA  FILES  CONVERTED  TO  SECONDS  FOR  KALFN2 ) 
C 

CR2tll*FPI(16)/C 
CR2(  2  )=FPI<  17  )/C 
CR2I3  )=FPI(  la  1/C 
CR2I  4  )=FPI (  19  1/C 
C 

C  LI  CARRIER  DOPLER  PHASE 
C 

DOPK  1  )*FPII  20  ) 

DOPK2»FPI(21) 

DOP1I 3  )=FPI1 22  ) 

OOPll  4 )*FPI(  23  ) 

C 

C  L2  CARRIER  DOPLER  PHASE 

C 

D0P2I  1 )=FPI( 24  ) 

OOP2( 2  )=FPI! 25  ) 

OOP2I3)=FPI<26) 

OOP2I  4  )=FPI<  27) 

C 

C  SV  PRN  OF  EACH  TRACER 

C 

IS ATI  1  l*INTG( 1  I 
ISATl 2  )=INTG( 2  ) 

ISAT(  3  )=INTG( 3  ) 

ISATl 4  l=INTG<  4  ) 

C 

C  L1.L2  QUALITY  FACTOR  I  TRACER, FREQUENCY  ) 

C 

MOVE  LI  1  )*INTG(  9  l-INTGI  13  1 
MQVELI  2  )*INTG<  10  l-INTGI  14  ) 

MOVE  LI  3  )=INTG(  11  l-INTGI  IS  I 
MQVELI  4 l=INTG< 12  l-INTGI 16  I 
C 

C  IF  ANY  ERRORS  FLAGS  THE  DATA  OF  THE  SATELLITE 
C 

00  120  INDEX* l.NSAT 

C 

C  ERRORS  IN  SIGNAL 

r» 

v 

IF!  MQVELI  INDEX). EQ.O)  GOTO  114 
MQVELI INDEX  1  =  15 
IH8+IN0EX 
I2=12+IN0EX 

WRITE!  6,110  I  MONTH, DAY  .YEAR  .HOUR  .MIN.SEC  , INDEX, INTGI II I  .INTGl  12  I 
110  FORMAT!/, 1  **-**-*»  WARNING  *#-*#-##  AT  ’, 

9  2112, IX), 14, IX, 2112, ),F5.2,/, 

*  '  DATA  FROM  TRACER  • , 

9  12, •  GESAR  BAO  STATUS  ',15, IX, 15) 

C 

MQVELI INDEX  1*15 

C 

C  IF  THE  DATA  IS  ALREADY  BAO  NO  NEED  TO  TEST  THE  SATELLITE 
C 

GOTO  120 
C 

C  NON  EXISTENT  SATE LITE  DATA  FOR  THAT  TRA'SR 
C 

114  CONTINUE 

IF!  ST  AT ( INDEX  I . EQ . ISAT ( INDEX  1 1  GO  TO  120 
C 
C 

WRITE (6,115)  MONTH , 04Y , YE AR , HOUR , MIN , SEC . INDEX , 

#  ISAT  I  INDEX), ST  AT  I  INDEX) 

115  FORMAT  </,  '**««*  WARNING  *«  »  »*  AT  ’  , 


*%*!>  v?  jomcvcwc >  1'  < 


nmmammwMUUM'U™*™ . . 


120 


ISO 

C 

C 

C 


*  2X, 12. IX, 12,  IX,  14, 2X, 12, ' : ' ,12. ' : ' ,F5.2,/, 

*  •  TRACKING  TRACER  1  ,12.  '  SV  ‘,12, 

«•  EXISTING  NAV  DATA  FOR  SV  •,I2,'  FLAGED  BAD  QUALITY  DATA'  1 
MQVELI  INDEX  U15 
CONTINUE 

COMPUTS  THE  SIGMAS 

00  150  IN0EX=1,NSAT 

IFI  MOVE L < INDEX  ).NE. 01  GOTO  150 

CALCULATE  LI  ANO  L2  SIGMA  RANGE 

BLC=TLDLL 

BP0C=PD8PC 

SGRK  INDEX  )  =  (  (  BLC*M  l/<  10 .  **(  CN1(  INOEX  )/10 .  1>  > 

9  »  1  .S*BPDC/(  10.  »*(CN1(  INOEX  1/10.  )))*K1+K2 

SGR21  INOEX  )*(  <  BLC»M )/( 10 .  *«l  CN2(  INDEX  )/10.  )  I  ) 

9  *  I .5»BPDC/( 10 . **( CN2I INDEX )/10.  )  )  )*K1*K2 

CALCULATE  LI  ANO  L2  SIGMA  OOPPLER  OATA 

BLC=TLPLL 

BPOC=POBC 

SG01(  INDEX  )  =  (  (BLC*M)/(  10. **( CN1( INDEX )/10.  ) )  ) 

9  *  11.  ♦(  ,5*BP0C/<  10.  **(CN1<  INOEX  1/10.  >)>  )»K31*K4 

SG021  INOEX  )=(  (  8LC#M  )/l  10.**<CN2<  INDEX  1/10.  )  )  ) 

*  »  ( l.+t  .5»BP0C/<  10.**(CN2(  IN0EX)/10.  ))  )  >»K32*K4 
CONTINUE 

OUTPUT  THE  DATA  TO  THE  FILE  IF  THERE  IS  NAV  OATA 
TYPE =3 

WRITE! 7)  TYPE 

WRITE17)  TT,  ISAT,  CR1,  CR2,  OOP1,  D0P2, 

9  SGR1,  SGR2,  SG01,  SGD2,  HOVEL 

GO  TO  READ  ANOTHER  BLOCK 

GOTO  10 


100 


CONTINUE 

IF  ( BLOCK. NE. 104)  GOTO  200 


C*» 
Ch 

c 


NAVIGATION  OATA  (AS  TRANSMITEO ) 

TRACER=INTG( 1 ) 

SAT  sINTGl 2 ) 

GO  TO  READ  ANOTHER  BLOCK 
GOTO  10 


200 

C 
C* 

C«H* 

c« 
c 
c 
c 

c»» 

c 

c 

c 


CONTINUE 

IF  1  BLOCK . NE . 9 )  GOTO  300 


NAVIGATION  OATA  (DECODED  SUBFRAMES  1  TO  3  I 

DAY  OF  KEEK 

IW<NOS=OIWT<  FPK  4 )) 

IH<NOS= I  AND  ( IWKNOS  ,17778) 

SV  HEALTH 

IEDATS*OINT(  FPI(  9)) 


89 


w-V. 


nwwnro ■  wwwi»w»JIIJ1 


mcrrm  r***w»w*  - 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CLOCK  EPOCH  ( GPS  SECONDS  OF  WEEK  ) 

T0CS*FPIU3) 

CLOCK  MAS  (SEC> 

ASI1)*FPKX6) 

CLOCK  DRIFT  < SEC/SEC ) 

ASI  2  )*FPII  IS  ) 

CLOCK  DRIFT  RATE  <SEC/SEC**2> 

AS(  3 )=FPIt 14  ) 

AGE  OF  DATA  (CLOCK)  (SEC) 

A0C*FPIt 10) 

RADIAL  SINE  CORRECTION  (DIVIDED  BY  1000.  TO  GET  KM) 
CRSS*FPI1 271/1000. 

CORRECTION  TO  MEAN  MOTION  1 RAOIANS/SECONOS ( 

0NS*FPI( 28) 

MEAN  ANOMALY  AT  EPOCH  (RADIANS) 

M0S=FP1<  29) 

IN  TRACK  COSINE  AMPLITUDE  (RADIANS) 

CUCS=FPI( 30 ) 

ECCENTRICITY 
ES*FPIl XI ) 

IN  TRACK  SINE  AMPLITUDE  ( RADIANS ) 

CUSS*FPI(  32 ) 

SQUARE  ROOT  OF  SEMI -MAJOR  AXIS  CONVERTED  TO  SQR  OF  KM 

SQAS*FPI (  J3  )»FPI (  3 J  ) 

SQAS*SGAS/1000 . DO 
SQAS=0$QRTl SQAS ) 

TIME  Of  EPOCH  (GPS  SECONDS  OF  WEEK) 

T0ES«FPII34) 

INCLINATION  COSINE  CORRECTION  (RADIANS) 

CICS«FPI ( 46 ) 

RIGHT  ASCENSION  NODE  ( RADIANS  ) 

0M€GS*FPIt  47 ) 

INCLINATION  SINE  CORRECTION  (RADIANS) 

CISS*FPI( 48  ) 

INCLINATION  (RADIANS) 

IOSsFPI ( 44  ) 
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L 


Wi  VL'Al  Vi 


.  -,\a\  Km 


KIT*  w»»rw  »  " 


C  RADIAL  COSINE  ADJUSTMENT  (DIVIDED  8Y  1QOO  TO  GET  KM  ) 

C 

CRCS*FPI( 50 1/1000. 

C 

C  ARGUMENT  OF  PERIGEE  (RADIANS) 

C 

WSsfPK  51 ) 

C 

C  RIGHT  ASCENSION  OF  ASCENDING  NOOE  (TIME  DERIVATIVE) 

C  I  RADIANS/SEC ) 

OM£QSsFPI  (  52 ) 

C 

C  AGE  OF  DATA  (SEC) 

C 

A0E*FPI<26) 

C 

C  INCLINATION  TIME  DERIVATIVE  I  RADIANS/SECONDS  ) 

C 

IDOTS=FPI( 54 ) 

C 

C  SETS  THE  FLAG  FOR  EXIST ANCE  OF  NAV  DATA  FOR  THIS  SATELITTE 
C 

STATI TRACER )*SAT 
C 

C  TESTS  IF  THERE  IS  ENOUGH  NAV  DATA 
C 

SAVE.. TRUE. 

00  210  IN0EX=1,NSAT 

IFI ST AT  1  INDEX ) . EQ . 0 )  SAVE* . FALSE . 

210  CONTINUE 
C 

C  OUTPUTS  THE  DATA 
C 


C 

c 

c 

c 


TYPE«2 

WRITE! 7)  TYPE 

WRITE) 7  I  TRACER, SAT ,AS( 1 ) , AS( 2  )  ,ASI 3  )  ,CICS,CISS 
#  .CRCS.CRSS, CUCS, CUSS,DNS,ES,IDOTS, IDS, MOS 

9  ,OMEDS,OMEGS,SGAS,AOE,  TOCS,  TOES,  IS*NOS,WS 

9  .IEOATS.AOC 

INFORMS  THAT  NAVIGATION  DATA  WAS  RECEIVED  FOR  THIS  SATELLITE 


WRITE ( 6 , 220  )  TRACER ,SAT .MONTH .DAY , YEAR , HOUR .MIN ,SEC 

9  ,  ’  **  NAVIGATIONAL  **',/, 

#  •  #*  DATA  FROM  TRACER  '  ,12, '  SV  ’,15,’  RECEIVED 

9  *  #»  •  ,7X,IZ,1X,IZ,1X,I4,2X,I2,’  ;  •  ,12,'  ;  •  ,F5. 2  ,7X> '  **>  ,/, 

PRINTS  OUT  THE  WARNING  WHEN  THE  HEALTH  IS  NOT  GOOD 

IF<  FPI(  9 )  .NE.  0.0  )  WRITEI6.201)  TRACER, SAT, FPK  9) 

TRACER  ’’,13,'  SV  M3,'  HEALTH  STATUS  1  ,F1S.6, 


GO  TO  READ  ANOTHER  BLOCK 
GOTO  10 


300  CONTINUE 

IF  ( BLOCK. NE. 101)  GOTO  400 


C*»  BLOCK  1  PLUS  INPUT  OATA 
C» 

C 

C  PRESSURE 
C 
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Lw. 


PB*FPI< 1 ) 


C 

C  TEMPERATURE 

C 

TPsFPIl 2  ) 

C 

C  HUMIDITY 

C 

HD=FPI( 5 ) 

C 

C  TESTE  FOR  SAO  HEATHER  OATA.  CHANGE  TO  STANOARO  VALUES  BAO  DATA 

C 

IF  II  PB  .LT.  800.  I  .OR.  (  PS  .GT.  1200.  ))  PB  *  980. 

IF  l(  TP  .LT.  -99.  )  .OR.  I  TP  .GT.  »99.  I)  TP  *  IS. 

IF  II  HD  .LT.  X.  )  .OR.  I  HD  .GT.  100.  )  )  HO  *  75. 

C 

C  OUTPUTS  THE  DATA  FOR  THE  FIRST  TIME 

C 

IF( .NOT. FIRST  )  GOTO  10 
C 

TYPE*1 

WRITE! 7  )  TYPE 

WRITE! 71  RL10FF  >RL20FF  »PB ,TP ,H0 .X , Y  ,2 
C 

C  SETS  THE  FLAG 

C 

FIRST*. FALSE. 

C 

C  GO  TO  READ  ANOTHER  BLOCK 

C 

GOTO  10 
C 

400  CONTINUE 

IF  ( BLOCK. NE . 3  )  GOTO  460 

C 


c 

C  SOLUTION  USER  EPOCH  TIME 

C 


TT*FPI ( 2  I 
C 

C  IF  NO  EPHMERIS  DATA  DOES  NOT  KEEP  THE  DATA 

C 

IF!  . NOT . SAVE  1  GOTO  10 
C 

C  COORDINATES  OF  RECEIVER 

C 

X=FPI(51 
Y*FPII 4  ) 

Z*FPU5> 

C 

C  COMPUTES  THE  DATA 

C 

CALL  SUB! TT. DAY, HOUR .MIN, SEC) 

DAY  *OAY  +STDAY 
C 

C  ONLY  PRINTS  POSITIONS  IF  UNCOMENTED 

C 

C  WRITE! 6 .450  )  MONTH. DAY .YEAR .HOUR ,MIN, SEC ,X ,Y ,Z 

450  FORMAT  (/,  *  GESAR  SOLUTION  AT  1 1*1  DAY  YYYY  HH :  FBf :  SS .  SSS  ) '  . 

9  2X.I2.1X.I2.1X.I4.2K.I2.' : ' .12.' :  ' »F5. 2 . 

«/,'  PSN  COORD  X,  Y,  Z  1KM )  '.FI4.J,'  ’.FI4.3,'  ’.F14.5) 

C 

C  CONVERTS  TO  LATITUDE  /  LATITUDE  /  HT 

C 

CALL  XYLLH!X,Y,Z,AG,ESG,OLAT  .OLOG.DHT.PI  ) 

C 

C  PRINTS  THE  RESULTS 

9: 


c 

C  WRITE! 6.455)  DLAT .DLOG.OHT 

455  FORMAT!  •  LAT ,  LON.  HT  (DEG,  DEG.  KM)  1  ,2F17. 12 ,2X,F12 . 7 ) 

C 

C  CONVERTS  LATITUDE  TO  DEGREES  HIFAJTES  AND  SECONDS 
C 

CALL  DMSIOLAT.LAOEG.LAMIN.LASEC) 

C 

C  CONVERTS  LONGITUOE  TO  DEGREES  MINUTES  AND  SECONOS 
C  IN  A  RANGE  0  TO  ISO 

IF  (OLOG.GT.1BO.OO)  0LOG=S60 .DO-OLOG 
CALL  OMS! OLOG,LOOEG,LOMlN, LOSEC  ) 

C 

C  WRITE  I  6 ,456  )  LAOEG  .LAJ1IN ,  LASEC  .LOOEG ,  LOMIN ,  LOSEC 

456  FORMAT!  '  LAT,  LON  < OEG  MIN  SEC,  DEG  MIN  SEG  )  '  ,2! 214 , IX ,F7 . 4 ) ,/ ) 
C 

C  OUTPUTS  THE  DATA  TO  A  FILE 
C 

WRITE  1 S  )  MONTH .DAY .YEAR .HOUR , MIN , SEC , 

8  LAOEG , LAMIN , LASEC  .LOOEG , LOMIN , LOSEC 

C 
C 

C  GO  TO  REAO  ANOTHER  BLOCK 

C 

GOTO  10 
C 

460  CONTINUE 

IF  ( BLOCK. NE. II )  GOTO  470 
C 


C 

WRITE! 6,461  1 


8  5X ,  '  **  TAPE  HEAOER/TRAILOR 

IF!  INTG!  I  l.EQ.l  )  WRITE  16,4621 
IF!  INTG  ll.EG.Z)  WRITE  (6,463) 

IF!  INTGl  1  I.EQ.4  I  WRITE  (6,4641 
IF!  INTG!  1 ) .  EG .  8  I  WRITE  16,465) 


462 

FORMAT! ' 

BEGINING  OF  DATA  SET' 

./! 

4o3 

FORMAT! • 

ENO  OF  DATA  SET' ,/ 1 

4o4 

FORMAT!  ' 

BEGINING  OF  CASSETTE' 

,/) 

4b5 

FORMAT!  • 

ENO  OF  CASSETTE ' ,/ ) 

WRITE! 6 

,466)  INTG!  2),  INTG!  I) 

466 

FORMAT! 1 

CASSETTE  ',15,'  IN  DATA  SET',/, 

8  •  CASSETTE  SEGUENCE  NUTBER  ',15) 

C 

C  GO  TO  REAO  ANOTHER  BLOCK 
C 

GOTO  10 
C 

470  CONTINUE 

IF  (BLOCK.NE.il)  GOTO  480 

C 

C  IF  NO  EPHMERIS  DATA  DOES  NOT  KEEP  THE  DATA  SO  NO  NEED  TO 
C  INFORM  THE  USER 

C 

IF! .NOT. SAVE )  GOTO  10 


C 


WRITE! 6,471  )  MONTH , DAY .YEAR .HOUR  , MIN , SEC 

8  '  •*  ERROR  BLOCK  RECEIVED  AFTER  :  ** ' , 

8  1  «*• ,9X,2( 12, IX  ),I4,2X,2I 12, ’ ; 1  ),F5. 2,8X,  '»*' ./ 


TOW  w  M.  ■ ji  wwmrrjTtrsr^ 


gvnwwp  w  ^  *"  HRWHUHATTU-v 


IV.IVn^DW-  -VlfW-.-V  ■*'. 


Li 


>■ ,/ I 
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WRITE  <  6  ,472  )l  IMTGI I )  .1*1  .NIX  ) 

FORMAT! '  ERROR  LOG  MESSAGE  TI  ICO  FORMAT  «  1' 


,18,/, 


a  2 
a  l 


ERROR  LOG  OVERFLOW  COUNT 


,18,/, 
.18,/, 
*  4’  ,18,/, 
«  S' ,18,/, 
, 18  ,  /  , 
, 18 ,/ , 
«  8‘  ,18./, 

' ,18,/  ) 


a  o 
a  7 


GO  TO  REAO  ANOTHER  BLOCK 
GOTO  10 


C 

440 


CONTINUE 

IF  (BLOCK. N£. 8)  GOTO  S00 


IF  NO  EPHMERIS  DATA  DOES  NOT  KEEP  THE  DATA  SO  NO  NEED  TO 
INFORM  THE  USER 


IF!  .NOT. SAVE  )  GOTO  10 


C 

C*»  TRACKING  CONFIGURATION 
C* 

C 

C  COMPUTES  THE  OATE 

C 


CALL  SUB! FPII 2  I  ,DAY .HOUR , MIN, SEC ! 
DAY*DAY*STOAY 


481 


>'  ,/, 


WRITE! 6 ,481  I  MONTH, OAY , YEAR, HOUR, MIN, SEC 
FORMAT!/,  ' 

a  •  **  TRACKING  CONFIGURATION  AT  **',/, 

a  •  **'  ,4X,2(  12,1X1 ,14,2X,2(  12,'  :  '  '  ,F5.2,8X,  *  ** '  ,/, 
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WRITE  16,4821  FPI!  1  l,(  FPII  I  ),!=!, NFI  ),!  IMTGI  I  ), 1*1, Nil  ) 

FORMAT!'  PSEUOORANGE  FTF  OF  VALIDITY  SEC  '.E22.14,/, 

5  '  PREDETECTION  BANDHITHS  HZ  FREOUENCY .TRACER ',/ , 

a  2I4E20. 14, /),//, 

8  '  TRACKING  LOOP  BANDHITHS  HZ  FREQUENCY .TRACER ’,/ , 

8  4!  4E20.14,/!,/, 

a  •  LOOP  ROUND  TRIP  CALIBRACION  DELAYS  IN  SECONDS  ’ ,/, 

INCLUDING  PRE-ANTEF8U  ‘  ,/, '  LI  '.E20.14,'  L2  ,E20.14,//, 
LOOP  ROUND  TRIP  CALIBRACION  DELAYS  IN  SECONDS  ',/, 

(INTERNAL  TO  RFM I ' ./ >  '  LI  ',E20.14,'  L2  ',£20. 14,//, 
SV  PRN  ID  OF  EACH  TRACER  ',/,4IlS,//, 

ANTENA  INOICATOR/TRACER ‘ ,/,4IlS,//, 

CODE  INDICATOR  FLAG/TRACER './ .4115 ,// ) 


440 


445 


TESTE  IF  THE  CONSTEIATION  MATCHES  THE  EXISTING  EPHMERIS 

DATA 

DO  445  IN0EX*1,NSAT 

IF( STAT! INDEX ) . EQ. INTG!  INDEX  I !  GO  TO  445 
WRITE! 6,440  I  INDEX, STAT! INOEX  )  ,INTGI INOEX  ) 

FORMAT (////,  1  <hhhhw»  WARNING  •***»»*  ******’,//, 

*  '  ACTUAL  NAVIGATIONAL  DATA  FOR  TRACER  ',12,'  IS  FROM  SV  ',12,/, 
a  '  TRACKING  CONFIGURATION  SHOWS  SV  ',12,//) 

CONTINUE 


C 

500 


GO  TO  REAO  ANOTHER  BLOCK 
GOTO  10 
CONTINUE 


44 


-1 


> 

% 


. . 


C»*  SLOCKS  THAT  THIS  PROGRAM  IS  NOT  PREPERED  TO  DECODE 


C 

501 


520 

550 

540 

550 

C 

c 

c 


WRITE! 5,501) 

FORMAT!//,' 

#  /,'  •* 

*  /,'  *» 

*  /, 

WRITE! 5,520)  8LK , SLOCK ,NFI , Nil , NCI 
IF  (NFI.NE.O)  WRITE  (  6  >550  )  (  FPU  I  )  ,1  =  1  ,NFI  ) 
IF  (NII.NE.O)  WRITE! 6, 540) 

IF  INCI.NE.O)  WRITE!  6,550) 

FORMAT!  1X,A5,4I5) 

FORMAT!  IX, 4E 20. 141 
FORMAT! IX, 6112) 

FORMAT! IX, 10 I A8 )) 

GO  TO  READ  ANOTHER  BLOCK 

GOTO  10 


FOUND  A  BLOCK  THAT  THIS  PROGRAM  IS  NOT 
ABLE  TO  DECODE  .  BLOCK  DATA  IS  DUMPED  NEXT 


!  1NTGI I )  ,1  =  1  >NII  ) 
!  CHI! I  ), 1*1, NCI ) 


*#  ’ 
**  • 


C 

C 

C 

C 

2000 


END  OF  PROGRAM 

CONTINUE 

STOP 

END 


THIS  SUBROUTINE  WAS  TAKEN  FROM  THE  KALMN  PROGRAM 


SUBROUTINE  XYLLH 

G  (  X  ,  Y  ,  Z  ,  AG  ,  ESQ  , 

R  DLAT  ,  DLOG  ,  DHT.PI) 

C  THIS  SUBROUTINE  CONVERTS  POSITION  FROM  THE  EARTH-CENTERED  CARTESIAN 
C  COORDINATE  SYSTEM  TO  THE  GEODETIC  COORDINATE  SYSTEM. 

C  X,Y,Z=ESTIMATE  OF  RECEIVER  POSITION  IN  KM 
C  AG-SEMIMAJOR  AXIS  OF  REF.  ELLIPSOID! RADIUS  OF  EARTH)  IN  KM 
C  ESQsECCENTRICITY  SQUARED  OF  REF.  ELLIPSOID3!  2 . -1 . /OBU/OBL 
C  DLAT  3GEOOETIC  LATITUDE  OF  RECEIVER  POSITION  IN  DEGREES 
C  DLOG=GEOOETIC  LONGITUDE  OF  RECEIVER  POSITION  IN  DEGREES 
C  DHT sGEODETIC  HEIGHT  OF  RECEIVER  POSITION  IN  METERS 
C  PI  VALUE  OF  PI  COMPUTEO  IN  DOUBLE  PRECISION  IN  THE  MAIN  PROGRAM 
C  SINCE  LATsLAT! X»Y ,Z,LAT  )  AND  HT=HT< X , Y , LAT  ) ,  AN  ITERATIVE  PROCEDURE 
C  IS  NECESSARY  TO  DETERMINE  THE  VALUES  OF  LATITUDE  ANO  HT. 

C  THE  WORKING  EQUATIONS  ARE  AS  FOLLOWS: 

C  HT  =  ( R/COSI LAT  )  l-AGl  WHERE  AG1=AG/! SQRT! 1-ESQ»SIN( LAT  )*»Z  )  ) 

C  TAN! LOG  1=Y/X 

C  TAN! LAT  )  =  ( ( Z+<  ESQ*SIN( LAT  )*AG1 )/R  ) 

C 

DOUBLE  PRECISION  LT  ,  LT1  ,  LG  ,  R  ,  EP  ,AG1 
DOUBLE  PRECISION  X  ,  Y  ,  Z  ,  AG  ,  ESQ  ,  DLAT  ,  OLOG  ,  OHT  .  t 
C 

C  FIND  LONGITUDE 

LG  *  DATAN21Y.X) 

IF  (LG  .LT.  0.00)  LG  3  LG  ♦  (2.  *  PI) 

R=OS««T( X**2+Y»*2 ) 


A  FIRST  GUESS  FOR  LATITUDE  WOULD  BE  ATANI Z/R  ) 
LT1=DATAN(Z/R  ) 

SOLVE  FOR  LATITUDE  BY  ITERATIONS 
DO  10  in. 5 

AG1-AG/I  OSQRT!  1 .  DO-ESQ*OSIVf  iTl  «; 
LT=DATAN1  I  Z*l  ESQ»OS:ni  L'i  *451  » 

EP=OABSl LT-LT1  I 

IF!  EP.  LE  .  .  OOOCOOOO  100  <  GCTC  >> 


in  «  lt 

10  CONTINUE 

5  AGX-AG/I OSGRTI 100-ESQ»0SIN< LT1 )**2 ) ) 

C 

C  FIND  HEIGHT 

OHT*( ( R/DCOSI LT ) l-AGl ) 

C  CONVERT  OUT  ANO  OLOG  FROM  RADIANS  TO  OEGREES 
DLAT«LT»100. 00/PI 
olog-lg*X8o . 00/PI 
RETURN 
END 
C 

CIHHMHHHHHHHHMHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHMHHHWMHHHHHHt 

C  THIS  SUBROUTINE  CONVERTS  DECIMAL  OEGREES  IN  OEGREES  MINUTES  SECONDS 

. . . 

c 

SUBROUTINE  OMStOEC , LDEG, MIN, SEC ) 

C 

DOUBLE  PRECISION  DEC  ,»AJM  ,XMIN 
C 

XNUM-OASSC DEC  I 
LDEG-OINT(XNUM) 

XMIN-I XNUM-OFLOATl LDEG  )  1*60 .00 
MXN-OINK XMIN ) 

XNUM-I XMIN-DFLOAT(  MIN ) 1*60 . DO 
SEC-SNGLI XNUM ) 

IF  ( LDEG. GE. 560)  LOEG-LDEG-560 
RETURN 
END 
C 

C  ********************************************************************** 
C  THIS  SUBROUTINE  CONVERTS  TIME  TAG  IN  OATS  HOURS  MINUTES  SECONOS 


C 

SUBROUTINE  SUBI TT ,OAY, HOUR ,MIN, SEC  ) 

INTEGER  OAY  >  HOUR,  MIN 
OOUBLE  PRECISION  TIME  >TT 
REAL  SEC 
T1ME-TT 

OAY  *  DINT) TIME  /  (26.00*5600.00)) 

TIME  *  TIME  -  DFLOATI DAY )  *  26. DO  •  5600.00 
HOUR  ■  OINK  TIME  /  5600.00) 

TIME  *  TIME  -  DFLOATI HOUR )  *  5600.00 
MIN  s  OINK  TIME  /  60.00) 

TIME  *  TIME  -  DFLOAT ( MIN )  *  60.00 
SEC  *  SNGL( TIME ) 

RETURN 

END 

/» 

//GO . FT06F00I  00  SYSOUT«* 

//GO. FT07F001  00  OSN-MSS. 30812. CVFICA. OUTPUT, OISP-SHR 
//GO. FT08F001  00  OSN-MSS. S0S12.GESAR.P0S.QISP-SHR 
//GO. FT5XF001  00  DSN-MSS. S0B12.CV0ATA, OISP-SHR 
//GO. FT50FOOI  00  OSN-MSS. S08X2. GPS. TAPES . DATA! 0L28S  )  ,OISP«SHR 
//  00  OSN-MSS. S0B12. GPS. TAPES. 0ATAIDL286), OISP-SHR 

//  00  OSN-MSS. S0812. GPS. TAPES ■ OATAI 0L286  >, OISP-SHR 

//  00  OSN-MSS. SC812. GPS. TAPES. DATAIDL287), OISP-SHR 

//  00  DSN-MSS. S0812. GPS. TAPES. DATA! DL288 ), OISP-SHR 

//  00  DSN-MSS. S08X2. GPS. TAPES. DATA! DL286 >, OISP-SHR 

//  00  0SN-MSS.SO8X2. GPS. TAPES. DATAl 0L29O ), OISP-SHR 

//  00  0SN-MSS.SO8X2. GPS. TAPES. DATA! 0L2H ) .OISP-SHR 

//  00  DSN-MSS. S0812. GPS. TAPES. OATAI 0L262 >, OISP-SHR 

// 
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APPENDIX  F 

PROGRAM  TRANSDUC  GPS.  SOURCE  LISTING 


//J08TRANS  JOB  (  0812 , 9999 )  > '  EZEGUIEL 1  .CLASSIC 
//•MAIN  0«G»NPGVM1.0812P,  LINES*  (  99 )  .CAROS*  1  99 ( 

//•FORMAT  PR  .DDNAME  *G0 ■ FT06F001 > 

//•FORMS «SEP1 

//  EXEC  F0BTVCL8 

//FORT . SYSIN  00  • 

C 

C  PROGRAM  TRANSOUC  GPS 

C 

C  AUTHOR:  AUGUSTO  EZEQUXEL 

C  DATE:  MARCH  30 ,  198? 

C 

C  DESCRIPTION: 

C 

C  THIS  PROGRAMS  TAXES  EACH  POSITION  OF  THE  ANT  EMU, 

C  THE  COMPUTED  COURSE  <  FILTERED  OR  NOT  )  .THE  PITCH 

C  ANO  ROLL  DATA,  AND  COMPUTES  THE  POSITIONS  OF  THE 
C  TRANSDUCER,  APLYXN6  THE  SEVEN  PARAMETER  TRANSFORMATION 
C  TO  THE  OFFSETS  OF  THE  TRANSDUCER  IN  RELATION  TO  THE  AMTEM4A 

C  CONSIDERED  THE  CENTER  OF  COORDINATE  SYSTEM  AND  THEN  AOOXNB 

C  THE  CORRECTIONS  TO  THE  COORDINATES. 

C 

C  THE  PROGRAM  RUNS  IN  MVS 

C 

C  I/O  SPECIFICATIONS:  SEE  END  OF  THIS  JOG. 

C 

C  THIS  PROGRAM  IS  RESTRICT  AS  IS  TO  THE  SEAFLOOR  BENCHMARK 
C  EXPERIMENT 

C 

C  THE  PROGRAM  MILL  RUN  NITH  ANY  AMOUNT  OF  DATA  ONLY  LXHITEO  TO  DISK 
C  SPACE 

C 

C  ANY  BLANC  LINES  HILL  TERMINATE  THE  PROGRAM  IN  ERROR 

C 

c 

c 

DOUBLE  PRECISION  XPOS ,YPOS,TINEl , TIME Z, TIMES, SECD,SfCH, RATE 

DOUBLE  PRECISION  TIMER,  TIMES 

DOUBLE  PRECISION  TT,OLAT,OLOG,OHT,X,Y,Z 

REAL  COURSE .PITCH! 2 )  ,ROLL( 2 )  ,SEC1<  2  I  ,SEC2,PIT ,R0L,SEC31 2  ( 

REAL  OX.DY .DZ.OFFX.OFFY.OFFZ.HEAOI 2 ( ,LASE6, LOSES 
INTEGER  M0NTH1!  2 )  ,0AY1(  2 ( ,  YEAR11 2 )  ,H0UR1(  2  )  .MINK  2 1 
INTEGER  M0NTH2  ,DAY2  >YEAR2  .HOUR 2  ,MIN2 
INTEGER  MONTHS! 2  1 ,0AY3I 2  ) , YEARS! 2 ( .HOURS! 2 ) .MINS!  2 ) 

INTEGER  LA0E6 » LAMXN , L00E6 , LQMXN ,  STDAY  ,  NSAT 
CHARACTERS  TITLE!  10 1 
C 

C  INITIALIZATION  OF  CONSTANTS 

C 

SCCH«3*00.00 
SECO>M .  00*SECH 

C 

0FFX--4.2M 
0FFY»-2. 9SS 
0FFZ»-9.502 
C 

C  READS  THE  TITLE  FROM  THE  FILE 

C 

READ!  SO  ,5 ,EN0*100  (TITLE 
HRITE! 9,5 )  TITLE 
READ! 30,5,END*100 (TITLE 


non  non  onn  onno  n  non  non  n  n  n  non  n  non  o  nor* 


MUTE!  4. 5)  TITLE 
S  FORMAT!  10AS ) 

RSA0iM.*,EM«100>  MNTHS.STOAY  .YEARS 
MRXTEI 4,4)  MONTHS, STOAY, YEAR 
4  FORMAT!///. 

•  '  OATS  Of  STARTZNB  CPS  MEEK  (MONTH  DAY  YEAR)  ' .IS. IS. IS  I 
READS  TNO  SETS  OP  PITCH  ANO  ROLL  DATA 

READ!  51  .*>EN0*10O )  MONTHK 1 1 .0AY1(  1 ) ,  YEAR11 1 1  .MURK  1 )  .MINK  1 1 . 

•  SEC11 1 ) .PITCH!  1 ) .ROLL! I > 

TIME1*0FL0ATI  OAYK 1 1 )  *SE  CO  »0  FLOAT!  HOUR  If  1 )  >*SECM 

4  ♦  OFLOATI  MINK  11  )«40.  DO*ORLEI  SSCll  1)) 

READ!  SI  ,»  »IM«100 )  MONTHK  S )  .OAYli  1 ) ,  YE  ARK  S )  .MURK  S )  >MIN1(  S ) . 
4  SEC1I  SI.  PITCH!  S).  ROLL!  S) 

TIME 2*OFLOAT !  OAYK  1 1 ) •SECOND FLOAT  I  HOUR II  S  >  )«3ECH 

•  ♦  DFLOATIMXNK  S 1 )*40.00«OSLI(SCCK  Z  > ) 

READS  TNO  SETS  OP  COURSE 

READ! 52 ,SN0«100 )  MONTHS! 1 1  .OAVSI 1 )  .YEARS! 1 )  .HOURS! 1 ) .MENS!  1 ) . 

4  SECS!  1 )  .MAO!  1 ) 

TIME*»OFLOATI OAYSI 1 )  )«SECO«OrLOATI HOURS! 1 1  )«SECH 

•  ♦  OFLOAT!MZNS( 1 )  )H40.00*0SLIISRCSI 1 ) ) 

RE AO! 52 ,EN0«100 1  MONTHS! S I  .OAYSI S I  .YEARS! t  \  .HOURS! 1 ) .HXMSt 2 ) . 

4  SECS! 2  I  .NEAOI Z  > 

TIME5-0FL0AT! OAYSI  S )  )«SECO«OFLOAT( HOURS! S  >  )«SECH 
4  ♦  OFLOATIMZNSI Z  )  )«40.D0«0RLE( SECS! 2  1 1 

READS  ONE  POSITION  OP  THE  SHIP 

20  READ! 50 .EM-100)  TT  .  X  .Y.Z,  OUT  .  OLOR  .  DMT  .  MAT 

COMPUTES  THE  OATE 

CALL  SURITT.0AY2. HOURS .MINS .SECS) 

OAY2«OAY2*STOAY 

CONVERTS  LATITUOS  TO  DECREES  MINUTES  ANO  SECONOS 

CALL  ONE!  OUT .LAOER. LANIN, LASEC I 

CONVERTS  LOMCtTUOE  TO  OECREES  MINUTES  AMD  SECONOS 
IN  A  RANEE  0  TO  ISO 
IF  I OLOO . ST . IRQ . 00 )  OLOR ■ S40 . DO -OLOR 

CALL  OMSI OLOR. LOOCR.LOKZM, LOSER) 

MRITEI4.S4I  LAOEC.LAMIN, LASEC, LODES, LOHIN. LORES 
24  FORMAT!'  UT,  LON  IDES  MIN  SEC,  DEC  MI N  SEC)  '  ,21  SX4.1X.P7. 41 ,/ ) 

TIME5-0PL0ATI  DAYS  H*ECO«OFLOATt  HOURS  IMECH 
4  «  OPLOATIMSNS  H*40 .00«0SL(<  SEC2  I 

PQRZTtON  IN  TIMS  OETMBEM  THE  TNO  SITS  OF  PITCH  ANO  ROLL  OATA 
AMR  THE  TNO  SETS  OP  COURSE  OATA 
25  CONTINUE 

POSITION  TO  EARLY  FOR  PITCH  ANO  ROLL  OATA 

IF! TIMES. SI. TIMED  SO  TO  40 
MRITEI 4.50 )  MONTHS .DAYS .YEARS .HOURS .MINS .SECS 
SO  FORMAT!'  PSN  R  '  .IS. IX. IS. IX. 14, IX, IS, IX. IS. IX, F4.1, 

4  ‘  REJECTED.  M  PITCH  ANO  ROLL  OATA'  I 
SOTO  SO 

POSITION  TO  EARLY  FOR  PITCH  AM  ROLL  DATA 
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40  CONTINUC 

IFITIHiS.SC.TIW4)  SO  TO  *0 
HNXTIli.M)  MOUTHS »0AY2 » YCAN2 .HOUR! .MINI iSCCI 
so  pocnati1  psm  a  ',i2,ix,ii.ix,I4,ix,iz,ix,ii,ix,f4.i, 
t  •  MJKCTtD.  NO  COURSC  DATA* ) 

SOTO  CO 

POSITION  LATCH  THAN  THC  SCCONO  SCT  Of  PITCH  AND  NOLL  DATA 
40  CQMTXMJf 

ZflTXHCS.Lf.TZWC)  SOTO  70 

HOVfS  THC  PITCH  ANO  NOLL  OATA  TO  THI  PI NOT  SIT 

MONTH!)  1  >aMONTH!i  t ) 

0AVXIX»OAV1I2) 

YCANXI 1  )»VlARll S ) 

HOUNXI X  >“H0UR1!  2 1 
MZNl!X»MZNXI2> 

SICK  2  I “SIC 11  2 1 
PITCH! 1 )  «PITCH( 2 1 
NOLL! X ) “POLL I  2 ) 

TIW1«TIW2 

MAOS  A  SCCONO  SCT  Of  PITCH  ANO  NOLL  OATA 

NO AO I SX  ,p .IND“100  I  MONTH!! 2 1 .OAYXI 2 ) .YCANXI 2  > .HOUNXI 2 ) .MINI! 2 1 
•  SCCXI 2). PITCH!  2). NOLL! 2 1 

TIW2-OFLOATI 0AV1I  2  )  )«SCCD*OF  LOATl  HOUNXI  2  )  IPSCCH 
4  ♦  0PLOATIMINXI2)  >*40.00*OOLI(SCC1I  2  > ) 

SOTO  2S 


POSITION  LATIN  THAN  THI  SCCONO  SCT  Of  COURSC  OATA 
70  CONTINUI 

IfITIWS.Ll.TIWB)  SOTO  00 

MOVCS  TW  COURSC  OATA  TO  TW  FIRST  SIT 

MONTHS! X  ) -MONTHS I  2 ) 

DAYS!  X  l-OAYSI X  I 
YC ANSI  X  l-YCANSI X ) 

HOURS! X  I “HOURS!  2  I 
MINS! X  I -MINI!  2  I 
SICS!  2  I  “SICS!  2 1 
HCAO!  X  l-HCAO!  2  ) 

TIW4-TIWS 

MAOS  A  SCCONO  SCT  Of  PITCH  ANO  NOLL  OATA 

MAO!  S2  >CNO“XOO  I  MONTHS!  2 )  .DAYS!  2  >  .YCARSl  2 )  .HOURS!  2 1  .MZNSI 2 ) . 

•  SCCSI  *  I  .HCAO!  2  ) 

TIWO“Of  LOATl  SAYS!  2  I  MSCCO*Of  LOATl  HOURS!  2  )  )«SCCH 

•  •  Of LOATINSM!  2 1  )*40.00»OSL(I SCCSI  2 )  I 

SOTO  20 

POSITION  MTTHIN  TW  TMO  SITS  IN  TXW 

M  RATC“I  TTWS-TIW  X  )/l  TIW  2-TIHC 1 1 

PIT  “PITCH!  X  1*1  PITCH!  2  1-PITCHI  X  )  l“RAT« 

NOLWOLLI  X  )•!  NOLL!  2  I -NOLL!  X  I  l“RATC 
RATIal  TXWS-TIW4  I/I  TIWS-TIW4  I 
COURSC  >HIAO<  X  1*1  MAO!  2  l-HCAO  I  X  I  )“NATI 

CQMPUTS  TW  CORMCTZONS  ANO  THC  POSITZOW  Of  TW  TNANOOUCCN 

CALL  PTH! Of FH .Of FY .Of F2 .COURSC , PIT ,ROL .OX ,0V .02 1 


C0NPUT1S  TW  UTM  COONOZNATCS 


nnn  o  n  n  o  onoo  o n n n n  n  nnn  on  nnn  nnn  non 


CALL  flPUTNI  LA0C6.LAHZN.LASM.LOOC6.LOMZN.LOSE8.YPOS  .XPOS ) 

C 

xpob«xpos*ox 

VPOB«YPOS*OY 
CONVERTS  SACK  TO  6P 

CALL  UTHBPlYPOS.XPflS.LAOfe.LAMN.LASM.LOOM.LONIN.LOfiM) 
OUTPUTS  THC  RESULT 

HAITI ( 71  MONTHS ,0AYZ , YCAR2 .NOUNS  »MXNS .SIC 2 > 

•  LAOCe.LAMZN.LASK6.LOOC8.LOMZN.LOfiC8.NSAT 

HRZTI(t.RO)  MONTHS .OATS. TSARS. HOURS .MZNS. SICS. 

•  LAOCfi.  LAMZN.LASM.  LOOM.  LOMZN.  LOSS*  .NSAT 

TO  EORMATtStlX.ZS).lX.Z*.S(lX.ZS>.lX.E«.X,S(lX,IS.lX.IS.lX.P7.«l,ZS> 

RKAOS  ONE  MOM  POSZTZON 
80T0  SO 

CND  OP  PROORAM 

ZOO  CONTZNUC 
STOP 
CND 

SUBROUTZNK  PTH4H0PE .YOPE.SDPE .C0UR8S .PZTCH.ROLL .OK.OV.OZ  I 

THZS  SUBROUT ZNC  CONPUTIS  THC  CORMCTZONS  TO  THC 
COOROZNATSS  OUC  TO  THC  OPSCT  OP  ANTCNA  IN  RCLATZON 
TO  THC  TRANSOUCIR 

PZ*ARCOS( -1.0 1 
ROLL«ROLL*PI/XBO. 

PITCH-PITCMPZ/180. 

COURSE ■COURSE •PZ/100 . 

CHANBC  SZQNAL  OP  0RZ6ZNAL  OATA  TO  NATCH  THC  RZOHT  HAND 
SYSTEM  CONVICTION 

ROLL*-ROLL 

COSPHZ*  COSI ROLL) 

SINRHI*  SIN4  ROLL ) 

COSM  COS(  coursei 
SIM*  SIN(  COURSE  I 
COSM*  COBI PITCH  I 
SIM*  SZN4  PITCH  I 

OX  •  XOPP  •  COfiPNZ  •  CORK  • 

•  TOPE  •  <  COON  *  SZM  «  SIM  •  SZNPHZ  •  CORK  I  ♦ 

•  XOPP  •  I  SZM  m  SZM  -  COSH  •  SZNPHZ  •  COSR  ) 

OY  ■  HOPE  •  CQBPMZ  *  SZM  •  C -X.  »  ♦ 

•  TOPE  •  I  cam  m  CORK  •  SINN  *  SZNPHZ  •  SZM  I  ♦ 

I  2DPE  •  I  SZM  •  COSR  ♦  COSH  •  SZNPHZ  •  SZM  I 

OZ  ■  NOPE  •  SZNPHZ  - 
t  TOPE  *  I  SZM  ■  COSPHZ  )  * 

•  COPE  •  <  COSM  ■  COSPHZ  ) 

RETURN 
I  NO 


C 


SUBROUTINE  0PUTM4  LADCO .  LAMIN . LAS! 8 ,  LOOCfi . LQMZN , LOBES .NORTH .EAST  I 
OOUBLI  PRECISION  A.R.N.AR.SP.CP .OP ,EP .S.R1 .CSB.ISRP.RH.RP.KO 
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OOUiLl  W3CXSX0N  U .11  >M.IM  ,>t ,U ,U .SXNKC 
OOtOU  MCCISXON  0441, OLAH.NORTH. CAST. 0.044IMIN, PI, ION 

un  um.iflSM 

this  MMuriNi  camrrs  tni  uth  coomxnatcs  oo  or  in  nob  n 
xn  am  10  ccntoal  mmioxan  m  oo  oo  n 

CH**ltl.flOO 

PNXaOflOATI  LAMS  )*Of  LOATI IAHXN  )/*0 . 00*004.1 4  LASCO  1/1*00.00 
LONaOf LOAT  I LOOOO  4*07L0AT4 UWXN  >/*0 . 00 ♦004.1 1  LOS* 0  4/3*00 . 00 
LON* -LON 

OUtMal  L0N-O4  4*3*00 . 000 
C 

*■*370131.000 

Mm.ttto 

c 

KOaO. *00*00 

c 

•  a  **4  0-1.004/0 

C 

N  ■  (A-0I/4A«0) 

AO  •  A*4  4  1.0O-m*».BO/*.DO*IN— t-W**3  1*01 .00/**. 00*4 W**0  10*011 
BO  a  3  ■  00/» .  00*A*4  4  W  M**2  I  *7 . 00/0  ■  00*4  M**3  -W»0)»BO.  00/** .  00aN**0  4 
CO  ■  10. 00/1*. 00*A»4M— 2-W**3*3 .00/0.00*4 W— 0  M**»4 1 
00  a  30. 00/00. 00»A»4H*»3  WO*  1 11 . 00/1* ■  OOlOWiO  I 
CO  a  310.00/m. 00*A*IW— »  M**»4 
OMXHXN  a  0MS**0 . 00*1 . *0000<00***e-4 
C 

OXaOAOCOS4 -1.00 1 

c 

PHI  "0441/100. 00*01 

C 

3  a  AQwQHXMXN  00*00X10  *■  00*0*41 1  ♦coaoopw*.  00*0441 1 

•  •  DO*OOXNI  *  .00*0441  )*IO*OSXN(  *.00*0441  4 
01  a  440*3 

C 

SXNBCC  a  4 1.00/3*00.00 4/100. 00*0t 
SINOCCaOOXN4  SINSCC 4 

c 

CM  a  tAa*2-0**2)/A**C 
CSOO  a  ISO/4 1.00-IMI 

m  a  A* I  I.OO-IM 4/4  000RT4 1 . 00-f S0*0SXM4 ONI  4***  4  4**3 

C 

00  a  W1*l  1 .  PO*«333*OCOB4  P44X  4**«  I 

02  a  00*OSXN4  ONX  4*0C004  ONX  )*SXNB(C**2/t .  004*40*1 . 000 

OS  a  3XHB0C*M*B0  *03X111 OWI  >*0C00l ONX  4— 3/2*. 00*4 ». 00-OTAW4  044X  4—2 

•  *  * . 0O*CSQO*OCOO4 0441  )*»I«*.OQ*«3S3*CSBO*BCOS40NX 4**» 1*440*1.01* 

M  a  MO*0C004  ONX  4*OXM0CC*4(O*l .  0* 

03  a  SINBOC— 1*00*OC004  ONI  )*«l/« .  00*4 1 .  OO-OT AH4  ONX  l**2 

•  ♦  IS00*0C0B4  ONX  )**2  )*K0*1 .012 

c 

0  a  0OOlOO*OUN 
02  a  044*2 
03  ■  0**1 


A*  >  0**BP— C— **NO*00IN4  0NX  l*0C004  0NX  l**0/720.00 

•  •  4  *1.00-30. 00*OTAH4  ONX  144***0TANI  0441  H*** 

•  ♦  270  ■  00*3 300*0C0O4  ONX  )— I -330 . 00*0300*031144  ONI  4**2  4140*1 . 02* 
00  •  OMSXNOOC— S*00*tC0340NX  4— 3/1X0 .00*1 3. 00-10. B—0TAN4  ONI  4—2 

•  *  OTAPNONX  4***alO.0O*«300*0C00l0NI  (**« 

•  -  30.  OOf 300*03X444  ONX  4**»  4*440*1  ■  020 
N0*TNa01*«t*0t*03«0*aA* 

(AST ■) 0**OaOS*OS«00  4*300000 . 00 
MCTU3N 
■NO 
C 

SUOOOI/TXNC  4/TN004  NORTH ,  C AST .  LAOCO ,  LANIN .  LASIO .  LOOCO .  LOHIN .  LOUS  I 
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OOUBLE  PRECISION  A,R,N,5,  BP.CP,OP.EP,S,R1,ES8,ES8P,RM,RP,KO 
DOUBLE  PRECISION  R7,R8,R4, £5, 0,92, 03,04.05, 96, 06, AP,SZNSEC 
DOUBLE  PRECISION  RIO, 0LAM.N0RTH, EAST .RPHI.DPHI  .PHIHIN 
DOUBLE  PRECISION  EPRIME,0ELTA,MNUH,CH,PHI,LON,PI 
REAL  LASE6, LOSES 

THIS  SUBROUT: -C  CONPUTS  THE  UTM  COORDINATES  Of  OP  IN  NOS  72 
IN  ZONE  10  CENTRAL  HERIDIAN  123  00  00  M 

CH«-123.000 

*>6378135.000 

*>248.2600 

K0>0 . 444600 

PI«OARCOS( -1.00) 

•  •  **l  R-1.00  l/R 

N  >  (A-0I/IA«0) 

AP  ■  A*H1.00-N)»8.00/4.00«IN—2-M— 3 1*01. 00/64. 00*1 N**4  N**5I) 

M  >  3.00/2.00*4*1  IN-N**2l*7.D0/8.D©*(N**3-N**4)*88.00/64.0O«N**S) 
CP  •  15.00/16 . 0O*A*l N**2-N— 3*8 . BO/4 ■  00<M N*J4  N**8)l 
DP  •  35. 00/48. D0*A*IN**3-N**4*11. 00/16. DO*N**S> 

EP  •  315. 00/512 .00*** IN* *4  N— 5) 

fIRST  APROWMATION  OP  PHI 

PHI aN0RTH/30. 800/3600. 00 

COMPUT  TRUE  HERIOIONAL  OZST  AND  APROKZNATS  PHI 

00  100  I>1,8 

PMIMN»PHI*60 . 00*2 . 408882006660-4 
•PHI >PHZ/180 . 00*PI 

S  •  AP*PHIWN-0P*08INI  2  ■  00*RPHI  >*CP*08IN1 4  ■  BO*RPHI I 
•  -0P*0$IN(  6 . 00*RPHI  )*EP*OSXN<  0.04*RPHt  l 

ri  >  non 

OSLTA*NORTH-R1 
PMIX  OS LT A/30. 8/3600.  I*PHI 
100  CONTINUE 

RPHI>PHX/180.00*PI 


SIN8EC  •  1 1 . 00/360000  )/180 . 00*PI 
SIN8EC  •  OBXIH  SZNBIC  I 

C 

(SR  >  I A**2-0**2 l/A**f 
ESQP  «  ISO/ I 1.00-fSRI 

RN  >  6*1 1  ■  00-EBB  I/I  BOORTI 1 .  QB-EWB8IMI RPHI  )**2 » >— 3 

C 

RP  *  RW*I 1 . 00*E80P*8COBI  RPHI  )**2» 

R7>l  OTAfH  RPHI  )/<  2 .  BOUP— 2»*INB2C  I  J*t  1 . 00»ES0P*0COBI  RPHI  )**2  I 
•  *1 . D12/KO**2 

RB*I  DTA4H RPHI I/I  24 . 00*RP**4*SINSEC  I  >*i 5 . 00 ♦! . D0*0TANt  RPHI  l**2 
4  *6 .  80*ES8P*PC08I  RPHI  l**2  -6 . 00*«S8P*OSINI  RPHI  )**2  -3 . 00*ES8P— 2 

4  *00001  RPHI  X*4-4. 00*E80P**2*QC08I  RPHI  )** 2«0SINI  RPHI  )**2  ) 

4  *1.024/K0**4 

R4>1 .  OO/OCOBI  RPHI  I/I  RP*SZN8EC  1*1  ■  06/K0 

R10>1 . 00/0C08I RPHI  I/I  4 . 00*RP**3«»INB2C 1 . 00*2 . 00*OTAH! RPHI  l**2 
4  »ISBP*OCOBI  RPHI  l**2  1*1 . 018/KO**S 

C 

(PRIME *t AST -500000 . DO 
0  *  . 00000 100*CPRXHB 
02  • 

OS  •  0**S 
04  •  0**4 
05  •  0**» 


io: 


04  a 

D4NUM*94  *0T  AN! RPHI ) 

04-1  04NUH  )/!  720 . 00*RP**4*SINSEC  )*( 41 .00+90. D0*DTAN(  RPHZ  )**2 

•  +45 . 00*0 TAN  I RPHZ ) **4+107 . D0*ES9P*0C0SI RPHZ  )**2 

•  -1*2 . 00*CS9P*OSIN! RPHZ  )**2-4S . 00«ES9P*OTAN! RPHZ  )**2 

•  *OSZN(  RPHZ  )**2  1*1 . 03A/K0**4 

tKI «9*1. D0/DC0S!  RPHZ  >  >/<  120 . D0*RP**5*SINS£C  >*15.00+28.00 

•  *OTAN( RPHZ  )**2+24.00*0TANI RPHZ )*«4+4. 00*ESOP*0COSI RPHZ  )**2 

•  +8 . DO«ESQP«OSXNI RPHZ  >**2 >*1 . 050/KO**5 
C 

OPHZ«( -R 7*92+88*94-04  1/5400 . 00 
OLAM«<  R 9*9-810*95+85 1/3400 . DO 
PHZ*PHX«OPHX 
LON*CM+OUH 

CALL  OmiPHI, LA 08 6 . LANIN . LASES  > 

CALL  OMSI LON, L00E6.L0MXN. LOSES! 

RETURN 

END 

C 

SUBROUTINE  OMSI  OEC.LDBO.MXN, SEC) 

C 

DOUBLE  PRECISION  DEC  >XNUM  .W1IN 

C 

XNUM-OABS! DEC ) 

loeg«ointoooo 

XMIN-I  XNUM-DFLOAT!  LOES ) )*40 . DO 
MIN-OXNT! XNZN ) 

XNUH-I XMIN-OFLOAT! HZN >  >*40 . 00 
SEC+SNCLIXNUM) 

IF  I L0E6.6E . 540 1  L0EG-LDEG-340 
RETURN 
END 
C 

SUBROUTINE  SUB( TT, DAY. HOUR, HZN, SEC  > 

INTEGER  OAY,  HOUR.  HZN 
DOUBLE  PRECISION  TIME .TT 
REAL  SEC 
TiME*rr 

DAY  *  OXNTI  TIME  /  124.00*5400.0011 

TIME  >  TIME  -  0 FLOAT! OAY )  *  24.00  •  3400.00 

HOUR  «  0INT1TIME  /  3400.00) 

TIME  »  TIME  -  OFLOATIHOURI  *  3400.00 
MIN  «  DINT!  TIME  /  40.00) 

TIME  *  TIME  -  DFLOATI MZN )  «  40.00 
SEC  *  SNGLI TIME  ) 

RETURN 

END 

/* 

//GO. FT04F001  DO  SYSOUT** 

//GO. FT07F001  00  OSN*HSS.S0812.GPS.TRANSOUC.POS.OZSP*SHR 
//GO. FT50F001  00  OSN*MSS.SOR12.CVOATA,OISP»SHR 
//GO. FT50F001  00  OSN-MSS.SOB12.GPS.ANTEFMA.POS,DISP*SHR 
//GO. FT51F001  DO  DSN-MSS. S0812.PITR0LL .OISP-SHR 
//GO. FT52F001  00  OSN-MSS. S0812.S0AS. COURSE ,DISP»SHR 
// 
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APPENDIX  G 

PROGRAM  COMPARE  POSITION.  SOURCE  LISTING 


//JOBCOMPA  JOB  ( 0*12.9*99 )  > ' EZEGUIEL ' ,  CLASS *C 
//•MAIN  0RG*NPGVM1 . 0012P  .LINES*!  99 )  .CAROS* I  99  > 

//•FORMAT  PR .OONAME *00 . FT04F001 . 

//•DEST *NPGBACH 
//  EXEC  F0RTVCL8 

//FORT .SYSIN  00  • 

C 

C  PROGRAM  COMPARE  POSITION 

C 

C  AUTHOR:  AUGUSTO  EZEGUIEL 

C  DATE:  MARCH  27,  X9S7 

C 

C  DESCRIPTION: 

C 

C  THIS  PROGRAMS  TAKES  EACH  GPS  POSITION  ,  FINOS  THE 
C  CORRESPONDING  POSITION  USING  THE  MR  FALCON,  COMPUTES  THE 
C  X,  Y  ANO  RADIAL  DIFFERENCES  ANO  PRINTS  THE  DATE  TIME  TAG 
C  GEOGRAPHIC  POSITIONS  ANO  DIFFERENCES.  THIS  PROGRAM  READS 

C  THE  PERIOOS  TO  COMPARE  THE  DATA  FROM  A  SEPARATE  FILE. 

C 

C  THE  PROGRAM  RUNS  IN  MVS 

C 

C  I/O  SPECIFICATIONS:  SEE  ENO  OF  THIS  JOB. 

C 

C  THIS  PROGRAM  IS  RESTRICT  AS  IS  TO  THE  SEAFLOOR  BENCHMARK 
C  EXPERIMENT 

C 

C  THE  PROGRAM  HILL  RUN  HITH  ANY  AMOUNT  OF  DATA  ONLY  LIMITED  TO  OISK 
C  SPACE 

C 

C  ANY  BLANK  LINES  HILL  TERMINATE  THE  PROGRAM  IN  ERROR 

C 

C 

C 

DOUBLE  PRECISION  XPOS! 3 > .YPOSI 3 1.TIMEI 5>,SECD,SSCH,RATE 
DOUBLE  PRECISION  OIFX,OIFY,OIST,XNEH,YNEH 
REAL  SEC!  5 ) .LASEG!  3 )  ,LOSEGI 3 ) 

INTEGER  MONTH! 5 ) ,OAY( S ) , YEAR! 5 ) .HOUR! 5 ) .MINI S I ,NSAT 
INTEGER  LAOEG! 3 1 , LAMINI 3  ) , LOOEGI 3  ) , LOMIN! 3 1 , LINE , LOOP , INDEX 
C 

C  INITIALIZATION  OF  CONSTANTS 
C 

SECH*3400.D0 
SEC0*24. 00*S£CH 
OIFX-O.O 
DIFY*0 . 0 
01  ST  *0.0 
C 

C  IHVTS  THE  NUMER  OF  PERIOOS  TO  COMPARE  THE  DATA 

C 

REAOI 2S,*,ENO*100 )  LOOP 

C 

C  LOOP  OF  THE  PROGRAM  TO  THE  NUWER  OF  TIMES  DESIRED 
C 

DO  1000  INOEX  ■  1,  LOOP 
C 

LINE  *40 
REMIND t 30 ) 

RE  HI  NO!  SO  ) 

C 

READ! 2S ,* ,EN0*100  )  MONTH! 4  )  ,OAYI 4  I  .YEAR! 4 )  .HOUR! 4  I  .MINI  4  I  ,SECI 4  I 
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*  .MONTH!  5 1 ,0AY< 5 ) .YEAR! 5  J .HOUR! S )  .MINI  5 1  .SEC! 5 ) 

STARTING  TIME 

TIME! 4  )»OFLOAT( DAY! 4 ) )*SECD+0FL0AT<  HOUR! 4 1 )*SECH 
»  ♦  OFLOAT! MINI  4 )  )*60 . 00+DBLE!  SEC(  4 ) ) 

ENDING  TIME 

TIME!  S  )*OFLOAT(  OAYI  5 !  >*SECO*OFLOATI  HOUR!  5 )  )«SECH 
t  ♦  OFLOAT! MIN! 5 )  )*60. 00+D6LE! SEC!  5 ) ) 

READS  TWO  MR  FALCON  POSITIONS 

10  READ! SO >EN0*100  )  MONTH!  1 1  .DAY!  1 1,  YEAR!  1 )  .HOUR!  1 )  .MINI  1 )  ,SECI  1 ) . 

*  LADES! 1 )  .LAM1NI 1 ) .LASEGI 1 )  .L00E61 1 I  .LOMIN! 1 ) .LOSES!  1 ) 
TIME! 1  )*0FL0AT( DAY! 1 )  )*SEC0«-0FL0AT< HOUR! 1 )  >»SECH 

$  *  OFLOAT! MINI  1 )  1*60. 00+0SLE! SEC! 1 ) ) 

POSITIONS  THE  FIRST  POSITION  INSIOE  THE  PERIOO.  THE  REST,  OF  THE 
PROGRAM  MILL  ADJUST  SY  IT  SELF 

IF!  TIME 1 1 ) . LT . TIME ( 4 )  )  GO  TO  10 


READ! 30 ,EN0*100 )  MONTH! 2 ) ,0AYI 2  > .YEAR! 2 ) .HOUR!  2 1 .MINI  2 ) .SEC!  2  I , 

*  LADES! 2 ) .LAMXN! 2  ), LASEGI 2 ) .LOOEG! 2  )  .LOMIN! 2  )  .LOSES! 2 ) 
TIME! 2 )*DFLOATI DAY! 2  )  )»SECD+0FL0AT( HOUR! 2 ) )*SECH 

*  ♦  OFLOAT! MIN! 2  )  )*60.00+06L£( SEC! 2  ) ) 

COMPUTES  THE  UTM  COORDINATES 

CALL  GPUTM! LADES! 1 ) , LAMINI 1 1 .LASEGI 1 ) , LOOEG! 1 ) , LOMIN! 1 ) , LOSEGI 1 ) . 
«  YPOSI 1 )  ,XPOS( 1 ) ) 

CALL  GPUTM! LAOEGI 2 ) , LAMINI  2  I , LASEGI 2  ) . LOOEG! 2  ) , LOMIN! 2 ) , LOSEGI 2 ) , 

*  YPOSI 2  )  .XPOSI 2  ) ) 

READS  GPS  POSITION  OF  THE  SHIP 

20  REA0I50.EN0«100I  MONTH!  3  I ,OAY!  3 l.YEARI 3  l.MOURI 5  1  .MINI  3 > .SEC! 3 > , 

*  LAOEGI 3  1 , LAMINI 3 ) , LASEGI 3 ) , LOOEG! 3 ) . LOMINl 3  ) , LOSEGI  3 ) ,NSAT 
CALL  GPUTM! LA0E6I 3  ) , LAMINI 3 ) . LASEGI 3  ) . LOOEG! 3  ) , LOMIN! 3 ) ,  LOSEGI 3  > , 

$  YPOSI  3 ) .XPOSI 3  )  I 

TIME! 3  I'OFLOATI OAYI 3  )  )*SECO*DFLOAT ( HOUR! 3  I  l«SECH 

*  ♦  OFLOAT  I  MINI  3  1 1*40 .  D0*D6LE  ( SEC  IS)) 

POSITION  OUTSIOE  OF  THE  PERIOO 
PROGRAM  MILL  PICK  ANOTHER  PERIOO 

IF!  TIME! 3  I.GT . TIME! S )  1  GO  TO  1000 

POSITION  IN  TIME  U TWEEN  THE  TWO  MR  FALCON  POSITIONS 

25  CONTINUE 

POSITION  TO  EARLY  FOR  MR  FALCON  DATA 

IF! TIME! 3 1. GE. TIME!  1 ) I  00  TO  40 
»  WRITE!  *,30  1  MONTH!  3). OAYI  3), YEAR!  3), HOUR!  3), MINI  3), SEC!  3) 

»  30  FORMAT!'  GPS  PSN  9  ' ,12 .IX, 12 ,1X,I4, IX, 12, IX, 12 ,1X,F4. 1 , 

»  *  '  REJECTED.  NO  MR  FALCON  DATA'  ) 

SOTO  20 

POSITION  LATER  THAN  THE  SECONO  MR  FALCON  POSITION 
40  CONTINUE 

IF!  TIME!  3  I.LE.  TIME)  2))  GOTO  70 
MOVES  THE  SECONO  MR  FALCON  POSITION  TO  THE  FIRST  SET 
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MONTH! 1  )«M0NTH( 2 1 
DAY! X )*0AY( 2  ) 

YEARtl)«YEAR12> 

HOUR! 1  )*HOUR( 2 ) 

MINI  X  )«MINI  2 I 
SECI 2 )*SECI 2 I 
TIMEI X  )*TIM£(  2 ) 

LAOEGI 1  )«LA0EG<  2 ) 

LAMXNI X  )*LAMINI 2 ) 

LASEGI 1  )«LASEG(  2 ) 

YPOSI 1  )«YPOSI 2  ) 

LOOEGI X  I* LOOEGI 2 I 
LOMXNI 1  )*LOMIN< 2 ) 

L0SECI1  >»LOSEGI  2 ) 

XPOSI  X )>XPOSI 2  ) 

READS  A  SECOND  SET  OF  MR  FALCON  DATA 

REAOI 30 >ENO>100  I  MONTH! 2  )  .DAY! 2  )  .YEAR! 2  I  .HOUR! 2 ) .MINI  2 ) ,SEC< 2 I , 

*  LAOEGI 2  I  .LAMXNI 2 )  .LASEG1 2 ), LOOEGI 2 ) , LOMXNI 2 ) , LOSEGI 2  ) 
TIME! 2  l*OFLOATI 0AY1 2 )  IPSECO+DFLOATt HOUR!  2 ) )*SECH 

*  ♦  D FLOAT I  MINI  2  )  )*60 . 00+0BLEI SECI 2  )  ) 

CALL  GPUTMI LAOEGI 2 ) . LAMXNI 2 ) , LASEGI 2 ) .LOOEGI 2 ) .LOMXNI 2 ) .LOSEGI  2 ) . 

*  YPOSI 2).  XPOSI 2)) 

GOTO  25 

GPS  POSITION  NXTHIN  THE  TWO  MR  FALCON  POSITIONS 

70  RATE»I  TIME!  3  1-TIMEI  X  )  I/I  TIME!  2  l-TXMEI  1 1 1 

XNEHiXPOSi 1  )♦< XPOSI 2  I -XPOSI X  I ) PRATE 
YNEW*YPOSI 1  )♦( YPOSI 2 ) -YPOSI 1 1 IPRATE 

OIFX-XPOSI 3  l-XNEM 
OXFY«YPOSI3)-YNEM 
OIST-OSGRTI 0IFXpp2*0IFYpp2 ) 

OUTPUTS  THE  RESULT 


IF  ILINE.LT. *0)  00  TO  85 
WRITE  I  6.80) 

80  FORMAT) 1HX./.X7X. 'DIFFERENCES  BETWEEN  GPS  AND  MR  FALCON  POSITIONS' 

*  ,// ,  2W  > '  AVOIDING  SVXl'  , 

*  ,// ,4X , ' DATE ' , 6X , ' TIME ' , 09X , ' LATITUDE ' , 7X , ' LONGITUDE • » 

*  12X, ' I  METERS )' ,9X, '»',/, 2X,'M',2X,'D' ,XX, 'YEAR* ,1X, 'HH' , 

*  XX.'WI'.XX.'SS.S'.XX.'OO'.XX.'Ftr.XX.'SS.SSSS'.SX.'ODD'.lX. 

*  'MM' .XX.'SS.SSSS' ,5X,'0X' ,5X, 'DY • ,5X, 'DIST. ' , XX, 'SAT. • ,/) 

LINE«0 

85  LINE*LINE*1 

WRITE! 6 . 90  I  MONTH!  3  I ,OAYI  3 1 .YEAR!  3  )  .HOUR! 3  I  .MINI  3  I  .SECI 3  I , 

*  LAOEGI 3  ) , LAMINI 3  ) , LASEGI 3  I , LOOEGI 3  ) , LOMINI 3  I , LOSEGI  3  I . 

«  OIFX.OIFY.OIST.NSAT 

90  FORMATI 21 1X.I2 )  .1X.I4.1X  .I2.1X.I2  >1X,F4. 1 .1X.I3.XX.I2 .1X.F7.4. 

I  '  N', XX, 13, XX, 12, XX, F7. 4,'  W' ,1X,F6. 2 ,XX,F6. 2 ,XX,F7. 2.1X.I2 I 

READS  ONE  MORE  POSITION 
GOTO  20 

ENO  OF  PROGRAM 

XOOO  CONTINUE 
100  CONTINUE 
STOP 
ENO 


SUBROUTINE  GPUTMI LAOEG . LAMTN , L ASEG , LOOEG , LOMIN , LOSEO .NORTH .EAST ) 
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DOUBLE  PRECISION  A,R ,N,AP ,BP ,CP .DP ,EP ,S,R1 ,ESG,ESQP,RM,RP,KO 
DOUBLE  PRECISION  R2,R3 ,R4,R5,P,P2 ,P3 ,P4,P5,P6 ,A6 >BS, SINSEC 
DOUBLE  PRECISION  PHI.OLAM, NORTH, EAST .B.PHIMIN, PI, LON 
INTEGER  LAOEG, LANIN .LOOEG.LOMIN 
REAL  LASEG. LOSES 

THIS  SUBROUTINE  COMPUTS  THE  UTM  COORDINATES  OF  GP  IN  MGS  72 
IN  ZONE  10  CENTRAL  MERIDIAN  123  00  00  M 

CM*-123.0D0 

PHI=QFLOAT( LADES  )«OFLOAT( LAMIN >/60 . OO+OBLE ( LASEG  1/3600 . DO 
LON^O FLOAT! LODEG  ) +0FLOAT ( LOMIN  1/60 . OO+DBLE ( LOSEG  1/3600 . 00 
LON=-LON 

OLAM=(  LON-CM  1*3600 . 000 
C 

6=6378135.000 

R*298.2600 

C 

K0*0. 999600 
C 

B  *  A*( R-1.00  )/R 
C 

N  *  1  A-B 1/1  A*B  1 

AP  a  A*( ( 1.00-N1+5. 00/4. D0*(N**2-N**3H81. 00/64. DO*t N**4-N**5 1  > 

BP  a  3. 00/2. 00*A*((N-N**2>*7. DO/8. 00*1 N**3-N**4  1+55. 00/64. D0*N**5  1 
CP  *  15.00/16 .DO*A*<  N**2-N**5+5. 00/4. D0*( N**4-N**5 1 1 
OP  «  35. 00/48. 00*A*(N**3-N**4+11. 00/16. D0*N**51 
EP  *  3 15 . 00/512 . DO*A* I N**4-N**5 1 
PHIMIN  -  PHI*60.00*2.90888208666D-4 
C 

P1=DARC0S( -1.D01 
C 

PHI =PHI/180 . 00*PI 
C 

S  ■  AP*PHIMIN-BP*OSIN<  2 . 00*PHI 1 ♦CP*OSIN( 4 . D0*PHI 1 

*  -  0P*0SIN(6.00*PHIl+EP*OSIN(8.O0*PHI) 

R1  *  KO*S 

C 

SINSEC  *  (1.00/3600.00 1/180. 00»PI 
SINSEC=DSIN(  SINSEC 1 
C 

ESQ  =  ( A**2-B**2  )/A**2 
ESQP  =  ESQ/t 1.00-ESQ  1 

RM  =  A*( 1. OO-ESQ  )/( DSQRTI 1 . DO-ESQ*DSIN( PHI )**2 1 )**3 
C 

RP  a  RH*I 1 . DO+ESQP*DCOS( PHI )**2 1 

R2  =  RP*0SIN(PHI)*0CaS(PHIl*SINSEC»*2/2.D0*KO*1.0D8 

R3  a  SINSEC**4*RP*0SIN( PHI  1*OCOS( PHI  l**3/24.00*l 5. DO-DTAN! PHI  1**2 

*  ♦  9.00*ESQP*DCOS(PHI  1**2*4.00*ES<8P*ESQP*OCOS(  PHI )**4 1*K0*1 . 016 
R4  *  RP*OCOS(  PHI  )*SINSEC*K0*1 . 04 

R5  a  SINSEC**3*RP*OCOS(  PHI l**3/6 .00*1 1 .OO-OTAN!  PHI )**2 
$  ♦  ESQP*DCO$< PHI  1**2  )*K0*1 .012 

C 

P  *  .000100*0 LAN 
P2  «  P**2 
P3  ■  P**3 
P4  a  RM4 
P5  ■  P**S 
P6  a  p**6 
C 

A6  *  P6*SINSEC**6*RP*0SIN( PHI )*DCOS(  PHI )**5/720.00 

*  •  (61.00-58.00*07 AN! PHI )**2»0TAN( PHI 1**4 

»  ♦  270. DO*ES9P*OCOS( PHI )**2-330. DO*ESQP*OSIN( PHI )**2 )*K0*1 . 024 

B5  >  P5*SINSEC**5*RP*DC0S( PHI  1**5/120. D0*( 5.00-18. 00*OTAN( PHI )**2 
»  ♦  DTANIPHI  l**4»14.OO*ESQP*0COS(  PHI )**2 

$  -  58.00*ESGP*0SIN(PHI  )**2  l»KO*l.D20 

N0RTH«R1+R2»P2*R3»P4+A6 
EASTs( R4*P*R5*P3*B5 1+500000.00 
RETURN 
END 
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END 


//GO. FT06F001  DO  SYSOUT»» 

00  °SN«HSS.S0«12.PERI00S. COMPARE, DISP*SHB 

;~3SS  SS  SS«:SB£3SS3BSr 


[ 

I 

I 
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APPENDIX  H 

PROGRAM  COMPARE  PLOT.  SOURCE  LISTING 


//EZEGUIEL  JOB  I  0812,9999), 'EZEGUIEL 1  >CLAS5*C 

//  EXEC  FRTVCLGP 

//FORT.SYSIN  00  * 

PROGRAM  COMPARE  PLOT 

RUNS  IN  FORTRAN  VS 

AUTHOR:  AUGUSTO  EZEGUXEL 

DATE  :  31  MARCH  1987 

THIS  PROGRAM  MAXES  A  PLOT  OF  THE  POSITIONS  OF  THE  SHIP 
USING  GPS  AND  MR  FALCON  POSITIONS  OF  THE  TRANSDUCER 


DOUBLE  PRECISION  XPOS.YPOS .XLEFT .YLEFT  > TIME! 4 1  ,SECD ,SECH 
REAL  S£C( 3  ) ,XPLT , YPLT >BLXH ,8LYH, SCALE , VALUE, LASECiLOSEG 
REAL  YORG,XORG 

INTEGER  MONTH! 3 ) , YEAR! 3  I  .DAY! 3  I  ,HOUR( 3  )  ,MIN( 3  ) 

INTEGER  IPEN , LAOEG  >  LAMIN , LODEG >  LOMXN 
INTEGER  10, LOOP, INDEX 

DIMENSIONS  OF  SHEET,  LEFT  CORNER  AND  SCALES 

BLXH*15. 

BLYH«20.S 
SCALE »l./5000. 

XORGaO.O 

YORG»-2S.O 

INITIALIZATION  OF  CONSTANTS 

TIME ( 4 )*0 . DO 
SECH*3600.00 
SECD*24.00*SECH 

PLOTTER  INITIALIZATION 

CALL  PLOTS  (0,0,0) 

READS  THE  AMOUNT  OF  PERIODS  TO  COMPARE 
READ! 30, *)  LOOP 
OO  100  INDEX»1,L00P 

READS  THE  PERIOOS  AND  INFERIOR  LEFT  CORNER  OF  EACH  AREA 

READ! 30,* )  MONTH! 1 ) ,DAY< 1 )  .YEAR! 1 )  .HOUR! 1 )  ,MXN( 1 ) ,SEC( 1 ) 

*  , MONTH!  2), DAY!  2), YEAR!  2), HOUR!  2), MINI  2), SEC!  2),  XLEFT, YLEFT 

STARTING  TIME 

TIME! 1  )«OFLOAT<  DAY! 1 )  )*SECD*DFLOATI HOUR! 1 ) )*SECH 
$  ♦  OFLOATI MIN!  1 ) )*60. DO+OBLE I  SEC! 1 )  ) 

ENDING  TIME 

TIME! 2  )*OFLOAT( DAY! 2  )  )*SECO»OFLOAT( HOUR! 2 ) )*SECH 

*  ♦  OFLOATI MINI  2  )  )*60. OO+OBLEI SEC! 2  ) ) 
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POSITIONS  THE  PEN  IN  THE  ORIGIN 
YORGa-YORG 

XORC*2S.O*(  l-(  -1  )<h*INDEX  )/2 
IF! INDEX. EG. 2  I  X0RG»0 
CALL  PLOT  (X0RG.Y0RG.-3) 

GRID 

CALL  PLOT  10.1, 0.1,3) 

CALL  PLOT  ( 0. l,BLYH-0. 1,2 ) 

CALL  PLOT  <  8LXH-0 . 1 ,BL YH-0 .1,2) 

CALL  PLOT  (BLXN-0. 1,0. 1,21 
CALL  PLOT  10.1,0.1,2) 

SCALE 

CALL  METERI SCALE) 

TITLE 

CALL  SYMBOL! OS. 1,19. 2 ,0 . 20,29HSEAFL00R  BENCHMARK  EXPERIMENT, 0. ,29) 
CALL  SYMBOLIOB. 1,18. B. 0.20, 29H  PHASE  II  ,0.,29) 

CALL  SYMBOL! 0B. 1,18. 4, 0.20, 29H  14  AUGUST  1484  ,0.,29) 

CALL  SYMBOLIOB.  1, IB. 0,0. 20, 29MAVOIDING  SVU.  NO  CONSTRAIN  ,0.,29> 
CALL  SYMBOL! 01. 0,03.375, 0.15, 3, 0.,-l) 

CALL  SYMBOL! 1.20,03.3,0. 15.29HMR  FALCON  TRANOUCER  POSITIONS, 0. ,29) 
CALL  SYtWJLI  01 . 0 ,03 . 075 ,0 . 15 ,11 ,0 .  ,-l ) 

CALL  SYMBOL! 1.20 ,03. 0,0. 15, 29HGPS  TRANOUCER  POSITIONS  ,0.,29) 

SETS  SYMBOL  FOR  MINI  RANGER  POSITIONS 
ISM*3 

REWINDS  THE  FILES 

REWIND! 31 1 
REWIND!  32 ) 

LOOP  FOR  THE  TYPES  OF  POSITION 
DO  80  IOa31,32 
SETS  PEN  UP 
IPEN-3 

LOOP  WITHIN  ONE  TYPE  OF  POSITION 
10  CONTINUE 
READ  THE  OATA 

READ! 10, END-50  ) MONTH! 3 ) ,0AY! 3  !  ,YEAR( 3  )  .HOUR!  3 ) , 

*  MINI  3 ) , SEC 1 3 ) , LA0E6 , LAMIN , L ASEG , LOOEG , LOMIN , LOSEG 

COMPUTES  THE  UTM  COORDINATES 

CALL  GPUTMI LAOEG , LAMIN ,  LASEG , LOOEG , LOMIN , LOSEG ,YPOS ,XPOS ) 

COMPUTES  THE  TIME  IN  SECS 

TIME! 3 )*OFLOAT( OAYI 3 )  )»SECD*OFLOAT( HOUR! 3 ) )*SECH 

*  ♦  OFLOAT1MIN1 3  )  1*40. 00+DBLEf  SEC! 3  ) ) 
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c 

c 

c 


TOTS  IE  NXTHIM  THE  PERIOO 


C 

C 

C 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


IHTWttS) 

iFfTxwm 


.er.Tweizn  aoro 

.IT.TIlfai)  GOTO 


SO 

10 


COMPUTES  THE  PLOTTER  COORDINATES 


XRtT*l  XPOS-XLEET  >*SCALE«100 

mr«(  ypos-yleet  i*scale*ioo  . 

TESTE  IE  INSIDE  AREA 


.tI - - 

IEIYPLT.LT.O. >  GOTO  SO 
lEJXPLT.GT.SUOU  GOTO  SO 
IEITPLT.GT.OLVH I  GOTO  SO 

TESTS  IE  POSITIONS  ARE  AMAY  MORE  THEN  10  SECS  IN  TINE 

IEITIMBI 3  J-TZHE(AI.ST.  10.00  I  IPCN.t 
TZMe<«)„rZMEtSI  001  inmt 

PLOTS  THE  POSITION 


CALL  PLOT( XPLT 
IPIN*2 


•VPLT.IPEN 


) 


PLOTS  THE  TINE  EVERY  MINUTE 


C 

C 

C 


C 

c 

c 

c 


lr{  I"X!  SEC  I  3  ) )  .GT.O  )  00  TO  10 

CALL  MPMRlXPLT«0.  lS,yPLT»0  .0  it  w.lt_  , 


*-i  i 


WN  UP  NHZLE  THE  POSITIONS  ARB  OUT  OP  THE  SHEET 

SO  IPEN>S 
SOTO  10 


END  OE  PLOT 
SO  CONTINUE 


C 

C 


C 


C 


ISN«1I 
I PEN*! 

00  CONTINUE 
100  CONTINUE 

CALL  PLOT»0.,0. 

STOP 

ENO 


**"OVTIHR  METERt  SCALE ) 


XO»*.0 

V0-1.0 

CAU  PLOTCHO,YO,-S) 


00  10  J«1,12 
XP»  ELOATIJ 1-1.0 
CALL  PLOTIXP.O.O.Sl 
CAU  PLOT<XP,0.2S,t  J 
10  CONTINUE 


DO  *0  J*1,0 
XP«ELOATl  JIPO.I 
CALL  PLOTIXP.O.O.Sl 
CALL  PLOT! XP.O. 2.2  I 


e 


20  CQMTZfM 


call  HONO.  .0.0,3) 

C ALL  PUmi.,Q.,|l 
CALL  KOTI0..0  2.J) 

CALL  OtOfNll.,2.21 
C 

VALU0*1 . /tCALf/100 . 

CALL  MMCRl  -0  20,0. 27,0.  U.VALM  .O.O.-l  I 
CALL  »M«IRl  0.40,0.27,0. 11,0.0,0.0,-11 
VALUI ■* ■ 0/3CALI/100 . 

CALL  **•«<  S.  771, 0.27,0.11, VALUi. O.O.-l  I 
VALUI *10 . 0/3CAL1/100 . 

CALL  NMill  10. 7,0. 27,0. 10 .VALUI ,0.0,-l ) 

CALL  SV)«0LIS.S.-O.2S,.lA,ANHinn.0.0.A) 

C 

CALL  SV790LI 10. 0,12- 0,1.1,42.0.0, -l  I 
CALL  SVHMLI  10.1.12.4,0.2,01.0.0,-11 
XOa-MO 
TO«-VO 

CALL  OLOTINO,VO.-S> 

C 

MTURM 

(NO 

C 

C 

aumounm  mum lamp,  lawn,  lamp,  loom,  loud*,  loom .imrtii.iaot i 

c 

OOUM.I  OOICXSXON  A,A,N,A0,00  .CO.OO.CO.O.Rl.COO.fOOO.RM.AO.KO 
OOUOLl  OMCISZON  R2.R3, A4,Ri.P,02.0S, 04,01, 04.A4.M.0XMRIC 
OOUOLI  OOICIOXON  ONI .OLAH.NOOTN.IAOT .O.ONXHXN.OI ,LON 
RIAL  LAMO.LOMO 

C 

C  THIS  OUMOUTXNI  C0H0UT3  TNI  imt  COOOOXNATIS  00  00  IN  NM  72 
C  IN  ZOm  10  CINTRAL  NIRZOZAN  121  00  00  M 

c 

cm-123.000 

ONZ  >OOLOAT  I LA0I8 1  *OOLOAT  I  LAMM  I/M .  00*0»LI  I LAMO  l/XAOO .  00 
LON-OOLOATI  LOOCO  1*0010071  LOKZN I/M.  DO*OOLH  L0080  1/5400-00 
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APPENDIX  I 

CONTROL  FILE  FOR  KALMN2  PROGRAM.  AS  DEFINED  BY 

BROW  N(  1986) 

585590 .  700000.  16  10  4  1  1  000 

0  -2740.506  -4341.136  3772.0820  .8738E1  0.0  -3.0  0.0 

6373.135  298.26 

1.0  1.0  1.0  l.E-12  2.1  l.E-14  l.E-50  l.E-50  l.E-50  3.0  4.0  5.0 
1  0  1  0  0  0  0  1  0  0 

0.000005  0.0 

1  1  0  0  0  0  1  0  1 

980.  15.  75.  10. 

0.75E-4  .75E-4  .75E-4  .IE-4  1.0  1.0  1.0  2.0  .IE-03  .IE-04 

1.0  1.0  1.0  1.0  1.0  1.0  l.E-19  l.E-19  l.E-19 

-.031 


APPENDIX  J 

CONTROL  FILE  FOR  CVFICA  PROGRAM. 


The  first  two  lines  may  have  a  title  up  to  80  columns.  The  third  line  has  the 
month  day  and  >ear  of  the  starting  day  of  the  GPS  week.  The  fourth  line  has  the 
number  of  satellites  to  be  used,  and  other  two  constants  are  kept  at  zero  The  last  hr.e 
has  the  estimate  position  in  WGS  ~2  coordinates 
The  next  is  the  one  used  for  this  thesis: 

SEAFLOOR  BENCH  MARK  PHASE  II 
R  V  POINT  SL  R  16  AlGLST  1986 

8  10  1986 
4  0.  0 

■2*39.58  -434077  37^3.14 


APPENDIX  K 

CONTROL  FILE  FOR  COMPARE  POSITION  AND  PLOT  PROGRAMS 


The  first  lines  has  the  number  of  periods  to  compare  the  data.  This  line  must  be 
followed  by  the  same  number  of  lines  as  periods  are  to  compare.  Each  line  the  date 
time  tags  from  the  start  and  of  the  periods  has  the  left  inferior  corner  of  the  area 
corresponding  to  each  period. 

The  next  is  the  one  used  for  this  thesis: 


II 


8  16  1986  19  14  0.0  8  16  1986 

8  16  1986  19  20  0.0  8  16  1986 

8  16  1986  19  25  0.0  8  16  1986 

8  16  1986  19  30  0.0  8  16  1986 

8  16  1986  19  35  0  0  8  16  1986 

S  16  1986  19  4"  0.0  8  16  1986 

8  16  1986  19  5*0.0  8  16  1986 
8  16  1986  20  0*0.0  8  16  1986 

8  16  1986  20  13  0  0  8  16  1986 

8  16  1986  20  20  0.0  8  16  1986 

8  16  1986  20  30  0  0  8  16  1986 


19  20  3.0  567320.0  4038400.0 
19  25  3.0  566920.0  4039000.0 
19  30  3.0  566670.0  4039600.0 
19  35  3.0  5662*0.0  4040200.0 
19  40  3.0  566020.0  4040800.0 

19  53  3.0  564920.0  4040950.0 

20  05  3.0  564970.0  4040200.0 
20  13  3.0  565570.0  4039500.0 
20  20  3.0  566020.0  4039000.0 
20  25  3.0  566420.0  4038500.0 
20  37  3.0  565985.0  4038400.0 
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